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ABSTRACT 



Jet Quenching in Relativistic Heavy Ion Collisions at 

the LHC 



Aaron Angerami 



Jet production in relativistic heavy ion collisions is studied using Pb+Pb collisions at a center of 
mass energy of 2.76 TeV per nucleon. The measurements reported here utilize data collected with 
the ATLAS detector at the LHC from the 2010 Pb ion run corresponding to a total integrated 
luminosity of 7 fib^ 1 . The results are obtained using fully reconstructed jets using the anti-fct 
algorithm with a per-event background subtraction procedure. A centrality-dependent modification 
of the dijet asymmetry distribution is observed, which indicates a higher rate of asymmetric dijet 
pairs in central collisions relative to peripheral and pp collisions. Simultaneously the dijet angular 
correlations show almost no centrality dependence. These results provide the first direct observation 
of jet quenching. Measurements of the single inclusive jet spectrum, measured with jet radius 
parameters R = 0.2,0.3,0.4 and 0.5, are also presented. The spectra are unfolded to correct for 
the finite energy resolution introduced by both detector effects and underlying event fluctuations. 
Single jet production, through the central-to-peripheral ratio Rqp, is found to be suppressed in 
central collisions by approximately a factor of two, nearly independent of the jet pr- The Rqp 
is found to have a small but significant increase with increasing R, which may relate directly to 
aspects of radiative energy loss. 
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Chapter 1 

Introduction 

The subject of this thesis is a series of experimental measurements of jets in relativistic heavy ion 
collisions which are related to the phenomenon of jet quenching. The high energy densities created 
in these collisions allow for the formation of a hot, evanescent medium of deconfined quarks and 
gluons termed the quark-gluon plasma (QGP). The QGP is a fundamental physical system described 
by quantum chromodynamics (QCD), in which degrees of freedom of the QCD Lagrangian are in 
principle manifest. The transition from ordinary nuclear matter to the QGP phase may provide 
critical insight into the dynamics of the early Universe. Experimental searches of such a medium 
in nucleus-nucleus collisions have been made at the AGS, SPS, RHIC and most recently in Pb+Pb 
collisions at the LHC. Measurements of collective phenomena such as elliptic flow have led to 
the interpretation that the produced medium rapidly equilibrates after which it is well described 
by near- ideal hydrodynamics. Highly collimated clusters of particles produced from hard partonic 
scatterings known as jets are produced in these collisions. Jets, which are not in thermal equilibrium 
with the rest of the plasma, have long been thought to be sensitive to detailed medium properties 
as they provide an external probe of the system. The process by which a quark or gluon loses 
energy and suffers a modification of its parton shower in a medium of high color charge density is 
referred to as jet quenching. 

The measurement of the suppression of high transverse momentum (pp) hadrons and the struc- 
ture of dihadron correlations at RHIC provided indirect experimental evidence for jet quenching 
and was one of the early successes of the RHIC program. These measurements indicated that 
the formulation of QCD factorization used in the calculation of hard scattering rates in hadronic 
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collisions could not be applied to nuclear collisions, implying that dynamics at new scales are an 
essential feature of the medium. However, the utility of these results was limited as single particle 
observables do not provide the full details of the jet-medium interaction. To relate partonic energy 
loss to single particle measurements a fragmentation function must be applied. This procedure 
loses sensitivity to the angular pattern of medium-induced radiation and forces the assumption 
of vacuum fragmentation. Furthermore, the limited rapidity coverage and low rate at RHIC for 
jets clearly above the underlying event fluctuations prevented more detailed conclusions from being 
drawn about the quenching mechanism. 

More recently, fully reconstructed jets, which carry the full information of the partonic energy 
loss, have been measured at the LHC. Dijet events were observed to become increasingly asymmetric 
in more central collisions, while retaining an angular correlation consistent with back-to-back jets. 
This result pQ has generated a great deal of excitement, as it strongly suggests the interpretation 
that back-to-back jets traverse different path lengths in the medium resulting in events with large 
dijet asymmetry where one jet leaves the medium mostly unscathed while the other experiences 
significant quenching. An event display of a highly asymmetric event is shown in Fig. This 
observation is the first direct evidence for jet quenching and is one of two results presented in this 
thesis. 

The asymmetry is sensitive to the energy loss of one parton relative to one another, and ad- 
ditional insight into the effects of quenching can be provided by observables that are sensitive to 
the quenching on a per jet basis. Energy loss of single partons could result in jets emerging from 
the medium with significantly less energy, resulting in a different final pp spectrum for events in 
which the medium is produced. Due to the steeply falling nature of the unquenched pt spectrum, a 
systematic downward shift in jet energies due to quenching causes a reduction in the total number 
of jets at a given pp. The second result presented in this thesis is a measurement inclusive pro- 
duction of single jets. Ratios of jet spectra in different centrality intervals relative to a peripheral 
reference were used to construct an observable, Rcp, quantifying the extent to which jet produc- 
tion is suppressed by the presence of the QGP through quenching. The Rqp has been measured 
as a function of jet and size as well as collision centrality, which should provide the ability to 
distinguish between different models of the quenching mechanism. 

This thesis is organized as follows. Chapter [2] provides background on both experimental and 
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Figure 1.1: A highly asymmetric dijet event recorded by ATLAS during the early portion of the 
2010 Pb+Pb run. A view along the beam axis is shown on the left, with high-pT charged particle 
tracks indicated by lines and energy deposits in the calorimeter by colored bars. The center figure 
shows the calorimeter distribution in rj — (f), indicating a highly energetic jet with no balancing 
jet opposite in azimuth. A similar distribution is shown (right) for the inner detector charged 
particle pt, which is consistent with the calorimeter signal. 

theoretical developments leading to the conclusion that QCD is the correct theory to describe the 
strong nuclear force. These ideas are developed to motivate the use of heavy ion collisions as an 
experimental tool, and various formulations of the jet quenching problem are presented. Chapter [3] 
describes the experimental apparatuses, the LHC and the ATLAS detector, used to provide the 
measurements presented in this thesis. Chapter [4] discusses the experimental techniques used to 
separate underlying event contributions from the jet signal and perform jet energy scale calibrations. 
Analysis of reconstructed jets is discussed in Chapter [5j including unfolding of the jet spectrum 
to account for detector effects and estimates of systematic uncertainties. The final results of this 
analysis, the jet Rqp and dijet asymmetry, are presented in Chapter [6] with conclusions following 
m Chapter [7| 
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Chapter 2 

Background 

Nuclear physics deals with the study of the fundamental properties of nuclear matter and the 
character of the strong nuclear force. By the 1960's the concept of nuclear matter had been extended 
from the nucleons making up atomic nuclei to a vast number of strongly interacting particles known 
as hadrons. The flavor symmetries among various types of hadrons suggested that these particles 
are composite particles with their quantum numbers carried by common, sub-hadronic objects, 
and that these objects possessed a new quantum number known as color [2]. A group theoretic 
approach to hadron taxonomy led to the development of the Quark Model [HHZ] which predicted 
the Q~ baryon observed at BNL in 1964 [5]. 

Developments in current algebra led to a proposed scaling behavior of the structure functions 
in deep inelastic scattering. Bjorken Scaling [9] implies that at high momentum transfer, Q 2 , the 
hadronic structure functions are independent of Q 2 and functions of a single scaling variable, x. 
Scaling behavior consistent with this relation was observed in a series of measurements in deep 
inelastic scattering (DIS) at SLAC [HUH], leading to the interpretation of point-like constituents 
within the nucleon known as the Parton Model |12H14j . 

Advances in renormalization group techniques |15H18j led to the conclusion that only a theory 
possessing asymptotic freedom, a coupling strength that becomes arbitrarily weak at high energies, 
could be consistent Bjorken Scaling \19\ I20j . The symmetry considerations of the Quark Model 
suggested that non-Abelian gauge theories [21] would be potential candidates for such a theory. 
These theories were found to possess the required renormalization criteria with the /3-function in the 
calculation of the running of the coupling being negative for an SU(3) gauge symmetry [22j[23]. Fur- 
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thermore, explicit calculations of the anomalous dimensions of the DIS structure functions showed 
that these non-Abelian gauge theories possessed scaling properties consistent with experimental 
observations 



This description of the strong interaction as an SU(3) (color) gauge symmetry between spin-- 
quarks and mediated by massless spin-1 gluons is known as Quantum Chromodynamics (QCD). At 
short distance scales, asymptotic freedom ensures that the coupling becomes weak and perturbative 
techniques can be used to calculate observables. The validity of QCD as the correct theory of the 
strong interaction has been rigorously demonstrated by a variety of experimental results that are 



well described by perturbative calculations. A summary of this evidence is presented in Section 2.1 
Color charges are never observed as isolated particles but rather as constituents of color neutral 
hadrons, an aspect of QCD known as confinement. Furthermore, the QCD Lagrangian possess 
chiral symmetry which is dynamically broken through the interactions that bind the quarks into 
hadrons and is responsible for the vast majority of the mass of ordinary matteiQ However, at soft 
momentum scales the coupling becomes increasingly large and at hadronic energies the theory is 
non-perturbative. This leads to significant challenges in theoretical calculations and interpreting 
experimental results, and much of the work done since the inception of QCD has been to navigate 
these challenges to develop a comprehensive picture of the strong interaction. 

2.1 Experimental Evidence for QCD 

Initial experimental evidence for QCD as a theory of the strong interaction was provided in DIS 
through the observation of Bjorken Scaling. Measurements of the proton structure function F2 are 



shown in Fig. 2.1 and DIS structure functions in the context of the Parton Model and collinear 
factorization are further discussed in Section |2.3| Subsequent tests of the theory have been per- 
formed in e + e~ collisions and hadronic collisions, particularly through the Drell-Yan process, the 
inclusive production of lepton pairs in hadronic scatterings, A + B — > l + £~ + X [28]. These results 
have provided overwhelming experimental support for QCD [29, 30], and a brief survey is discussed 
below. 



1 This symmetry is only approximate as it is broken explicitly by the quark masses. Thus the pions are pseudo- 
Goldstone bosons. 
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Figure 2.1: The proton structure function i 7 ^ vs Q 2 at fixed x measured by a variety of experiments 
(see [30] and references therein) . F2 is approximately independent of Q 2 . The observed logarithmic 
violations of this scaling are well-described by the QCD phenomenon of parton evolution. 



The Callan-Gross [3T] relation between the DIS structure functions, 2xF\ = F2, arises due to 
scattering off spin--, point-like objects; the measured behavior of the structure function, shown in 



Fig. 2.2, is approximately consistent with this relation. In e + e collisions the ratio of cross sections 
for the production of inclusive hadrons to muon pairs is given to leading order in QCD by 



R(s) 



o~(e + e — > hadrons, s) 
a(e + e~ — > fi + fi~,s) 



No 



f 



(2.1) 



where s is the center of mass energy, N c is the number of colors (N c = 3 for QCD), and Qf is the 
charge of quark /. The sum runs over the quark flavors up to the heaviest quark capable of being 
produced at that s. This ratio, aside from resonance peaks and away from the Z threshold, is 
measured to be approximately constant with discrete jumps at the quark mass thresholds as shown 



in Fig. 
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Figure 2.2: The ratio 2xF\/F2 as a function of x [32]. The ratio is approximately independent of x, 
consistent with the Callan-Gross relation, which is the result of scattering off of spin-- constituents 
within the proton. 



Highly collimated, energetic sprays of particles balanced in momentum known as jets were ob- 
served in e + e~ collisions at SLAC [33], consistent with the hypothesis that hadrons are produced 
in these collisions through the production of qq pairs of opposite momentum. The observation of 
three jet events at PETRA through the process e + e~ — > qqg provided experimental evidence for 
the gluon as well as measurements of the strong coupling constant as [SSj- Further precision mea- 
surements of the strong coupling constant have shown remarkable consistency [36] and a summary 



is shown in Fig. 2.4 



2.2 Fundamentals of QCD 

QCD is the quantum field theory that describes the strong nuclear force, and is expressed in terms 
of quarks and gluons which contain color charge described by an SU(3) gauge symmetry. The 
spin-- quarks fields, ip a , transform under the fundamental, three-dimensional, representation of 
the gauge group, 

4>' a (x) = e ia °^Mx) = U ab (x)Mx) , (2.2) 
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Figure 2.3: World data on inclusive e + e~ — > hadrons cross section (top) and ratio to muon pair 
production cross sections (bottom) [30, 



with the subscript a describing the quark color. The t^ b are the generators of the gauge group, 
represented by eight 3x3 matrices (C = 1, • • • 8, a,b = 1, 2, 3), and obey the Lie Algebra^] 

[t A ,t B ]=f ABC t c , (2.3) 

with f ABC the structure constants of SU(3). The local gauge symmetry is imposed by replacing 
the derivatives appearing in the usual Dirac Lagrangian with the gauge-covariant derivative, D^^f,, 



(2.4) 



2 In this section the generator/adjoint indices are denoted by capital, Roman script, as a superscript with the 
summation over repeated indices is implied. The fundamental, matrix multiplication indices are denoted by lower 
case, Roman script appearing as subscripts. The Lorentz indices are indicated by Greek characters and the standard 
Einstein summation convention applies. 
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Figure 2.4: The strong coupling constant, as, as a function of Q 2 (left) and fixed at 



Tflr 



(right). Figure adapted from Ref. |36j . From these results, as was determined to have a world 
average of a s (m z ) = 0.1184 ± 0.0007. 



The color indices a and b on D are explicit here to indicate that the covariant derivative is Lie 
Algebra- valued (i.e. matrix- valued) . The covariant derivative introduces a non-Abelian coupling, 
of strength g, between the quarks and massless, spin-1 gauge fields A^ representing gluons. These 
transform under the adjoint, eight- dimensional, representation of the gauge group. Including the 
pure Yang-Mills action for the gauge fields, the full QCD Lagrangian is given by 

1 



(2.5) 



with the index q denoting the quark flavor with mass m q . The field strength tensor F^ u is expressed 



in terms of the gauge field via 



F% = d ll A2-d v A$-gf 



ABC A B aC 



(2.6) 



The QCD vertices are shown in Fig. 2.5 In addition to the analog of the QED gauge coupling 
shown on the left, the non-Abelian structure of the theory allows for coupling between the gauge 
fields themselves. Color factors, associated with different processes, arise from the group structure 
in the Feynman rules, 



j.B _ rp x 

t cd t cd - +F0 



AB 



ACD rBCD 



CaS B , tfctcd = CpSftd- 



(2.7) 
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Figure 2.5: Possible couplings allowed by the QCD Lagrangian. The quark-gluon (left) and three- 
gluon (center) vertices are proportional to g. The four-gluon vertex is proportional to g 2 . 

In QCD, these have values 

T F = \, C A = 3, C F = A/3, (2.8) 

and are typically associated with the splittings g — > qq, g —> gg and q —> qg respectively. 

The quantization of this theory is typically performed using the path integral formalism. In this 
technique, the gauge-fixing condition imposes a constraint on the functional integral which can be 
removed by the introduction of Faddeev- Popov ghost fields [37] . These unphysical fields appear as 
an artificial mechanism for preserving gauge invariance. Their contributions, and the contribution 
of the gluon propagator, will depend on the choice of gauge, with this dependence dropping out in 
the computation of any gauge invariant quantity. 

2.2.1 Ultraviolet Behavior 

The ultraviolet (UV) behavior of QCD was the main reason it was originally proposed as the 
theory of the strong interaction. The relationship among the UV divergences is encoded in the 
Callan-Symanzik equation. In particular, the derivative of the coupling constant with respect to 
the renormalization scale, //, is defined by the /3-function, 

jjp = M ■ (2-9) 
a In fx 

The one-loop /3-function for SU(iV) non-Abelian gauge theories was first computed by Wilczek, 
Gross and Politzer in 1973 [22l [23]. For N = 3 it is given by 

/3(«s) = - 5 (g/3 1 + (g) 2 /3 2 + •••), (2.10) 
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with as = g 2 /4:TT. The one and two- loop coefficients are 



2 38 
fa = 11 --N f , h = 102 - — N f . (2.11) 



To one- loop, the solution is 



«B(A) = g ^ = , w t A2 , (2-12) 

l + ^as(ti)HWo) ^ ln ^/ A QCD) 

Aqcd = f, oe -^/(^M) , 

where Aqcd is used to set the scale of the strong coupling. To higher order in ln^ 2 /A this is 

„ (,\- ^ /?2ln[ln(/, 2 /A 2 )] / 1 \ 

asW A ln(^/A 2 ) /^lnV/A 2 ) +U V l /3 1 3 ln 2 (^ 2 /A 2 )y' ' 1 J 

where the definition of A becomes dependent on the renormalization scheme. The result above is for 

the MS scheme [29]. As the renormalization scale is increased the coupling decreases, as indicated 

by the leading — sign in the f3- function. This property is asymptotic freedom. For general SU(iV), 

Pi = -(11N C — 2Nf). The leading, positive term is due to the gluons, which reduce the coupling at 
3 

large \x. There is competition, which weakens the asymptotic freedom, from the second term which 
is due to fermion loops and is proportional to the number of flavors. 

2.2.2 Infrared Behavior 

Although QCD is well behaved at large momenta, the infrared behavior must be handled with 
care in perturbative calculations. Since gluons (and to a good approximation the light quarks) 
are massless, any sensitivity to the long range behavior of QCD appears in perturbation theory 
as an infrared divergence. Fortunately, the configurations that introduce this sensitivity have 
been systematically analyzed and a formal procedure exists to define quantities which are infrared 
safe [38 - 40j . These divergences are associated with the contributions to the momentum integration 
where the massless particle has zero momentum (soft) or at zero angle (collinear). Observables that 
are safe from these divergences must be insensitive with respect to the emission of an additional soft 
or collinear gluon. Formally, this means that for an inclusive quantity Z, defined by the functions 
<Sn> 

X= \\ I dn2 §V 2 S2iP ^ P * ) + h /^3^S 3 (^M + ---, (2-14) 
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the quantity is said to be infrared and collinear (IRC) safe if 

WpM."- - ,(l-AK,A^ +1 )=5 n (^,^,-- - X). (2-15) 
This idea is expanded in the discussion of applications of perturbative QCD in the following section. 

2.2.3 Non-Perturbative Dynamics 

Inspired by Wilson's picture of renormalization, a discrete lattice formulation of QCD was devel- 
oped to supplement perturbation theory where the latter is inapplicable [H]. This procedure uses 
the path integral formalism in Euclidean space-time, where the field configurations are explicitly 
integrated over. The integration is accomplished by recasting the quantity of interest in terms 
of gauge links between adjacent lattice sites. The finite lattice spacing, a, serves as an ultravi- 
olet cutoff, and a valid lattice formulation must respect all of the features of QCD in the limit 
a —7- 0. Technical challenges arise as spurious fermion states can be produced in this limit as lattice 
artifacts 143] . This can be mitigated, although not completely removed, while maintaining chi- 
ral symmetry, by introducing staggered fermions [54J. Furthermore, the finite spacing introduces 
integration errors in the evaluation of the action. These can be minimized by using improved defi- 
nitions of the action [45J . Advances in computing power and improved actions have led to highly 
accurate calculations, such as decay constants, form factors and the spectrum of hadrons, which is 



shown in Fig. 2.6 Additionally, the Euclidean path integral formulation is amenable to calculating 



thermodynamic quantities in QCD, a topic discussed in Section [2.4.2 

A qualitative picture of confinement is to view the static qq potential as Coulombic at short 
distances, but growing at long range, 

4 a 

V QC D(r) = --- + kr. (2.16) 



A recent lattice calculation of the static quark potential is shown in Fig. 2.7 which indicates 
behavior approximately consistent with this analytic form. As the separation between the qq pair 
increases the force remains constant, and an increasing amount of energy is stored in the stretched 
gluon field called a flux tube. At some separation it becomes energetically favorable to break the 
flux tube and produce a new qq pair from the vacuum. This process, known as string fragmentation, 



is shown schematically in Fig. 2.8 Such a mechanism provides a heuristic construct for modeling 
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Figure 2.6: Lattice calculation of the hadron spectrum with comparison to measured values (left) 
and ratios (right) 




12 3 

R/r, 

Figure 2.7: Static quark potential as a function of r calculated with different values lattice spacings 
|46j . The distance is expressed in units of the lattice size, with the colored gradations indicating 
the scale in physical units. 
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confinement; no matter how much energy is applied only color- neutral objects can be created from 
the system. The flux tube picture has led to developments in hadron phenomenology by considering 



hadrons as relativistic strings |47] 148] . and the constant, k, in Eq. 2.16 can be interpreted as a string 
tension. Monte Carlo (MC) event generators for simulating the production of hadrons have had a 
long history of success using string models |49] 150] . 



Figure 2.8: Flux tube representing static qq potential (left) fragmenting into two color-neutral 
objects (right). 



2.3 Applications of Perturbative QCD 
2.3.1 The Parton Model and Factorization 

QCD has had great success in providing reliable perturbative calculations of experimental observ- 
ables. The most basic application of QCD uses the parton model, which is a tree-level, impulse 
approximation to the full perturbation theory. The scattering process is formulated in terms of 
point-like constituents integrated over probability distributions for a given parton to have momen- 
tum fraction x, 

&ab(pa,Pb) ~ / dx i dx j&( x iPA,XjPB)(f>t{Xi)(t>f{x j ), (2.17) 

Here A and B denote the colliding hadrons and i and j denoting partons of a particular type. The 
functions <j>f(xi) give the probability density of finding a parton of type i in hadron A with momen- 
tum pi = XiPA- While asymptotic freedom ensures that at sufficiently hard scales the partonic level 
matrix elements are calculable in perturbation theory, almost any real world observable will involve 
hadronic initial and/or final-states, where the theory becomes non-perturbative. The applicability 



of Eq. |2.17 in certain kinematic regimes suggests that calculations can be performed by separating 
the short distance behavior of QCD, encoded in the partonic matrix element a, from the long range 
behavior represented by the probability distributions. The formal apparatus through which this 
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is accomplished is known as factorization. Factorization theorems take a similar form to Eq. 2.17 



where each term in the integrand is also dependent on a scheme-dependent factorization scale, [if. 
This parameter sets the cutoff for which aspects of the dynamics are being included in the descrip- 
tion of the long and short range components. The form of the factorization ensures that there is 
no quantum mechanical interference between the long and short range behavior. Furthermore, the 
probability distributions, which are known as parton distribution functions (PDFs) in the general 
factorized formulation, are independent of the specific scattering process. Although they are not 
calculable by perturbation theory, once these functions are determined experimentally they can 
be applied to a calculation regardless of the details of the partonic level scattering, thus making 
them universal features of the hadrons. Factorization theorems have been proven for a variety of 
processes, A + B — > C + D among particles A, B, C and D, with the schematic form of the theorem 
being 

da(A + B -)• C + D) = da ® $a ® $B <8> Ac <g> A D + P-s.c. (2.18) 



Here $ and A denote the parton distribution and fragmentation functions (see Section 2.3.4) and 
<S> represents a convolution over parton momentum fraction, trace over color indices and sum over 
parton species. This form is traditionally referred to as a twist expansion, with the first term known 
as the leading twist. The remainder terms, the power suppressed corrections, are suppressed by 
powers of a hard scale and are referred to as higher-twist effects. Most factorization theorems have 
been proven using a collinear factorization scheme, where the transverse momentum of the parton 
has been integrated over. These theorems have been proven in the cases of DIS and e + e~ , but 
the only case involving hadron-hadron scattering to have been rigorously proven is the Drell-Yan 
process [5TH53] . 

Observables such as spin asymmetries in hadronic collisions provide access to the spin structure 
of the proton. However, the extension of the usual collinear formalism in these cases has been 
problematic [54-57J as the observables contain sensitivity to the transverse momentum dependence 
(TMD) of the PDFs. Without integrating over the transverse dependence of the PDFs, gauge links 
couple soft and collinear portions of diagrams to the hard sub-graphs, breaking the factorization. 
Attempts to construct TMD factorization schemes, such as /c t -factorization |58} I59|. have proven 
difficult and it has been shown explicitly that fc t -factorization is violated in high-px hadron pro- 
duction in hadronic collisions |60j. Building a framework to perform these types of calculations is 
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still a subject under active development [61] . 
2.3.2 Deep Inelastic Scattering 

The parton model was first developed in DIS, and this collision system will be used here to illustrate 
the utility of the factorized approach. In these collisions an electron interacts with a hadronic target 
of mass M via the exchange of a photon with large virtuality q 2 = —Q 2 . In addition to Q 2 , the 
system is typically described in terms of the variables x and y 

Q 2 E — E' , 

y = ^^, (2-19) 



2g-P' " E 

where P is the four-momentum of the target and E and E' are the initial and final energies of the 
electron as measured in the rest frame of the target. In the parton model, x corresponds to the 
fraction of target's momentum carried by the struck quark. The differential cross section is usually 
expressed in terms of structure functions F, 

da _ 4ttq;| m f 2 2 / 2Mxy\ ^ ^ 



ixiy zyQl { V^,Q') + [1- V -^)F 2M ^^ (2.20) 

The observation of Bjorken scaling [TUJ, [TT] , that for large Q 2 , F\ and F2 are independent of Q 2 , 
was the inspiration for the parton model. For spin-- (charged) partons, the structure functions 

obey the Callan-Gross relation 2xF\ = F2. The experimental confirmation of this result provided 

.1 

support for the interpretation of spin-- quarks. In collinear factorized form the structure functions 
have the form, 

F 1 (x, Q 2 ) = W 1 f f a /A& V 2 ,Vf)Ci a ^,a(fi)) + 0{Q- 2 ) , (2.21a) 

F 2 (x,Q 2 ) = Y^ f dtf a/A (Z,i?,iMf)Cka(j,%,%a(ri} +0{Q- 2 ). (2.21b) 

a Jx \Z V M / 

Here the functions f a /A are t ne PDFs for parton a in hadron A. The coefficient functions, C a , 
are calculable in perturbation theory and are IRC safe. They are dominated by contributions of 
order Q; propagators off shell by fi 2 will contribute to C a while contributions below this scale are 
grouped into f a /A- The result depends on both the factorization and renormalization scales, which 
need not be set equal. 
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Shortly after the computation of the /3-function, similar renormalization group techniques were 
applied to the DIS structure functions, by considering moments of the PDFs (specifically non-singlet 

PDF, f q — fg), 

(2.22) 



Jo 



These are related to the anomalous dimensions, 

d 

In 



-/ (n V), 



(2.23) 



d In fx " 

which appear to give dimensions to dimensionless, but scale-dependent quantities through renormal- 
ization (e.g. ln/( n ) ~ /u~ 7n ), and can be computed in perturbation theory [20j. The Q-dependence 
of the PDF moments can be used to determine the Q-dependence of the moments of the structure 
functions, 



J x^F^Q^dx = c{ n) (a(Q 2 ))f( n \Q 2 ) x exp 



dt ln {a{Qle 1 )) 



(2.24) 



This relation predicts power-law scaling, Fi(x,Q 2 ) ~ {Q/Qo)~ a °"'"^ 7r , for theories that do not 
possess asymptotic freedom and posed a problem reconciling potential theories with the observed 
scaling prior to the demonstration of asymptotic freedom in QCD. To one-loop in QCD the scaling 
is [21H27] 

"lnQ 2 /A 2 ' 



Fi(x,Q z ) oc 



-27n/4|/3i| 



(2.25) 



\nQl/K\ 

This weak Q-dependence correctly described the early measurements of scaling in DIS, further 
building support for QCD as the correct theory of the strong interaction. 



2.3.3 DGLAP and Parton Evolution 

When calculating C a , logarithms can arise due to the collinear emission at some scale A, which 
gives a factor of approximately 



a{Q z )\n 



A 2 ' 



(2.26) 



These logarithms can become large such that the above product becomes of order unity. The result 



is that C a will have large contributions from all orders and is shown schematically in Fig 2.9 This 



Q 2 

multiple emission is enhanced by a(Q 2 ) N \n N but only for the case where the virtualities are 



strongly ordered, q 2 , < q\_ 



< 



< 



A 2 ' 

with other orderings appearing with fewer powers of In 



A 2 ' 
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Figure 2.9: Schematic representation of multiple collinear splittings. The portion of the diagram 
shown in red is interpreted as part of the structure of the struck hadron instead of a correction to 
the hard matrix element. 



The large logarithms are a symptom of interactions far away from the scale at which the coupling 
was fixed. Fortunately, these collinear contributions can be resummed by renormalization group 
methods. As most of the collinear emissions are well separated in scale from the probe, qo, these 
emissions can be reinterpreted as modifying the hadron structure as opposed to corrections to C a . 
This results in a Q-dependence of the PDF that evolves the probe from the hard scale to lower 



momentum scales (indicated by the red sub-diagram in Fig 2.9). For the change Q — > Q + AQ the 



differential probability of an emission with energy fraction z and transverse momentum Q < p± < 
Q + AQ is given by 

" ' ,J± -P a ^ b (z) ~ -^P a ^ b (z), (2.27) 



<* d P± P l (A ^ aA Q 

27T p 2 ± ^ b[Z ' - 7T Q 



where P ai - b {z) is the splitting function for parton of b splitting into type a, and can be computed 



from the diagrams shown in Fig |2.10 Changes in the distribution of parton a at momentum fraction 
x can come from splittings of other partons at x' = x/z, and can be written as 

Af a (x,Q) = V / dx' ! dz a ^P a ^ b (z)f b (x',Q)5(x-zx') (2.28) 
^ Jo Jo n Q 

V n Jx Z Z 
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Figure 2.10: Collinear QCD processes used to compute splitting functions. 



a 



i 7T Jx z z 



dlnQ 



(2.29) 



See Ref. [30] for a complete tabulation of the splitting functions and full evolution equations. This 
result indicates that if f a is measured as a function of x at a given value of Qo, the PDF at any 



other scale Q can be determined by evolving the PDF according to Eq. 2.29 The Q dependence 



of the structure function Fi at fixed x is shown in Fig. 2.1 This dependence agrees well with the 
dependence predicted by the DGLAP equations over a range of Q 2 and x values. 



2.3.4 Fragmentation Functions 

The final-state analog of the parton distribution function is known as a fragmentation function. The 
fragmentation function Df(z, /x 2 ) encapsulates the probability that a parton of type i will fragment 
into a hadron of type h with momentum fraction z of the original parton. In e + e~ collisions, where 
there are no initial parton distributions, cross sections for the production of a particular species of 
hadron can be written as, 



da h ' 



dz 



C t (z,a s ^ 2 ),^)D^z,^) + 



1 



(2.30) 



dx j x z 

For situations with arbitrary numbers initial/final-state hadrons, the collinear factorized relation 



takes the form of Eq. 2.18 The e + e measurements provide the cleanest access to the fragmentation 



functions since there is no integration over PDFs. The fragmentation functions obey an evolution 
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equation in /i 2 identical to Eq. 2.29, with the inverse splitting functions (i.e. P a <—b(z) is replaced 
with P a ^{z)). In the parton model these obey the momentum sum rule, 

J2 f 'dzzD?(z,p?) = l. (2.31) 

The higher-order corrections to the splitting functions can see logarithmic enhancements at low 
x, an effect which causes the leading order approximation of the evolution to break down much 
sooner in e + e~ than DIS. This has led to development of calculations using a modified leading log 
approximation (MLLA), which include corrections for next-to- leading effects [B5-69J. In general, 
the behavior of these enhancements is to cause the "hump-backed plateau" behavior, 

x J D(x, S )ocexp|-^(e-^) 2 |, (2.32) 

with £ = lnl/x and the width and peak approximately given by [70J, 

e p ~^lns/A 2 , croc (lns/A 2 ) 3/4 . (2.33) 



This behavior is shown in Fig. 2.11, with the fragmentation functions fit with a Gaussian for a 
variety of y^s and Q 2 values shown on the right and a comparison of the extracted peak position 
to the MLLA value shown on the left. 

To span a large kinematic range and provide separation between quark and anti-quark contri- 
butions, the e + e~ data has been supplemented by DIS and hadronic collider data in recent global 



NLO extractions |71H73] (see 2.3.5 for discussion of fixed order calculations). A comparison of 



results of some such analyses are shown in Fig 2.12 where the various partonic contributions to 



the charged pion fragmentation function are shown. 



2.3.5 Fixed Order Calculations 



Perturbative calculations are typically performed at fixed order and it is useful to consider IRC 
safety in this context |30j . For an n-particle process, the cross section for an observable, O n , can 
be constructed at leading order (LO) as 

°o,lo = as~ 2 (w0 / dn n\Ml \(pi, ■ ■ ■ ,p n )O n ( Pl , ■ ■ ■ ,p n ), (2.34) 

where M 2 is the n-particle matrix element at tree-level, d£l n is the measure over the n-particle 
phase space and /xr is the renormalization scale. Non-perturbative corrections have been omitted 
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Figure 2.11: The hump-backed plateau structure of fragmentation functions through the measure- 
ment of single inclusive hadron cross sections (right). The cross sections are fit with a Gaussian 
and the extracted peak position compares well with the MLLA QCD prediction. Figure adapted 
from |30j. 



here. The tree-level amplitudes are finite but the integral will diverge in the soft and collinear 
regions of the momentum integration. Therefore this calculation is only safe when O n also vanishes 
in this limit (e.g. jet cross sections). At next-to-leading-order (NLO) the n + 1-particle matrix 
element must be considered along with the interference between the tree-level and 1-loop n-particle 
processes, 

O"0,NLO = <to,lo (2.35) 

+ «S _1 (W0 J d ^n+ll M n+l,ol(Pl) • • • ,Pn+l)O n+ i(pi, ■ ■ ■ ,P n +l) 

+ a™-\fi R ) J dn n 2Re [M nfi M* yl ] (p X) • • • ,p n )O n ( Pl , ■ ■ ■ ,p n ). 

Here M n \ is the 1-loop n-particle amplitude, which diverges in the soft and collinear limits. This 
divergence will cancel with the aforementioned tree-level divergence in the integration if the ob- 
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servable is IRC safe, 

lim O n+ i(pi, • • • ,pi, ■ ■ ■ ,p n +i) = O n (pi, ■■■ ,p n ) (2.36) 

lim O n+ x(pi, ■ ■ ■ ,Pi,Pj, ■ ■ ■ ,p n +i) = O n {pi, ■ ■ ■ ,pi, ■ ■ ■ ,p n ) . 
Pj\\pi 

A similar procedure applies in extending to NNLO calculations, which involve additional contri- 
butions the tree-level n + 2 contribution, interference between tree-level and 1-loop n + 1-particle 
diagrams, interference between the n-particle tree-level and 2-loop diagrams and the squared con- 
tribution of the 1-loop n-particle diagram. These calculations become increasingly difficult to 
compute, especially in an automated fashion, due to the intricate higher-loop diagrams. 

Often these calculations are supplied in the form of Monte Carlo event generators. These tools 
typically compute a hard scattering matrix element and then evolve the partons by means of a 
parton shower. This technique constructs a probability for a parton not split from the splitting 



functions defined in Eq. 2.27 The probability that a parton of type a will not split between scales 
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Q 2 and Qq, A a (Q 2 , Qq), known as a Sudakov form factor is given by, 

A a (Q 2 ,Q 2 ) = expl- ^ f dzf-^P^z)]. (2.37) 

MC sampling determines the scale of the first splitting, and this process is repeated down to some 
hadronization scale Q ~ 1 GeV. At this point a hadronization model |49} [71} 175] converts the 
evolved partons into final-state hadrons which can be used to simulate response of a detector to 
an input physics signal. Leading order event generators (e.g. PYTHIA) typically give the correct 
description of the soft and collinear emission, but fail to describe wide angle emissions due to the 
lack of multi-parton matrix elements. This can be improved by generators which use this higher 
order matrix elements and merge the parton showers in such a way as to avoid double counting the 
real and virtual contributions |76H79j . 

In general LO calculations for applications at hadron colliders are accurate to approximately a 
factor of two. A ET-factor, which represents the ratio of NLO to LO contributions, is typically used 
to tune the cross sections. 



2.3.6 Jets 

Highly collimated sprays of particles known as jets are a ubiquitous feature of high energy particle 
collisions. These objects are the experimental signature of partonic hard scattering processes and 
have been used to test the theoretical application of the theory since their initial observation at 
SLAC [M! • However, the multi-particle final states associated with jet production are sufficiently 
complicated that any experimental measurement must have a precise definition of an observable [38J . 
Furthermore, these definitions must be relatable to quantities that are well-defined theoretically, 
particularly with regard to the IRC safety issues discussed previously. 

Most early jet algorithms were based on associating particles or energy deposits that were nearby 
in angle, typically referred to as cone algorithms. The Snowmass workshop in 1990 [80] proposed 
a standard set of criteria for jet algorithms [81J, however implementations of these criteria were 
not immediately obvious. Many cone algorithms used at the time suffered from varying degrees of 
IRC unsafety, which limited their relatability to theory [82]. The effects of unsafety are illustrated 



in Fig 2.13 where two algorithms, one safe (left) and one unsafe (right) are compared. In the 



IRC safe algorithm, the emission of a collinear gluon does not affect the clustering (b), allowing the 
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divergence in this graph to cancel the 1-loop (a) divergence. However, this emission causes different 
clustering in the unsafe algorithm (d), resulting in a failure in cancellation of the divergences. In 
general, issues of split or overlapping jets are not handled in a natural way by cone algorithms, and 
some additional criteria, often introducing undesired artifacts, are required to handle such cases. 



Collinear safe jet alg. 

a) b) 



jet 1 jet 1 

O^X (-00 ) o^X (+oo ) 

Infinities cancel 



Collinear unsafe jet alg 

c) d) 



jet 1 



jet 1 



jet 2 

O^X (-00 ) o^X (+oo ) 

Infinities do not cancel 



Figure 2.13: Jet reconstruction using two algorithms on a system with and without additional 
collinear gluon radiation. The IRC safe algorithm reconstructs the two scenarios as the same jet 
(a,b), allowing for cancellation of divergences. The unsafe algorithm gives different results in the 
two cases (c,d). Figure adapted from Ref [83J. 



Sequential clustering algorithms, which perform pairwise clustering of particles |84l I85j . showed 
particular promise as they were IRC safe. This led to the proposal for the kt algorithm in e + e~ 
collisions [86] , which was formulated to be boost- invariant and IRC safe and was later adapted for 
hadronic collisions |87j . The inclusive formulation of this algorithm [88] uses the measure, 



dij = min(?4iiPTj)-fl^-> AR ij = (w - Vjf + 



(2.38) 



where y and <f> are the particles' rapidity and azimuthal angle respectively, to cluster nearby parti- 
cles. The algorithm works using the following steps: 

1. Calculate dij for all pairs and also diB = v\i- 

2. Find the minimum of the dij and diB- 

3. If the minimum is a dij, combine i and j and start over. 

4. If the minimum is a diB, call i a final-state jet and remove it from subsequent clustering and 
start over. 
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5. Finish clustering when no particles remain. 

The parameter R controls the size of the jets in analogy with the radius parameter in cone jets. 
In this inclusive form, all particles are clustered into final-state jets; thus some final discriminating 
criteria, often a minimum cut on or pjj must be applied. This algorithm has the feature that for 
soft or collinear pairs the distance measure is inversely proportional to the differential probability 
of collinear emission, 

dPk ^i „ ( 2 39) 

dEidOij mm(E j ,E j )6 ij , v ' ; 

which is easily adaptable to theoretical calculations |82j . This algorithm was not originally favored 
by experiment due to its slow processing time and geometrically irregular jets, which made experi- 
mental corrections more difficult. However, a faster implementation by the Fast Jet package(0(iV In N) 
vs. 0(N 3 ) ) has significantly improved the former issue [89j. 

The distance measure of the kt algorithm can be generalized in the following way: 

AR 2 

^•=min(^.4^_Jl. (2.40) 

The choice p = 1 corresponds to the usual kt algorithm, while p = corresponds to the Cam- 
bridge/Aachen algorithm |90| 191] . which uses only geometric considerations when clustering ob- 
jects. The generalization was used to define the anti-fct algorithm (p = —1), which clusters hard 
particles first and produces jets similar with regular, cone-like geometry [83]. Examples of the 



clustering behavior of the k t , ATLAS cone [92] and anti-&t algorithms are shown Fig. 2.14. The kt 
allows jets clustered from soft background particles to compete with the real jet signal, a feature 
which is particularly problematic in heavy ion collisions. This is contrasted by the behavior of the 
anti-/c t which preferentially clusters particles with the harder of two jets. 

Jets have become an essential part of high energy experiments, both as tools for testing QCD, 



but also as input into reconstructing more complicated physics objects. Figure 2.15 shows the 
data/theory ratio for the inclusive jet cross section as a function of jet for a wide range of 
experiments. The theoretical calculation is provided by f astNLO with NLDJET++ 

2.4 Phase Structure of Nuclear Matter 



Prior to the advent of QCD, attempts were made to apply statistical methods to the large particle 
multiplicities in high energy collisions \94\ 195] . This led to the analysis of the high temperature 
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Figure 2.14: Jet clustering on an event in rj — (j) using three different jet algorithms: k t (left), 
ATLAS Cone (center), and anti-/c t (right) all with R = 0.4. The cone algorithm does a poor job of 
resolving the split jet. The kt and anti-fct resolve this as two separate jets, but the kt clustering is 
strongly influenced by other soft particles in the event. 

limit of nuclear matter in the context of hadronic resonances. The growth in the number of these 
resonance states with increasing energy led Hagedorn to propose the Statistical Bootstrap Model, 
where resonances are thought of as being composed of excited lower mass resonance constituents. 
In the high temperature limit the density of resonant states and the thermodynamic energy density 
of states approach each other and the partition function diverges. The temperature cannot be 
increased beyond some limiting value, the Hagedorn temperature, as adding energy will only excite 
more resonances, not add to the kinetic energy of the system [96J. 

Asymptotic freedom indicates that this will not occur, but rather that at sufficiently high 
temperatures the coupling will become weak leading to nuclear matter with dynamics described by 
perturbative QCD (97] . This suggests that there is a phase transition where the degrees of freedom 
evolve from hadrons to quarks and gluons. This state of matter the Quark Gluon Plasma (QGP), 
represents the form of nuclear matter in the early stages of the Universe [98J . 

In the weak coupling limit the Stefan-Boltzmann relation between energy density, e, and pres- 
sure, p, applies: 

e - 3p = 0. (2.41) 
Each is proportional to T 4 with the constant of proportionality indicating the degeneracy of particles 
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Figure 2.15: Ratio of data to theory for the single inclusive jet cross section at a variety of ener- 
gies [301 [93]. 



obeying Bose-Einstein, <7be, and Fermi-Dirac, ^fd, statistics, 



7 



9be + o^fd 



7T" 



(2.42) 



8 J J 30 

At the lowest temperatures the system is hadron gas containing only the lowest state, the pions 
(^be = 3i S ospin) <?fd = 0). At high temperature the system is a gas of deconfined quarks and gluons 
(<?be = 8 co ior x 2 spin , <7fd = 3 co i or x 2 spin x 2 q g x Nf). The total energy density is 



30 



3 T ~ 

16+yiV/ T^OO 



(2.43) 



For two flavor QCD the difference in phases corresponds to a factor of 9 increase in e/T 4 , indicating 
a dramatic change in energy density. 
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The temperature range over which this transition occurs, specifically the possible divergences 
of thermodynamic variables, can characterize the transition. In the limit of zero quark mass, the 
QCD Lagrangian possesses chiral symmetry, which is broken by non-perturbative effects resulting 
in three massless pions. The high temperature limit exhibits both deconfinement and restoration of 
chiral symmetry, and it is unclear to what extent these two phenomena are interrelated and whether 
QCD may exhibit a distinct phase transition associated with each of these phenomena [99]. 

2.4.1 Thermal Field Theory 

The extension to finite-temperature field theory is accomplished by using the relationship between 
the path integral formulation in quantum field theory and the partition function of statistical 
mechanics (100|, UOlj . Explicitly, the path integral for a single quantum mechanical particle with 
Hamiltonian H is given by 

/•?'(*) 

(q'(t)\ e~ iHt \q(0)) = / 2%(i)K So dt ' L , (2.44) 

which denotes the amplitude for a particle at q at t = to go to q' at t, and L = -mq(t) 2 — V(q(t)). 
This functional integral can be analytically continued to imaginary time via t —> —it. This generally 
provides a more rigorous definition for the path integral and its convergence. In the imaginary time 
formulation the action in the exponent is replaced by the Euclidean action 

S E = f dr'L E {r'), (2.45) 
J o 

where L E = -niq(t) 2 + V(q(t)). This can be compared directly to the partition function of the 
canonical ensemble, with temperature T and /3 = 1/kT, Z((3) = Tr[e - ^]. In the non-diagonal 
basis \q) this can be written in the same form as the usual path integral, 

Z(fi) = f dq(q\e- p6 \q) = [ 2%(r)]e~ So dt ' L z (2.46) 

J Jq(0)= q (r) 

Thermal expectation values of time-ordered products of fields (• • • )g can be computed in the same 
fashion as zero-temperature field theory, but using the imaginary time/Euclidean action. For 
example, a generating functional, Z^j3,J), can be constructed by introducing a source J, and 
adding the term Jq to the Lagrangian. A two-point correlation function can be evaluated with 



CHAPTER 2. BACKGROUND 



29 



functional derivatives with respect to J , 

The cyclic nature of the trace imposes the constraint A(r — j3) = ±A(r), with the sign depending 
on whether the field obeys commutation or anti-commutation relations. Thus the periodic (anti- 
periodic) boundary conditions on the correlation functions enforce the Bose-Einstein (Fermi-Dirac) 
statistics and A can be described as a Fourier series 

oo 

A(r) = p Yl e~ iuJnT A(ioj n ). (2.48) 

n=0 

Here 0J n = —— for bosons and w n = — — — - for fermions and are known as the Matsubara 

P P 
frequencies. The full real-time gauge theory propagators at tree level are 

D F (k ) = k 2_^ + u + Wti - ^ 2 ) e ^ r y, (2.49a) 

S(p) = j^ + 2vnV -m 2 Kp + m)^^, (2.49b) 

for the bosons and fermions respectively. In both expressions, the first term describes the usual 
vacuum propagation while the second represents the disconnected part of the propagation which 
counts the modes absorbed by the medium. 

2.4.2 Lattice Thermodynamics 

The lattice formulation of QCD discussed previously can be naturally applied to calculate thermo- 
dynamic variables in QCD, since it is formulated as a Euclidean space-time evaluated path integral. 
The four-dimensional Euclidean space-time is represented by x N T lattice sites with spacing a. 
The system is described by volume V and temperature, T, 

V = (N fT af, T=J-, (2.50) 

and the thermodynamic limit is taken with N T — > oo, a — > but holding T fixed. Statistical errors 
can enter since the integration is performed over a finite number of gauge field configurations, 
thus making calculations with large iV r difficult. Recent developments in improved actions with 
reduced discretization errors allow for accurate calculations without resorting to large iVV, and have 
facilitated advances in lattice thermodynamic calculations. 
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Typically the action is separated into a contribution from the gauge fields and a fermionic con- 
tribution. The latter is expressed as a functional determinant in a way that removes or minimizes 
fermion doubling problems. The action and determinant are written in terms of link variables, 
which involve the evaluation of gauge field configurations at adjacent lattice sites. The numer- 
ical integration techniques require a probabilistic interpretation, and problems occur when the 
determinant becomes negative. For zero chemical potential, fj,, this is not a problem, however 
the determinant becomes complex- valued for ju 7^ 0. This makes calculations at finite chemical 
potential more difficult, although recent advances have been made on this subject. 

The QCD equation of state can be determined on the lattice by evaluating the trace anomaly 
(the trace of the stress-energy tensor Q^ U (T)) in the grand canonical ensemble, 

ew _ e - 3 P _ a f P \ 

T 4 - t a ~ qt VT 4 / ' 1 j 

Since only expectation values can be calculated, the variables must be re-expressed in this form. 
The pressure can be evaluated directly from integrating the trace anomaly. The entropy density, 
s, and speed of sound c s , can likewise be calculated from the relations 

j__e + p *_dp 

T 3 " T 4 > C s " d£ ■ ^ bZ > 



2.16 



The pressure and energy density divided by T 4 from a recent calculation |102j are shown in Fig. 

The lattice results indicate a rapid, but continuous transition in e at Tq ~ 170 MeV, which is not 

indicative of a first or second order phase transition, but rather a smooth crossover between phases. 

Additionally, order parameters describing the transition to deconfinement and the restoration of 

chiral symmetry have been formulated on the lattice. The deconfinement order parameter is the 

expectation value of the Polyakov loop, (L) ~ exp(—F q (T)/T) where F q is the static free quark 

d 

energy. The chiral transition is described by the chiral susceptibility v = — — — (qq) . These 

dm m=o 

quantities involve subtleties relating to their renormalization and corrections due to finite quark 



masses. A calculation of these two quantities [102] is shown in Fig. 2.17 The results indicate that 
both of these transitions occur smoothly over the same temperature region where the equation of 
state shows a rapid change in energy density. 
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Figure 2.16: Pressure and energy density divided by T 4 as a function of temperature in a lattice 
calculation using physical quark masses [102J. 
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Figure 2.17: Lattice calculations of the (L) (left) and x/T 2 (right) for light quarks as functions of 
temperature. These are order parameters for the deconfinement and chiral symmetry restoration 
transitions respectively. The vertical lines indicate the temperature range over which the rapid 
cross over is inferred to occur from the equation of state. Figure adapted from Ref. [102J. 
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2.4.3 The Hard Thermal Loop Approximation 

When using perturbation theory at finite temperature, including the effects of loops becomes prob- 
lematic as the finite scale introduced by T breaks the usual perturbative expansion. If the external 
momenta are hard, all loops will contribute to the same order as the tree level diagram. However if 
the external momenta are soft, P ~ gT, only the contribution from Hard Thermal Loops (HTLs) 
will compete with the tree level term. These are cases where the "hard" loop momenta are much 
greater than the "soft" external momenta and the dominant contribution to the loop comes from 
momenta of order T. In gauge theories, the HTLs can be identified by power counting and an 
effective theory of resummed HTLs can be derived without affecting the vacuum renormalizability 
properties of the theory [103^ I104| . These loops lead to the generation of thermal masses m and 
nif for the gauge bosons and fermions respectively. The static component of the gauge boson 
propagator takes the form 

D 00 = -2.^2 , (2-53) 
q z + 2m z 

where the thermal mass gives rise to Debye screening. In the static limit, uj/q — > 0, the magnetic 
behavior is similar to the above expression, but with 

2 irm 2 U) 

m l -)■ . (2.54) 

2 q 

pKUJ 

This can be interpreted as a frequency-dependent dynamical screening with cutoff m. — , which 

V 2 <? 

is in some cases effective in protecting infrared divergences [105] . The full thermal masses including 
the effects from a finite chemical potential, fi, are given by 

m 2 = \g 2 T 2 C A + ±-g 2 C F (T 2 + 4^ 2 ), (2.55a) 

O 12 7T Z 

m 2 = \ g 2 (T 2 + \v 2 )C F , (2.55b) 



with the color factors given previously in Eq. 2.8 



One issue that is still not fully understood is the screening of the chromomagnetic fields. Using 
power counting arguments it can be shown that a diagram containing I + 1 gluon loops, shown in 



Fig. 2.18, with £ > 3, has infrared divergences of the form [106] 
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Figure 2.18: Diagram containing £ + 1 gluon loops 

For the longitudinal components this is regulated by the electric mass and gives behavior ~ g i+3 T 4 , 
which differs from the g 2e behavior expected from perturbation theory. Based on arguments from 
theories that are infrared-equivalent to QCD, the magnetic mass is of order ~ g 2 T [107] . This 
leads to infrared divergences in the transverse part of the t + 1 gluon loop diagram to go as ~ <7 6 T 4 
and signals a total breakdown of perturbation theory. This remains an outstanding problem, with 
thermal field theory not on rigorous theoretical ground except at asymptotically high temperatures. 



2.4.4 Heavy Ion Collisions 

The phase structure predicted by the lattice indicates that at low baryon chemical potential and high 
temperature, hadronic matter undergoes a rapid, but continuous increase in the number of degrees 



of freedom. This transition is illustrated in a possible QCD phase diagram, shown in Fig. 2.19 In 
addition to this transition, abnormal forms of nuclear matter can exist at high densities [108t 109J. 
Astrophysical objects such as neutron stars are stabilized against gravitational collapse by the 
degeneracy pressure created by an extremely high density of neutrons. At the highest densities 
novel states such as quark liquids and color superconductors (color flavor locked) are thought to 
exist pl0|. 

It is believed that the QGP transition is experimentally accessible in collisions among heavy ions 
at high energies. In these collisions the nuclei appear as highly Lorentz contracted "pancakes" [111] . 
Particle production occurs as the two pancakes overlap leaving a "central" region of high energy 
density as they recede and carry away the net baryon density (leading baryon effect). If the system 
is formed after some time tj, the energy density in this region can be estimated following the 
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QCD Phase Diagram 



quark-gluon plasma 




Figure 2.19: The QCD phase diagram. The transition is believed to possess a critical point where 
the first-order phase transition between hadron gas to QGP changes to a smooth cross over. 



Bjorken procedure by considering the total energy per unit rapidity in nucleus-nucleus collisions, 
-. The energy density is then 



dy 



dE T A 1 
Az 



(2.57) 



For a small slice along the beam direction Ay = — , and V = irR A Az. R\ is the nuclear radius 
and behaves as Ra — 1.2 A 1 / 3 fm. The energy density can be expressed as 

1 dE T 



-KR?pt f dy 



(2.58) 



dE T 



For Au 197 collisions at RHIC, was measured by PHENIX to be 688 GeV at ^/s^ = 130 GeV [TT2j . 
If the formation time is on the order of 0.1 — 1 fm, this procedure gives an energy density estimate 
of 4.6 — 46 GeV/fm 3 . The Stefan-Boltzmann equation of state can be used to infer a temperature, 

'£(fr C )^ 1/4 



T 



a 



(2.59) 



30 21 

where a = — ^(16 + —Nf) is a constant of proportionality determined by the effective number 
of degrees of freedom. Since the lattice results indicate that this limit is not reached a more 
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appropriate estimate would be to use the lattice value of a = e/T 4 ~ 13, just above the transition 
region. This estimate gives a temperature of 128 < T < 229 MeV. The transverse energy density 
in Pb+Pb collisions at the LHC at t/snn = 2.76 TeV is larger than at RHIC by roughly a factor of 
5, indicating a temperature range of 191 < T < 342 MeV. As the temperature region of the phase 
transition is spanned by these estimates, it is likely that such a transition is probed in relativistic 
heavy ion collisions. 

If the mean free path in the interacting system is small compared to the system size a hydrody- 
namic description can be applied to the created nuclear matter [95] . The evolution of the system 
in the absence of viscosity or heat conduction is described by ideal hydrodynamics, 

where T^ v is the relativistic stress-energy tensor, is the local four-velocity of a fluid element 
and the system of equations is closed by a thermodynamic equation of state. The Bjorken picture 
considers the longitudinal expansion of this medium which leads to a slow decrease in the temper- 
ature, T ~ t -1 / 3 , where r is the local fluid proper time, and constant entropy inside the co- moving 
fluid volumes. 

Local anisotropics in the transverse density of produced particles can be produced through 
initial-state fluctuations or more commonly through the elliptical geometry of the collision zone 
when there is incomplete nuclear overlap. As the particle production is isotropic, the only way this 
can be converted to a momentum-space anisotropy, is if these density fluctuations form pressure 
gradients, which result in radial flow. Such an effect should be observable in the angular distribution 
of hadrons and may allow for a determination of the fluid's shear viscosity, rj. 

Calculational techniques in string theory have established a correspondence between quantum 
field theories and their dual gravity theories. This AdS / CFT correspondence facilitates a translation 
between quantities in a gravity calculation in the weakly coupled region to quantities in the strongly 
coupled regime of the appropriate dual [1 13(. 1114] , In particular the thermodynamics of black hole 
can be mapped to thermodynamics in QCD-like J\f = 4 supersymmetric Yang-Mills theory. Some 
time ago, arguments based on the uncertainty principle suggested that there is a quantum lower 
limit to the viscosity |115j . This has been demonstrated explicitly in the AdS/CFT picture where 
a limiting value of the viscosity to entropy density ratio was found to be rj/s < 1/4-7T [1 1117] . 
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This application of hydrodynamics illustrates how the transport properties of the system can 
provide further insight into the dynamics of the medium, in particular through transport coeffi- 
cients such as the shear viscosity, which are sensitive to the microscopic behavior. As a transport 
phenomenon, hydrodynamics coincides with the long-wavelength behavior of the system and the 
relevant transport and other phenomena such as radiation and diffusion can provide access to other 
aspects of the medium. 

2.5 Hard Processes in Nuclear Collisions 
2.5.1 The Glauber Model 

In collider experiments, the measured rates of quantities are proportional to the luminosity, or 
instantaneous flux of colliding particles. At high energies, the nuclear size is easily resolved, and 
a beam constituent not only sees a flux of nuclei but an instantaneous flux of nucleons/partons 
determined by the geometric overlap between the colliding nuclei. For two nuclei, A and B, with 
centers separated by impact parameter b, this flux is represented by the factor TAs(b). Many 
nuclear effects are expected to vary with b, and Tab is often used to normalize quantities to 
remove any variation due to trivial geometric effects. Experimentally, events are grouped into 
different centrality classes, which are determined to have a similar amount of geometric overlap 
based on a set of global event criteria. The average Tab is used to normalize quantities observed 
in those events to assess the centrality dependence of an observable. This procedure makes use 
of a framework known as the Glauber Model |118H12U] and is employed by nearly all heavy ion 
experiments [121J. 

The Glauber Model treats the incoming nuclei as smooth distributions of nucleons, each travel- 
ing on independent linear trajectories. It is formulated in the optical limit, in which the overall phase 
shift is evaluated by summing over all per-nucleon pair phase shifts. It relies on a parametrization 
of the nuclear density, p, typically a two-parameter Woods-Saxon distribution, 



where R and a are experimentally determined parameters, describing the radius and skin depth 
of the nucleus. The constant po is an overall normalization factor ensuring the distribution is 




(2.61) 
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normalized to the number of nucleons. For Pb 208 these parameters are R = 6.62 ± 0.06 fm 



a = 0.546 ±0.01 fm |122j and the distribution is shown in Fig 2.20 



Cl 




Figure 2.20: Woods-Saxon distribution for Pb 208 . 



In the Glauber formalism the transverse density, 



T A (b') = J dz A p A (b',z A ) 



(2.62) 



is the expected number of nucleons at position b'. Then Tab is defined as AB times the probability 



to simultaneously find nucleons in nuclei A and B at the same position b as shown in Fig. 2.21 

T AB (b) = J d 2 b'T A (h')T B (b' - b) . (2.63) 

Each nucleon in A can interact with a nucleon in B with probability pab = J^ N , so the total 



AB 



probability of n collisions is binomial and is given by 



(2.64) 



This allows for the definition of the expected number of collisions as N co \\ = Tab^nN' an< ^ * ne 
expected number of participants A^ part , the total number of nucleons that participated in any 
scatterings (sometimes referred to as wounded nucleons) as 



part 



+ 



d 2 b'T A (W) 
J d 2 b'T B (W) 



T B (b'-b) 
B 

T A (V - b) 
A 



B 



(2.65) 
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Figure 2.21: Views of the collision system transverse (left) and parallel (right) to the beam axis [121J. 



An alternative to performing the analytic integrals is to use Monte Carlo techniques. This 
has the advantage of including terms neglected in the optical approximation that incorporate local 
per-event density fluctuations [121^ 1123] . This method is performed by sampling the full Woods- 



Saxon distribution in Eq. 2.61 A times to populate positions for nucleus A. To prevent overlap, a 
position is regenerated if it is found to be within some minimum distance of a previously generated 
nucleon. Once the positions have been generated for both nuclei, a random impact parameter 
vector is generated defining an offset between the nuclear centers. The transverse position for all 
nucleons in nucleus A is compared to each of the analogous nucleons in B. If the distance between 



the pair is Ar < y °"nn / ^ the nucleons are considered to have participated. N co n is defined as the 
number of times this condition is satisfied, with iVp ar t defining the number of nucleons for which 
this condition was satisfied at least once. An example of an event generated with this technique 



is shown for an Au+Au collision in Fig 2.22 The distribution of A r part values for different impact 



parameters is shown in Fig. 2.23, as well as the relationship between N co \\ and N paiT t. 



Experimentally, the per-event impact parameter is not measurable so a procedure must be per- 
formed to relate a distribution of some measurable quantity to the Glauber parameters. Variables 
like the total charged particle multiplicity or transverse energy typically have distributions similar 
to the .A/part and N co u distributions. This feature, combined with the fact that these are global vari- 
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Figure 2.22: Distributions of nucleons generated with the MC Glauber procedure in the x — y 
(left) and z — x (right) planes for Au 197 . The two nuclei are shown in different colors, with the 
participants shown in a darker color |121j . 




b [fm] Npart/2 

Figure 2.23: Distribution of N pa _ Tt values as a function of b (left) and the correlation between N pa _rt 
and iVcoii (right) for Pb+Pb collisions. 
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ables and therefore less sensitive to detailed features and fluctuations make them excellent choices 
for centrality variables. The centrality determination procedure considers the minimum bias dis- 
tribution of such a variable, Q, and divides the range of observed values into sub-ranges where the 
integral of the distribution over that range is some percentage of the total; an example of this divi- 



sion is shown in Fig 2.24 These sub-ranges are called centrality intervals, and centrality-dependent 



observables are usually calculated by averaging over all events in the same interval. 
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Figure 2.24: Schematic diagram from Ref. |121j of a distribution of a centrality variable, in this case 
the number of charged particles in the interval \q\ < 1. The events contributing to the upper 5% 
of the integral of the distribution are the 0-5% centrality bin, with near complete nuclear overlap. 
The parameters from a Glauber simulation, (b), (-/V part ), are shown as different horizontal scales. 



The i centrality bin, Q < ( < Q+i, is typically defined in terms of percentages, a — 6%, such 



that: 



1 dN e 



c . dC: N cvt d( 

S>min 3 



a%, 



c™, i dNcvt 

d< ^M = b% ■ 



(2.66) 



The Glauber model parameters must then be related to the class of events in a particular centrality 
bin. This can be accomplished by constructing a new variable, £, from the Glauber variables N co \\ 
and iVp ar t such that the £ distribution is similar to the experimentally observed £ distribution. The 
£ distributions tend to be similar enough to the iV co n and iV part distributions that a simple linear 
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combination of the two, a two-component model, is a suitable choice for £, 

Z = to(x^f± + (l-x)N coi y (2.67) 

The parameters £o an d x can be determined from fitting the measured dN ev t/dC distribution with 
dN evt /dt; from an MC Glauber sample. The centrality bins for £ are defined in terms of integral 
fractions just as for £. In the i th centrality bin, the iV co ii and iV par t are averaged over all events for 
which £j < £ < These (N co u) and (iVp art ) values are then associated with the i th centrality 

bin in the data. 

2.5.2 Nuclear Modification 

The first evidence for modification in nuclear collisions at high energy was observed p+A collisions, 
where it was found that the production of particles at large transverse momentum ( 2 < px i$ 
6 GeV) was enhanced in these collisions relative to pp [1241 1125] . This effect, Cronin enhancement, 
was originally interpreted as additional transverse momentum, k 2 oc L, imparted by additional, 
independent, elastic interactions with multiple nucleons. Here L is the average path length, which 
scales with atomic number as L oc A 1 / 3 . 

The nuclear parton distribution functions (NPDFs), show both a suppression at lower values 
of x and an enhancement with increasing x with respect to the nucleon PDFs, termed shadowing 
and anti-shadowing respectively. This nuclear modification is quantified through the ratio, 

Rf = ^, (2.68) 

Ji 

where / are the PDFs for the nucleus A and the nucleon, and i is the parton species: valence quark, 
sea quark or gluon. When viewed in the rest frame of the target nucleus, the shadowing/ anti- 
shadowing is the result of multiple scattering processes that destructively/constructively interfere 
at the amplitude level [126j. In the collinear factorized approach coherent multiple scattering 
terms are suppressed by powers of 1/Q 2 , however for a large nucleus these contributions receive an 
enhancement of A 1 / 3 , leading to sensitivity to higher twist effects [127J. 

When viewed in the infinite momentum frame, the shadowing effects arise from recombination 
of low-x gluons in the nuclear wave function |128j . This phenomenon is known as saturation and 
sets in as unitarity requirements force the nominal evolution equations to be modified at small x. 
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The gluon distribution, xG(x, Q 2 ), represents the number of gluons per unit rapidity in a region of 
transverse size 1/Q 2 and grows as In 1/x for small x at fixed Q 2 without this modification. However, 
as Q 2 R 2 < xG(x,Q 2 ), where R is the radius of a nucleon, the gluons will begin to overlap and 
recombination of these gluons will limit the growth of xG(x, Q 2 ) with decreasing x [129\ Il30j. 

At larger values of x, away from the shadowing/anti-shadowing region (x > 0.5), the NPDFs 
exhibit additional suppression which is known as the EMC effect |131j . The origin of this effect 
may be related to non-nucleon degrees of freedom in the nucleus |132j and has led to increased 
interest in short-range correlations between nucleons. A sharp enhancement of the NPDF at the 
largest x values is thought to be described by the Fermi motion of the nucleons. Global analyses at 
NLO have been performed using data from DIS, Drell-Yan and d+Au collisions at RHIC to extract 
the x-dependence of R^, R^ and R^ at multiple Q 2 values [133 . These distributions, shown for 



Pb in Fig. 2.25, are a crucial input to any interpretation of any high-px phenomena observed in 
nucleus-nucleus collisions. 



Hi 



% 1.4 
$ 1.2 h 



7\ o.8 ^Htwwm4H^^^f«** l, *^ v 



M 

0.6 ^ 
& 0.4 

of 0.2 - 



_ 1.4 

> 12 

o 

2 0.8 
ii 

q. 0.6 

h 0.4 



0.2 

0.0 



Q=1.69GeV 



Rr, 



10 



4 io- 3 



10" 2 10"' 




10" 2 10"' 



Figure 2.25: Nuclear modification ratios R^ h , R^ h and i?g b for the valance, sea and gluon PDFs 
for two values of Q 2 |133j . For Pb the shadowing region is x < 0.05, the anti-shadowing region is 
0-05 < x < 0.5 and the EMC region is x > 0.5. 
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2.5.3 HIJING 

The phenomenology underlying the implementation of pp event generators such as PYTHIA, re- 
quires extensions to be appropriate for nucleus-nucleus collisions. An important step in achieving 
this was the development of the HIJING MC event generator [134] . To correctly model the multiplic- 
ity fluctuations the soft production mechanism combines elements of different models. It introduces 
multiple, sea qq strings as in the dual parton model |135H137j . but also allows for induced gluon 
bremsstrahlung by introducing string kinks as in the Lund FRITIOF procedure [1381 1139] . Addi- 
tional transverse momentum kicks are applied that are dependent on the particle's proportional 

to d4 + a 2 )- 1 d4+pg)- 1 . 

The production of "minijets" with pi > pq ~ 2 GeV is expected to play an important role in 
the total energy and particle production [140| 1141] . These processes take place at an intermediate 
momentum scale, lower than that associated with typical jet production, but still describable by 
perturbative QCD. Production of multiple jets in HIJING is implemented through a probabilistic 
model of independent jet (pair) production. The average number of minijets at impact parameter 
b is given by a^ e tT^(b), where <Xj e t is the inclusive cross section for jets in nucleon-nucleon collisions, 
integrated above some threshold po- The probability for multiple independent minijet production 
is given by |142j . 

g . {b) = MMf^Affl j> lt {2m) 

go (b) = [l - e - CTs ° ftTN W] e~ a ^^ , (2.70) 

where <7 so ft is the non-perturbative inclusive cross section for soft processes and ?n(&) is the partonic 
overlap function between two nucleons separated by impact parameter b. An eikonal picture is used 
to relate the inelastic, elastic and total cross sections and to define the multi-jet probabilities for 
minijets in nucleon-nucleon collisions, 

G = — Jo 00 d 2 b [1 - e - 2 xs(M] e -2xn(M ( 2 . 71 ) 
G = — XT d 2 6 [2XH(M] J e - 2xH(M (2 72) 

<7in r- 

Here the eikonal factors are related by xs + Xh = Xo(l + Cjet/o'soft)) and ensure geometric scaling. 
The independence of the multi-jet production is an appropriate ansatz if cjj e t < 2A~ 2 /' i (po-RN) 2 Cmei- 
The overall jet cross section is reduced from the nucleon-nucleon case to include the effects of 
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nuclear shadowing, including the impact parameter dependence which has not been measured by 



experiment. The modification factor defined in Eq. |2.68 can be generalized to include impact 
parameter dependence, 

R A (x,Q 2 ,e) = R° A (x,Q 2 ) + a(r)R s A (x,Q 2 ), (2.73) 



where a(r) oc (A 1 ^ — l)y 1 — r 2 /R A , models the impact parameter dependence |134j . Then the 
effective jet cross section is a function of the transverse positions of each of the colliding nucleons 
in the binary system, 

0jrt(rA, tb) = o£t + a A (rA)ojrt + "b^b)^ + a A (r A )a B (rB)cT^f . (2.74) 

The event generation proceeds by using an MC Glauber setup to determine the set of colliding 
nucleon pairs. For each of these binary collisions, the probability of scattering and number of jets 
is determined, along with whether the collision is elastic or inelastic. Hard scattering partons are 
treated separately, and their energies are subtracted from the nucleons with the remaining energy 
used in processing soft string excitations. The resulting scattered gluons are ordered in rapidity and 
color-connected to the valence quark/di-quark of the nucleon. The correlated semi-hard particle 
production mechanism is the key feature of HIJING and is why it has remained a useful tool long 
after its inception. 

2.5.4 Jet Quenching 

A significant form of nuclear modification occurs in heavy ion collisions where the byproducts of hard 
scatterings can interact with the QGP through the phenomenon of jet quenching. The character of 
this interaction provides key insight into the dynamics of the medium. The most important is the 
identification of the relevant scales for the jet-medium interaction and whether it can be described 
by perturbative QCD. If so, it provides a key testing ground for thermal, perturbative QCD and 
the HTL formulation as well as an example of the transport phenomena of radiation and diffusion 
in a fundamental physical system. 

The potential for jets as a tool to study the plasma was first recognized by Bjorken [143], who 
suggested that events with back-to-back jets would be sensitive to differential energy loss if the 
two partons had different in-medium path lengths. In extreme cases a highly energetic jet may 
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emerge while its partner deposits all of its energy in the medium, which would represent a striking 
experimental signature. 

The energy loss mechanism originally proposed by Bjorken was through elastic collisions with 
the medium constituents. This was originally given in Ref. [143] as, 

i^SO-f ■ 

with more intricate forms of the term inside the logarithm due to improvements in the collision 
integral given in Refs. |144H146] . In a QCD plasma, the Debye screening mass, mo, is given by, 

m D = (1 + l -N f )g 2 T 2 . (2.76) 

O 

Some of the first predictions of jet quenching signatures involve a modification of the dijet 
acoplanarity distribution in heavy ion collisions |147l 1148] . These calculations were formulated in 
terms of the probability distribution for dijet pairs to have momentum imbalance 

dP 1 1 da , m . 
= (2.77) 

dK v ct (pt) Pt dcj) ' 

The acoplanarity was expected to show a temperature-dependent modification due to elastic colli- 
sions with the medium constituents. It was shown that similar effects could be produced through 
collisional energy loss in a hadronic resonance gas |149j . and thus an observed modification would 
not prove the presence of a QGP phase. 

Interactions with the medium can also induce radiative energy loss through the emission of 

bremsstrahlung gluons. In QED, energy loss for high energy electrons is typically in the Bethe- 
dE E 

Heitler regime — — = — — , where L is a characteristic length. Classically, the total energy loss per 
ctx Li 

scattering is the integral of the bremsstrahlung spectrum 

AE= [ dojdk t 2 u dI (2.78) 
J duodkt 

with multiple scatterings adding incoherently to give a total energy loss 



dl 

duidk t 



AE tot = N J du}dk 2 U) ____^ ( 2 79 ) 



where N = L/X is the medium opacity and A is the mean free path between scattering centers. 

Most models use a formalism that treats the medium as a series of static scattering centers 
with the parton and radiated gluons with energies, E and ui respectively, traveling along eikonal 
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trajectories. This is combined with the kinematic limits q± -C uj <C E, where q± is the momentum 
transfer with the medium. This kinematic regime is referred to as the soft eikonal limit [150J. 
The radiation spectrum is typically derived by using the single gluon emission kernel as a Poisson 
probability for multi-gluon emission. 

At high energies, significant interference occurs between the parton and quanta emitted at small 
angles. This results in a finite formation time for the radiation and suppresses this contribution 
relative to incoherent radiation, known as the LPM effect |151| 1152] in QED. The QCD analog of 
this phenomenon has been proposed as an important feature of the quenching mechanism [153, 154J 
and was shown by BDMPS [155J and independently by Zakharov |156j to give an energy loss that 
grows quadratically with path length in the medium AE oc L 2 . 

As the interference suppresses the coherent radiation, the emission spectrum will be dominated 
by those quanta which have decohered. These are gluons which have acquired a phase, up, of order 
unity [T571fT58] . 

and thus appear with a characteristic energy, 

w c = \qL 2 . (2.81) 

Here the transport coefficient q has been introduced, which in this picture represents the mean 

squared transverse momentum imparted to the parton per unit length, q = (q±)/L. 

In this picture the energy loss is determined by soft multiple scattering. In the coherent limit 

the parton undergoes Brownian motion with a Gaussian elastic cross section oc — e~ q ±^ L , and the 

qL 

scattering centers behave as a single source of radiation. The bremsstrahlung spectrum is given by 

dl 2a s C R 



UJ 



duj 7T 



1 

Y2 V uj 



UJc. 

UJ < UJ C 

(2.82) 



UJ > UJc 



which results in energy loss AE ~ asuj c = -asqL 2 . The L 2 dependence is qualitatively different 

than both the Bethe-Heitler energy loss and the LPM effect in QED (uj— ~ 

duj 

Other models frame the energy loss as an expansion in opacity. These models are not re- 
stricted to small momentum transfers, including the power-law tail in the scattering cross section, 
but the coherence effects of BDMPS-Z formulation must be enforced order-by-order in opacity. 
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The GLV |159} 1160] model is an example of this approach, which models the scattering centers 

as screened Yukawa potentials with screening length, [i. The gluon spectrum is constructed by 

integrating first over the longitudinal direction to enforce the LPM interference at a given order. 

The momentum transfer q± is then averaged over giving the double-differential gluon distribution, 

uj ^ • This can be analyzed to give information about the transverse pattern of radiation before 

dujdkt 

integrating over kt and uj to give the full energy loss. 

The ASW formalism |158} 1161] is a path integral formulation that can be applied to both the 
multiple soft scattering (MS) or single hard (SH) dominated scenarios. In the SH approximation 
the gluon distribution agrees exactly with the GLV formula to first order in opacity, however the 
assumptions of the two models and ranges of integration differ, giving different results for the 
total energy loss. In the limit where these differences can be neglected (ignoring the kinematic 
constraints) the radiation spectrum is given by [158J, 

TTUJr. 



dl „, 2a s C R L) 4- 



UJ < UJ C 



w r-— t n > ( 2 - 83 ) 

d " * A I ln^ uj>uj c 

UJ 

where uj c = ls a characteristic frequency, different than the BDMPS-Z case. In the SH 

scenario, the radiation is dominated by uj > uJ c , also different from BDMPS-Z. The total energy 
loss, 

AE^kln^ (2.84) 

7T A UJ C 

E 

is enhanced by In — relative to the region uj < uj c . Despite the differences from BDMPS-Z, this 

UJ C 

limit also gives AE ~ L 2 . 

The ASW approach allows for the calculation of quenching weights, which give the probability 
for gluon splitting in the medium. With the assumption of an ordering principle in the virtualities of 
the fast parton, these weights can be combined with the vacuum DGLAP splitting functions. These 
modified splitting functions can be used to define a medium-evolved fragmentation function [162 
and can be used to calculate the effect of quenching on jet shape observables [163J . 

Some of the difficulties occurring in these analytic calculations can be alleviated by considering 
an MC approach. These typically supplement the generation of hard scattering processes used in 
generators by an additional step to include the quenching effects, and allow for the application 
of medium effects to the entire jet, not just the leading parton (democratic treatment) [164J. 
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Generators like Q-PYTHIA [165J, PQM fTHT] and JEWEL [168J, replace the vacuum parton 
showers implemented in PYTHIA with medium modified parton showers by altering the Sudakov 
form factors to include the quenching weights. The PYQUEN model imposes radiative energy loss 
of the BDMPS-Z type as well as collisional energy loss before applying vacuum fragmentation [169J. 

The approach that most directly incorporates the thermal of the system was developed by 
Arnold, Moore and Yaffe (AMY) [170H175] . This relaxes some of the assumptions used by other 
models, and includes thermal partons which are dynamical scattering centers. Furthermore, the 
requirement that the emitted gluons be softer than the initial parton energy (u <C E) is not 
required. The energy loss occurs through elastic scatterings with differential collision rates of the 
form 

dT _ 1 g 2 Tml, 



d 2 q± 2vr 2 q^(q2_ + 



in 



(2.85) 



D/ 



Multiple soft interactions are resummed using the HTL techniques discussed in Section 2.4.3 , leading 

to a differential emission rate for each parton species — . These are used to construct rate equations 

dk 

describing the evolution of the momentum distributions. This setup treats collisional and radiative 
energy loss in a consistent picture, but loses the distinction between radiated and thermal gluons. 
Additionally, finite length effects which result in finite formation times for the radiated gluons, are 
not included in the formalism. 

Although these formulations differ on a number of features, they do have some key elements 
in common. The most important of these is the operating assumption of factorization between 
the production cross section, energy loss process and fragmentation. While there is no rigorous 
proof of factorization in heavy ion collisions, experimental evidence supports the interpretation that 
quenching is not an initial-state effect and that final-state interactions are partonic in nature [150J. 
In all of these models, the quenching mechanism is not applied to the radiated quanta, thus the 
calculation is limited to the energy loss of the leading parton, except in the MC implementations. 
Furthermore, emitted quanta themselves are radiation sources. Coherence effects between these 
sources establish the angular ordering in the vacuum parton shower and allow for the cancellation 
of infrared and collinear divergences. Such effects may play an important role in medium-induced 
parton showers as well [176J. 

It should also be noted that all models use repeated application of an inclusive single gluon 
emission kernel, extending to multiple emissions either inferring a Poisson distribution or using 
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rate equations. However, this is not the same as applying the exclusive multi-gluon emission. 
Implementing this type of emission forces the parton kinematics to be updated or dynamical changes 
in the medium. Including these local effects as well as the global effects from the LPM-type 
interference poses a serious challenge [150] 

There are also issues with these approaches relating to the treatment of large angle radiation. 
The BDPMS-Z formalism receives its dominant contribution from radiation that is outside the 
valid kinematical range of the approximation |177j . In ASW and GLV, the transverse distribution 
of emitted gluons is calculable. For energy loss, this distribution is integrated over k% and the 
behavior at large angles is linked to how the kinematic constraints on the parton are enforced and 
the assumption of collinearity. It has been shown that the choice of maximum opening angle # max 
can change the results appreciably, leading to at the least a large systematic error \lbi)\ 1178] . The 



effects of different transverse cutoffs are shown in Fig. 2.26 The different physical pictures have 





dN 



calculated in the GLV formal- 



Figure 2.26: The double differential gluon distribution (left), . 

dxdk t 

ism at fixed x. Here, x = co/E, is the fraction of the quenched parton's energy carried away by the 
radiated gluon. The effect of different kinematic cutoffs in the kt integration is illustrated by the 
vertical lines. The effects of these cutoffs on the integration is shown on the right. Figures adapted 
from Ref. [150] . 



made it difficult to consistently fix parameters for a direct comparison of models. However a recent 
comparison was performed by applying the GLV, AMY and ASW (both MS and SH) formalisms 
to a "brick" of QGP matter of fixed length |15Uj . A comparison of medium-induced gluon emission 
spectra and nuclear suppression factors as a function of q for each of the different models are shown 
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Figure 2.27: Single gluon emission spectrum — as a function of gluon energy oj for the AMY, GLV, 

duo 

ASW-MS and ASW-SH formalisms. The calculations are for a 20 GeV parton passing through QGP 
bricks at T = 300 MeV with size L = 2 fm and L = 5 fm, shown on the left and right respectively. 
Figure adapted from Ref. [150]. 



Depending on the formalism, the quenching mechanism is sensitive to the system size as well 
as some combination of the intensive parameters q, \x and A, which provide information about the 
microscopic medium dynamics. The soft multiple scattering approximation is only sensitive to q, 
which is defined as the mean-squared momentum transfer per unit length, 

f ,o o da _ f qmax 
/o 

dr e \ 



2 2 dT e \ 
d q±q ± - 



(2.86) 



where 



d 2 qi 



J9 — / fill J9 5 

d 2 q T Jo d 2 q ± 
is the differential elastic scattering rate for a hard parton in a thermal medium. This 



1 



behaves as — r at high temperatures, but is screened at low temperatures by mi,, given by Eq. 

qi 

The form, 



2.76 



dT, 



el 



Cr 



d*q± - (2*) 2 qi(qi+f/if,)' (2 ' 87) 

interpolates smoothly between these limits |150j . and the number density N has been intro- 
duced [T79] . 

M= C ^ ) (l + \N f )T\ (2.88) 
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C(3) 

The leading coefficient is the ratio of Riemann zeta function values, — — ps 0.731. The value of q 
is given by, 

where g max is the largest transverse momentum transfer allowed in the elastic scatterings and is a 
function of T and in principle the parton energy as well. This introduces a logarithmic dependence 
of q on E, which is an effect of approximating the collision integral in the evaluation of the elastic 
scattering rate. This value is sometimes taken as <? max = g 2 ET although <? max = g 2 Vet* may be 
more appropriate [18(3] • The value, 

q « — - = A V ' oc T , (2.90 

A 47T 

and A calculated below, is a commonly used, energy independent expression and is equivalent to 
omitting the logarithmic variation of q. These estimates can differ by up to 40% for large parton 
energies. 

The opacity can be calculated as 



which in turn, determines A. However the form of the scattering rate here is slightly different 



than the one given in Eq. 2.87, and is chosen to be consistent with the model used in the opacity 
expansion framework [150] . This gives 



47rm? 



which is consistent with Eq. 2.90 



Motivated by new experimental capabilities at the LHC, extensions of energy loss calculations 
to full jets have recently been made |181H183] . These function by expressing the differential jet 
cross section as, 

where P q (x) is the probability a jet will lose a fraction x of its energy. The Jacobian |j7(x)| relates 
the initial and final parton E^ by E' T = \J q {x)\E^, and is given by 

\J q (x)\ = [l-x(l-f q (R,pf n )T 1 ■ (2-94) 
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Figure 2.28: Comparison between different models of the single hadron suppression factor, Rj, of 
a p^ 7 spectrum as a function of q. respectively. Figure adapted from Ref. |150| . 



The factor f q (R,p™ m ) is the fraction of the emitted radiation above p™ in and inside the jet cone 
defined by R, 



rE-T ^ dig 

'Pt 1il jiudr 



"Ft 00 j_ rE T j, , dI q 



(2.95) 



duidr 



While this form still only considers the energy loss of the leading parton, it does take into account 
the extent to which radiated energy is not lost, but only redistributed within the jet. Furthermore, 
it facilitates the usage of NLO and NPDF effects in the pQCD calculation of the unmodified jet 

d(Jn 



cross section, . The factorization of medium-induced radiation from the production cross 

dE T dy 

section, which is used in Eq. |2.93| has been proven in the context of an effective theory. This uses 
Soft Collinear Effective Theory (SCET) techniques to decouple the hard and soft components of 
the gluon fields and has been shown to reproduce the GLV expression when expanded to first order 
in opacity |184| . This factorization is also applicable to multi-jet systems, allowing calculations of 
differential energy loss in dijet systems. Calculations of the inclusive jet suppression factor and 
dijet asymmetry in this framework are shown in Figs. 2.29 and 2.30| > respectively. 



The AdS/CFT correspondence principle was discussed in Section 2.4.4, in the context of esti- 
mating the shear viscosity. Similar techniques have also been applied to study the energy loss of fast 
partons. While many of these results have focused on the specific case of heavy quarks [185-188J, 
some work has been done to provide an estimate of q [189H191] . These results agree qualitatively 
with the range of q estimates from experiment. 



CHAPTER 2. BACKGROUND 



53 



Q. 

if 



cc 



Q. 

if 



cc 



1.5 
1.25 
1 

0.75 
0.5 
0.25 


1.25 
1 

0.75 
0.5 
0.25 




1 1 1 1 1 1 1 1 

3med = 1 - 8 ' 

9med 

- 3m e d = 2 - 2 

min _ 

■ P T =o 


1 1 1 1 1 1 1 

ALICE ch. h., central Pb+Pb 
NLO + med. effects, R = 0.6 _ 
NLO + med. effects, R = 0.2 

Central Pb+Pb - 






V - 

1 I 1 I ■ I 1 I 1 I ■ I 1 I 1 


without CNM 

- with CNM 

No nuclear effect 


NLO + med. effects, R = 0.8 . 
NLO + med. effects, R = 0.6 _ 
NLO + med. effects, R = 0.4 
NLO + med. effects, R = 0.2 




9med 2 
1 , 1 


LHCs 1/2 = 2.76TeV - 
I 







50 



100 150 
E T [GeV] 



2(H) 



Figure 2.29: Calculation of jet Raa for various jet radii as a function of jet E^ (bottom). The top 
panel shows a comparison to single particle Raa- Predictions with and without NPDF effects are 
shown in the solid and dashed lines. Figure adapted from Ref. |183j 

2.6 Relativistic Heavy Ion Collisions: RHIC to LHC 

With the commencement of the RHIC program, many of the speculated features of heavy ion 
collisions became experimentally accessible. The collective behavior of the system can be studied 
by measuring the angular distribution of particles, ^— . As discussed previously, in the absence 
of hydrodynamical expansion, this distribution is expected to be isotropic. The anisotropy can be 
quantified by considering the Fourier decomposition of the distribution, 



dN 



f ll + J2 2v nCOs[n(<j) - * n )] 



(2.96) 



n=l 



where ^ n are the event plane angles, and v n describe the magnitude of the modulation. In all 
but the most central collisions, the direction of the impact parameter vector between the nuclear 
centers defines a reaction plane. The overlap region is ellipsoidal and symmetric about this reaction 
plane. The second Fourier coefficient, V2, is an observable with sensitivity to how the initial-state 
anisotropy is converted to the final-state particle distribution, a phenomenon known as elliptic flow. 
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Figure 2.30: Calculation of jet dijet asymmetry for various model parameters using kinematic cuts 
chosen to match measurements from ATLAS (left) and CMS (right). Figure adapted from Ref. [183 



Measurements from the RHIC program show that the elliptic flow, in particular the magnitude of 
the azimuthal modulation of particles as a function of their pi , is well-described by ideal, relativistic 
hydrodynamics. The near-vanishing viscosity indicates that the system is strongly coupled. This 
came as a surprise when compared to the prediction of asymptotically free quarks and gluons and 
which has led to monikers such as the strongly-coupled QGP (sQGP) and "the perfect liquid" 
applied to the medium created at RHIC. Since early RHIC results, significant advances have been 
made in the formulation and numerical implementation of viscous hydrodynamics. RHIC and LHC 



V2 values are compared with a recent viscous hydrodynamical calculation in Fig. 2.31 



More recently, it was recognized that the event-by-event fluctuations in collision geometry could 
drive higher flow harmonics. These provide additional information about the initial-state geometry 
and may be able to provide additional constraints on the hydrodynamical formulation or the value 



of the shear viscosity. The first six harmonics as measured by ATLAS |195j are shown in Fig. 2.32 



2.6.1 Hard Processes 

Indirect experimental evidence for jet quenching was first established by two important measure- 
ments. The single hadron spectrum was found to be heavily suppressed at high-pT in central 
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Figure 2.31: Values of i>2 as a function of pt (left) and centrality (right) measured by the STAR 
(Au+Au, ^/snn = 200 GeV) [TM| and ALICE (Pb+Pb, v / s N ¥ = 2.76 TeV) [US] collaborations are 
compared with values calculated from the VISHNU viscous hydro model with hadron cascade [194J. 
The elliptic flow is measured using a four-particle cumulant technique. 



collisions [196H199] . Furthermore, the away-side correlation in the dihadron angular distribution 
was found to be heavily modified in central Au+Au relative to pp [2001 ? . 1201] . These two measure- 
ments provide slightly different handles on the phenomenon as the former is an inclusive measure- 
ment of the total effect of the suppression on the pt spectrum. The latter provides a differential 
measurement of the quenching as it is sensitive to the quenching of one jet relative to another. 

The interpretation of these observations as evidence for quenching was further supported by 
results from the d+Au run which established a critical baseline. The lack of suppression at high- 
Pt |202l 203J combined with the unmodified dihadron correlation indicated that modification of 
the NPDFs was not responsible for the observations in Au+Au. Subsequent measurements of the 
direct photon Raa were consistent with no suppression and provided further systematic control on 
the quenching effects [204J. 

The degree of suppression is typically quantified by the nuclear modification factor, Raa, 



E- 



d 3 71- 



cent 
jet 



Raa 



TaaE- 



dp 3 

d 3 alt p ' 



(2.97) 



jet 



dp 3 



CHAPTER 2. BACKGROUND 



5G 



10- 1 



10 



i-2 



> 

o 

Q. 

CM 



= 1-2 
> 

0. 1 
LD 

0.8 



I ' ' ' ' I 



1<p*,p^<2 GeV 



11 11 1 " 11 1 

ATLAS 



Pb-Pb ys NN =2.76 TeV L int = 8 |xb"' 



gjifliiiiji^ 

^AAAAA A AAAa 



EP 2PC 
. -e- V 



« * ' * 5 « 

t * £ 



v 4 I 

V 6 





+ + 



2<p a .,p^<3 GeV 



1 1 " " I 



e s 



8e ®«®®®® 8(e 



H a h b a 

5 A AAA A 



h ■ ■ B 

A A A A 



® 

9 
■ ■ 

A A 



5 * * * * ♦ ♦ 



* * * $ 
+ 



*4. 



O O O O O O g 

--ti°- S 




80 60 40 20 80 60 

Centrality [%] 



40 



20 



Figure 2.32: Values of vi—vq as a function of centrality measured by ATLAS |195j . Two techniques, 
two-particle correlations (open markers) and the event plane method (solid markers), were used 
and were found to be in good agreement. The ratio of the results from the two methods is shown 
in the bottom panel. 



A summary of the Raa measurements from PHENIX is shown in Fig. 2.33| These results established 
that the production of high-px particles is strongly influenced by the produced medium in heavy 
ion collisions. However, the single particle measurements (with the exception of the photon), are 
limited in their utility. The medium effects involve a high momentum parton, not a final-state 
hadron. Thus single particle observables can only be connected to the process of interest through 
a fragmentation function. This necessity forces the interpretation of the results in the context of 
strict factorization between the medium effects and fragmentation. While this factorization may 
ultimately prove to be an appropriate assumption, it is desirable and more objective to work in 
a more general paradigm for the jet-medium interaction, namely one that does not enforce strict 
separation between medium effects and the jet fragmentation. Furthermore, the single hadron 
observables are only linked to the parton level quantities in an average sense; there is no guarantee 
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Figure 2.33: PHENIX -Raa measurements for 7r°, rj, 0, w, J/-0 and direct 7 [205J. 



that the highest energy hadron in an event came from the leading jet, or the jet that suffered the 
least energy loss. This limitation is removed if per-jet fragmentation distribution is eliminated from 
the observable. 

It is less restrictive, although more experimentally challenging, to construct observables from 
fully reconstructed jets. These quantities provide direct sensitivity to quenching effects. Further- 
more, the possibility of using jets as input objects into physics analyses opens many new possi- 
bilities, which will be discussed in detail in Chapter [7j Measurements of full jets at RHIC have 
been attempted, but are limited by a number of factors and as of this time no such measurement 
has been published. The main complication is due to the production rates of jets that are easily 
detectable above the medium background. Measurements are further constrained by the limited 
acceptance of the PHENIX and STAR detectors. Although highly collimated, particles from jets 
can still be distributed over a substantial angular range. Typical sizes for the cone radii and R 
parameters in the sequential clustering algorithms used in pp experiments are on the order of 0.5. 
Thus larger acceptance detectors are preferred for jet measurements. 

Many of these issues are not present at the LHC, in particular the ATLAS experiment, discussed 
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at length in Chapter [3| At these high energies, the first Pb+Pb run was at <tnn = 2.76 TeV, the 
rate for producing jets well above the background from the underlying event is much higher; the 
differential single jet inclusive cross section da/dE^dij ~ 1 nb for 100 GeV jets in pp collisions at 
these energies. Furthermore, the high quality calorimetry covering 10 units in r\ enables precise 
measurements of jets and their properties. 
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Chapter 3 



Experimental Setup 



3.1 The Large Hadron Collider 

The Large Hadron Collider (LHC) is a particle accelerator at CERN outside of Geneva, Switzerland. 
Although it was primarily designed to collide protons, the machine is also capable of colliding heavy 
ions. This program began in November 2010 with the first lead ion collisions. 

The LHC machine [206J was constructed using the tunnel originally used for the LEP experi- 
ment. It consists of two parallel beam lines circulating particles in opposite directions, intersecting 
at designated interaction points (IPs). The ring is 26.7 km in circumference and contains eight 
arcs and straight sections. The particle orbits are primarily controlled by the 1232 dipole magnets, 
while strong transverse focusing of the beam constituents is maintained by alternating field gradi- 
ents supplied by 392 quadrupole magnets. Almost all of these devices use superconducting NbTi 
cabling operating in a cryogenic system maintained by superfluid He II at 1.9 K. 

Pb 208 ions are extracted from a source and processed by a sequence of injection chain elements 



before being injected into the LHC. This sequence is shown in Fig. 3.1 consists of: Linac3, the 
Low Energy Ion Ring (LEIR), the Proton Synchrotron (PS) and the Super Proton Synchrotron 
(SPS). In this process the ions are stripped of electrons, squeezed into longitudinal bunches via the 
application of a radio frequency electric field (RF) and accelerated to a beam energy of 177 GeV 
per nucleon. 
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Figure 3.1: LHC ion injection chain. 



3.1.1 Performance and Luminosity 

The most important factor in determining a collider's performance is the luminosity. The instan- 
taneous luminosity of a particle beam, ££ , is the flux of scattering particles per unit area per unit 
time. For a process with cross section a the interaction rate is given by 

(31) 

Therefore to enhance the rate of rare processes, it is a design goal of an accelerator to maximize 
this quantity. The colliding beams consist of bunches of ions with densities n\ and ni. Along the 
beam direction, s, the bunches can be organized into various patterns according to an injection 
scheme in which each colliding bunch pair collides regularly with frequency /. In the transverse 
direction, the beams have profiles characterized by a x and a y . For iV& bunches the instantaneous 
luminosity can be computed via 

J? = N b f-^. (3.2) 
47ra x a y 

The beam's transverse focusing is controlled by applying alternating gradient fields forcing the 
beam to converge. The transverse motion of single particles within the beam is then a sinusoid 
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with an envelope modulated by the beta- function, f3(s) 

x(s) = A^/f3(s) cos(ip(s) + 5) (3.3) 

where A and 5 are integrals of motion, describing the area and angle of an ellipse in the phase space 
(x, x' = dx/ds) inhabited by beam particles |30j. To maximize luminosity in the neighborhood of 
the interaction zones, the beams are squeezed in the transverse direction by focusing magnets [207J. 
The squeezing is parabolic in the longitudinal displacement and is controlled by the parameter f5* 

P*2 



a(s) = a(0) ( 1 + i ) ■ (3.4) 



In addition to /?*, the other parameter affecting luminosity is the emittance, e. In the absence 
of beam losses, the phase space ellipse has fixed area ttA 2 . For a beam with a Gaussian transverse 
profile, the emittance is the area containing one standard deviation a x : 

6x = 7T f x - (3 ' 5) 
In terms of these parameters, the luminosity can be expressed as 

* = Nbf A • ( 3 - 6 ) 

The longitudinal beam structure is governed by the bunch injection scheme |208| . The SPS 
is capable of injecting proton bunches at 25 ns spacing or 40 MHz into the LHC, corresponding 
to a total 2,808 bunch crossing slots (BCIDs). As of yet, this design limitation has not been 
reached. The end of the 2011 proton run injected 1380 bunches, 1331 of which were brought into 
collision for a luminosity of 3.65 x 10 33 cm s . The 2010 Pb ion run saw a maximum luminosity 
of 2.88 x 10 25 cnr 2 s _1 corresponding 121 bunches per beam with 113 colliding bunches and a 
500 ns bunch spacing shown in Fig. |3.2| For the entire 2010 Pb ion run, the LHC delivered a 
total integrated luminosity of 9.69 fib^ 1 . This quantity, as well as the total integrated luminosity 



recorded by the ATLAS detector as a function of time is shown in Fig. 3.3 Although not analyzed 
here, the 2011 Pb ion run recorded a total integrated luminosity of approximately 140 /ib _1 . 

In addition to the discrete bunch structure, RF cavities are used to focus the bunches in the 
longitudinal direction. The 400 MHz oscillating field causes particles to clump together around the 
ring in each of the 35640 RF buckets. Superimposing the RF structure with the injection bunch 
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500ns_1 21 b_1 1 3_1 1 4_0_4bpi31 inj IONS 




Figure 3.2: Injection scheme for the 500 ns bunch spacing fills in 2010 ion run. The blue and red 
points correspond to the filled BCIDs in beams 1 (clockwise) and 2 (anti-clockwise) respectively. Of 
the 121 filled bunches 113 are configured to collide. This scheme was used to achieve the maximum 
instantaneous luminosity for Pb ions in 2010: 2.88 x 10 25 cm s . 

spacing causes 10 RF buckets to be associated with each BCID. Ideally all particles in the injected 
bunch are forced into a single RF bucket, however a few stray particles become trapped in adjacent 



buckets as indicated in Fig. 3.4 Collisions among these particles, called satellite bunches, will 



appear out of time with the LHC clock in discrete 2.5 ns intervals. 

3.2 ATLAS Overview 

The ATLAS experiment is a multi-purpose particle detector situated at interaction point 1 (IP1) 
of the ATLAS ring [209] . It is forward-backward symmetric covering the full 2ir in azimuth. 
Charged particle tracking is provided by the inner detector, covering rj < 2.5, immersed in a 



2 T solenoidal magnetic field which is shown in Fig. 3.6 Energy measurements are provided 
by a combination of electromagnetic and hadronic calorimeters covering \rj\ < 4.9 enclosing the 
inner detector. A dedicated muon spectrometer is positioned beyond the calorimeter utilizing a 
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Figure 3.3: Total integrated luminosity delivered by the LHC (dark blue) and recorded by the 
ATLAS detector (light blue) as a function of day during the 2010 Pb ion run. The total delivered 
luminosity was 9.69 fib^ 1 . 



toroidal field maintained by a barrel and two end-cap toroidal magnets. This system allows muon 
measurements over the range \rj\ < 2.7 utilizing a variety of subsystems. Forward detectors such as 
the Zero Degree Calorimeter (|ry| > 8.3) and Minimum Bias Trigger Scintillators (2.09 < \i]\ < 3.84) 
provide minimum bias event triggering and event selection capabilities. A graphical representation 
of ATLAS highlighting the various subsystems is shown in Fig. |3.5| 



This analysis makes extensive use of the calorimetry, which is described in detail in Section 3.4 



The remaining systems are described in less detail in Sections 3.5, 3.6 and 3.7, except the muon 
spectrometer which is not used at all in this analysis and is mentioned for completeness. 



3.3 Trigger 

A trigger is a combination of hardware and software elements designed to select which collision 
candidates are recorded by the data acquisition (DAQ) system. When few colliding bunches are 
circulating in the machine, a simple requirement that the given pair of bunches is intended to 
collide in the filling scheme is the most basic criterion for triggering. However this is insufficient 
as many of the recorded events may contain no real collisions or detector signals from beam and 
cosmic backgrounds. Additionally, at higher luminosity the DAQ typically cannot record events at 
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Figure 3.4: A schematic view of the RF bucket structure corresponding to a single bunch fill. The 
red point indicates the filled bucket in time with the LHC clock. The blue points indicate satellite 
bunches which can cause out of time collisions. 



this rate and a trigger selection must be employed to ensure that only events with desired physics 
signals are recorded. In general, implementing a sophisticated trigger strategy to solve this problem 
is a significant experimental effort, requiring multiple triggers sensitive to different physics signals 
as well as using prescales to reduce the rate of less interesting triggers relative to the rarer ones. 
However, the luminosity of the 2010 ion run was low enough that a set of minimum bias trigger 
items could be selected without need for prescale. 

The ATLAS trigger system is composed of trigger items on three different levels. The Level 
1 (LI) trigger is entirely hardware-based. In addition to the normal data readout path, detectors 
integrated into LI have parallel paths for data readout, often involving coarser and faster signal 
sums. The specialized trigger signals are sent to the Central Trigger Processor (CTP), which 
combines the information into a set of LI bits. The other two levels Level 2 (L2) and Event Filter 
(EF) are both software based and are collectively referred to as the High Level Trigger (HLT). LI 
regions of interest (ROIs) are used to seed algorithms run as part of L2. Finally, L2 items feed full 
scan EF algorithms, which use higher level detector signals from the offline readout path and are 
not constrained to the LI ROIs. A trigger chain is defined as an EF item that is seeded by various 
LI and L2 items. Multiple chains can be seeded by the same LI and L2 items and prescales can 
be applied at any of the three levels. Trigger chains are grouped together by type to form data 
streams; an event is recorded and reconstructed if it has been selected by one or more of its chains 
after all prescales have been applied. 

During the 2010 ion run, a single minimum bias stream, physics_bulk, was used for physics 
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Figure 3.5: A diagram of the ATLAS detector showing the major detector systems. 

analyses. While some chains using HLT items were used for efficiency and background studies, this 
stream was primarily composed of raw LI items passed through the HLT without prescale. 



3.4 Calorimetry 

The ATLAS calorimeter system is composed of electromagnetic and hadronic calorimetry using 
sampling calorimeters based on two distinct technologies: liquid Argon (LAr) and scintillating 
tiles. The r\ coverage and segmentation of each subsystem is summarized in table [37Tj Energies of 
particles well above the ionization regime are measured by a sampling technique. Electromagnetic 
and hadronic showers are initiated when the particle strikes an absorber. The details of these 
showering mechanism are discussed in Sections |3.4.1 and 3.4.2 respectively. As the showers develop 



the energy of the incident particle is increasingly spread among an ensemble of particles of lower 
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energy. An active material is placed behind the absorber which collects some fraction of energy of 
the lower energy particles, either through ionization (LAr) or scintillation (tile). Alternating layers 
of absorber and active material are placed in succession and the shower-sampling is repeated. A 
pictorial representation of the different calorimeters and the sampling of the different shower types 



is shown in Fig. 3.7 



The EM barrel and end-cap sub-detectors constitute the first sampling layers over the central 
portion of the detector (\r]\ < 3.2) and possess fine segmentation for high-precision measurements. 
Both of these sub-detectors use sheets of steel- reinforced lead as absorbers, between 1.53 and 1.7 mm 



thick, folded into an accordion shape as shown in Fig. 3.8 The sheets are stacked and interleaved 
with readout electrodes positioned in the middle of the gaps on honeycomb spacers. The barrel 
consists of two separate half-barrels (z < and z > 0) each containing 1024 sheets with the 
accordion pattern extending radially, stacked to form cylinders covering full 2tt in azimuth. The 
cylinders each have a length 3.2 m with inner and outer diameters of 2.8 m and 4 m respectively. The 
end-caps consist of two co-axial wheels, of inner and outer radii 330 mm and 2098 mm respectively, 
composed stacked of absorber sheets and interspaced electrodes with the accordion pattern running 
parallel to the beam direction. 

The absorber structures are sealed inside cryostats filled with liquid Argon at 88.5 K. The ab- 
sorbers initiate electromagnetic showers, with the lower energy particles in these showers producing 
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Type 


Sub-detector 


Absorber 


Layer 


-range 


At] x A0 


Channels 


LAr 


Barrel Presampler 


None 


1 


|r?| < 1.52 


0.025 x 0.1 


7808 


EM Barrel 


Steel 


1 


\r)\ < 1.475 


0.003 x 0.1 


57216 


2 


0.025 x 0.025 


28672 


3 


|r,| < 1.35 


0.050 x 0.025 


13824 


End-cap Presampler 


None 


1 


1.5 < M < 1.8 


0.025 x 0.1 


1536 


EM End-cap 


Steel 


1 


1.375 < \r}\ < 3.2 




28544 


2 


0.025 x 0.025, 0.1 x 0.1 


23424 


3 


1.5 < M < 2.5 


0.050 x 0.025 


10240 


Hadronic End-cap 


Copper 


1 


1.5 < M < 3.2 


0.1 x 0.1, 0.2 x 0.2 


3008 


2 


2624 


FCal 


Copper 


1 


3.1 < |r?| < 4.9 


~ 0.2 x 0.2 


1008 


Tungsten 


2 


500 




254 


Tile 


Tile Barrel 


Steel 


1 


M < io 


0.1 x 0.1 


5760 


2 


3 


0.2 x 0.1 


Tile Extended 


Steel 


1 


0.8 < |r/l < 1-7 


0.1 x 0.1 


4092 


2 


3 


0.2 x 0.1 



Table 3.1: Description of coverage and segmentation of each calorimeter sampling layer. The 
Ar] x A(f>, correspond to the segmentation applying to most of that layer, although not necessarily 
constant over the full layer. The segmentation of the EM end-cap, with 77-dependent segmentation in 
Ar] but constant A(j) = 0.1, and the FCal which has non-projective geometry onto r\ — (f> coordinates, 
are not given. Additional tiles (not shown) are interspersed irregularly in the gaps between the 
cryostats and support structures. 
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ionization electrons in the LAr. In addition to the sampling layers, separate "presampler" layers 
of LAr are placed in front of the first barrel and end-cap layers. These modules allow for the col- 
lection of energy from electromagnetic showers that start early due to the material in front of the 
calorimeter. The amount of material, as indicated by the number of radiation lengths as a function 



of 77 is shown in Fig. 3.9 The sharp peaking in the upper left figure is caused by the edge of the 
cryostat, which projects to a narrow range in rj. 

The electrodes are held at high voltage and collect the drifting ionization electrons as a signal. 
This signal is processed by cold electronics mounted inside the cryostat. These signals are passed out 
of the cryostat through dedicated signal feed-throughs to a front-end crate, which provides signal 
amplification, pulse shaping and analog-to-digital conversion (ADC), and output with optical links 
to the main DAQ system. The front-end electronics possess an independent readout path for the 
LI trigger where analog sums of adjacent cells are taken before ADC. The trigger readout possesses 
a dedicated readout path. This has the benefit that a correction for "dead cells" in the offline 
readout due to a failure of the optical link is possible using the coarser trigger information (dead 
OTX correction). 

The high luminosity requires short drift times (~ 450 ns). In the barrel the geometry allows 
for gaps of fixed size between the absorber and electrodes, 2.1 mm, which corresponds to an op- 
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Figure 3.8: Accordion structure of the EM barrel (left) shown as a cross sectional slice transverse to 
the beam direction with particles incident from the left. The region between two sheets is zoomed 
in to show the position of the electrodes and the liquid Argon gap. A photograph of the accordion 
structure of the electromagnetic end-cap is shown on the right. 



erating voltage of 2000 V. The end-cap geometry does not allow this, thus the applied voltage is 
77-dependent, varying from 1000-2500 V. The signal current is proportional to the electron drift 
velocity, which is both a function of the applied voltage and the temperature of the cryogenic 
liquid, thus the temperature and purity are monitored closely. Measurements of the drift velocity 



performed using cosmic muons prior to the start of LHC operations are shown in Fig 3.10 good uni- 
formity was found in the measured drift velocities at the operating voltages and temperature |210j. 



The hadronic end-cap calorimeter (HEC) also uses liquid Argon as an active medium, however 
it uses copper plates to increase the number of interaction lengths and thus functions as a hadronic 
calorimeter. The end-caps are ±z symmetric, and each side is composed of a front and rear wheel 
module. The copper plates are oriented perpendicular to the beam axis and there are 24 plates each 
25 mm thick in the front wheel modules, and 16 plates each 50 mm thick in the rear modules. The 
wheels each have a front plate of half the nominal plate thickness, 12.5 and 25 mm respectively. 
The LAr gaps are fixed at 8.5 mm, and are divided by three electrodes into four drift zones of 
1.8 mm. The readout cells are defined by pads which are etched on the surfaces of the electrodes. 
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Figure 3.9: Material budget as a function of rj in front of various layers of the calorimeter in terms 
of the radiation length Xq. The total amount of material before the calorimeter and presampler 
is shown on the upper right. A complete breakdown of the number of radiation lengths before the 
presampler, before accordion sampling and after accordion sampling is shown in the upper right. 
The number of radiation lengths broken down by calorimeter sampling layer is shown for the barrel 
and end-cap in the bottom left and right plots respectively. 



A schematic drawing of the HEC including dimensions is shown in Fig. 3.11 



The final element of LAr-based technology are the FCal modules. These modules cover 3.1 < 
\rj\ < 4.9 and are divided into three longitudinal layers. All three layers use absorber plates 
perpendicular to the beam axis. These plates have tubes bored into them in the beam direction. 
Each tube is filled with a rod containing a coaxial cathode and anode electrodes. The space between 
electrodes is filled with LAr forming a small, 0.249 mm gap. The first layer is an electromagnetic 
calorimeter and uses copper for the absorber. The second and third layers have similar geometry 
but are constructed from 2.35 cm copper front and back plates which support tungsten rods. The 
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Figure 3.10: Distribution of electron drift velocities in the second electromagnetic barrel and end- 
cap layer (left) as functions of ij. Black points denote the mean drift velocity. The //-averaged drift 
velocities for each of the electromagnetic barrel layers is shown on the left. 
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Figure 3.11: Schematic diagram of the HEC showing views in the R — (left) and R — z (right). 
Dimensions are in mm. 



region in between rods is filled by filleted tungsten slugs to maximize the number of interaction 
lengths. The layout of the absorber matrix for the first FCal layer is shown in the left panel of 



Fig. 3.12 The FCal, HEC and EM end-cap modules are all situated inside the same cryostat with 



the positioning of the FCal modules shown on the right of Fig. 3.12 



The remainder of the calorimetry in ATLAS is provided by the tile system, which uses alternat- 
ing steel plates and polystyrene scintillating tiles. This system is composed of central (|r/| < 1) and 
extended (0.8 < < 1.7) barrel sections. Each consists of 64 self-supporting, wedge-shaped mod- 



ules, partitioned in azimuth. A diagram of one such tile drawer is shown in Fig. 3.13 Scintillation 
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Figure 3.12: The structure of the first FCal module is shown on the left. The red circles indicate 
the LAr gap. The Moliere radius is indicated by a pink circle, showing the scale of a single shower 
relative to the sampling matrix. The positioning of the FCal modules within the end-cap cryostat 
is shown on the right. 



light is collected in wavelength-shifting fibers and fed to photomultiplier tubes mounted at the edge 
of the modules. The fibers are grouped to provide three radial sampling depths of thicknesses 1.5, 
4.1 and 1.8 interaction lengths (at r] = 0). In addition to the nominal central and extended barrel 
regions, a series of tile scintillators are placed in the geometrically-irregular gap regions between the 
cryostat. The segmentation of the tile modules and positioning of the gap scintillators is shown in 



Fig. 3.14 A summary of the amount of hadronic sampling as a function of calorimeter layer an 77, 



as indicated by the number of nuclear interaction lengths (see Section 3.4.2) is shown in Fig. 3.15 



3.4.1 Electromagnetic Showers 

In the low to moderate energy regime, energy loss is described by the Bethe formula which describes 
the interaction of charged particles and material through the excitation and ionization of atoms. 
Typically the lost energy is extractable by collection of ionization electrons or scintillation light. 
Above some energy, called the critical energy E c , the energy loss is dominated by bremsstrahlung 
radiation, or pair production in the case of photons. The length scale associated with the energy 
loss rate is known as the radiation length, Xq, and describes the mean distance over which a particle 
will lose 1/e of its initial energy through radiation. This quantity depends on the detailed structure 
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Figure 3.13: Diagram of a tile drawer. 



of the material including the atomic charge screening effects on the bremsstrahlung cross section 
as well as the density of scattering centers. It has approximate dependence on the nuclear charge 
and atomic number of the material as Xq ~ AZ~ 2 . The photon pair-production cross section is 
controlled by physics at a similar scale, and at high energies Xq is approximately 7/9 of the mean 
free path for a photon to travel before pair production. 

When a high energy electron or photon enters material, the induced radiation or splitting results 
in the initial energy being distributed among two lower energy particles. These two particles lose 
energy through the same mechanism, and this process continues developing an ensemble of particles 
among which the initial energy is shared. This process is known as an electromagnetic shower or 
cascade. Once particles in this ensemble become lower than the critical energy, they no longer 
contribute to the shower's development, which eventually dies off. The longitudinal development 
of the shower can be described by a simple model. Since each step in the shower development 
is binary, after a number of interaction lengths t = x/Xq the shower will contain approximately 
2* particles, and the average energy is E = Eo/2 l , where Eq is the incident electron energy. The 
system should freeze out when the mean energy is of order E c , which corresponds to a shower depth 
of 

X ~ Xo In 2 • (3 ' 7) 
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Figure 3.14: The segmentation of the central (left) and extended (right) tile barrel calorimeters. 
The segmentation is chosen to be approximately projective, and lines of constant 77 are shown. 



More sophisticated models describe the average evolution by a gamma distribution, which allows for 
more parameters in describing the medium dependence of the shower. In the transverse direction, 
the shower width is well described by the Moliere radius, Rm = XqE s /E c , where E s k, 21 MeV. 

3.4.2 Hadronic Showers 

Electrons and photons, as well as particles that decay into them such as the 7r°, are measured 
exclusively through electromagnetic showers. Occasionally 7r =t will undergo quasi-elastic charge 
exchange, resulting in an electromagnetically detectable 7r°. The bremsstrahlung rate for heav- 
ier particles is suppressed by the large particle masses. Energy measurements of other particles 
typically occurs through hadronic showers, which are considerably more complicated than their 
electromagnetic counterparts. The exceptions are muons and neutrinos, neither of which interact 
hadronically. The former are typically detected as minimum ionizing particles and often require a 
dedicated muon tracking system to measure. The neutrino energy can be inferred from so-called 
missing E-j variables in hermetic detectors. 

Inelastic collisions of hadrons incident on atomic nuclei result in a wide range of byproducts 
resulting in large fluctuations in the measurable energy in a hadronic shower. As the shower develops 
some fraction of energy becomes electromagnetically visible in the form of produced 7T°'S, which 
start their own electromagnetic showers. While some of the collisional debris may be detectable 
through ionization, many of the struck nuclei become fragmented resulting in spallation neutrons, 
which are difficult to detect. Although some neutrons thermalize producing photons which can 
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Figure 3.15: Number of nuclear interaction lengths as a function of r\ with the contribution from 
each sub-detector and layer shown as a different color. The contributions without labels (blue and 
light brown) are due to material. An increase is seen near the cryostat boundaries at \r\\ ~ 1 and 



\rj\ ~ 2.8. 



be measured through late neutron capture, much of the shower energy in spallation neutrons and 
nuclear recoil cannot be collected. A material's ability to initiate hadronic showers and appear 
opaque to hadronically interacting matter is given by the nuclear interaction length Ai, which 
describes the mean free path for undergoing inelastic hadronic interactions in a medium. 

Hadronic showers typically develop a substantial electromagnetic component. Thus hadronic 
calorimeters have separate responses to electromagnetic and hadronic processes. If these responses 
are significantly different, an effect known as non-compensation, fluctuations in the electromag- 
netic/hadronic composition in a shower will significantly increase the energy resolution and cause 
a response that is not linear in the incident particle energy. Both of these features are undesirable 
and represent key constraints when designing a hadronic calorimeter. Typically these features are 
minimized by decreasing the electromagnetic sensitivity (lower Z) while increasing the hadronic 
(higher A). 
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3.5 Minimum Bias Trigger Scintillators 

The ATLAS Minimum Bias Trigger Scintillators (MBTS) detector is used to both select collision 
candidates online and reject background events in offline event selection. This detector is composed 
of 32 modules, each a 2 cm thick polystyrene scintillator embedded with wavelength shifting fibers 
for readout. The modules are located ±3.6 m from the detector center with a total coverage of 
2.09 < \rj\ < 3.84 and 2ir in azimuth. Each side (A and C) contains 16 modules in two groups in 
77: 2.09 < \rj\ < 2.82, 2.82 < \r/\ < 3.84. Each group contains 8 wedges which together span the full 
azimuth. 

A trigger readout is implemented by applying a leading-edge discriminator to the signal pulse 
sending a hit to the CTP for each of the 32 modules that is over a threshold. Additional LI items 
are built from these bits including the coincidence triggers, MBTS_N_N, which are fired if N or more 
modules fired on each of the A and C sides. 

Aside from the LI trigger, the MBTS is primarily used offline for timing. From the pulse 
sampling, a time measurement relative to the LHC clock can be determined for each side. The 
MBTS AiMBTS = t-A — tc can be used to reject out-of-time signals corresponding to non-collision 
background or collisions between satellite bunches. 

3.6 Zero Degree Calorimeters 

Positioned ±140 m from the interaction points are the ATLAS Zero Degree Calorimeters (ZDC). 
The detector measures neutral particles with \rj\ > 8.3, in particular spectator neutrons in heavy 
ion collisions; all charged nuclear fragments are swept away by magnetic fields before reaching the 
ZDC. Each side contains four rectangular modules, the first of which is an electromagnetic module 
with finer readouts. Each module contains 11 plates perpendicular to the beam direction made 
from a combination of tungsten and stainless steel absorbers. Between the plates are 1.55 mm 
quartz strips. Hadronic showers of incident particles are initiated by the absorber. Cerenkov 
light from the showering particles is produced in the quartz and fed through air light-guides into 
photomultiplier tubes situated on top of the detector. Some of the modules contain specialized 
quartz rods to provide position measurements, although this aspect of the detector capability is 
not used in minimum bias event selection for heavy ion collisions. 
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For the 2010 ion run, the ZDC was used primarily as a minimum bias trigger. Each side is 
capable of producing a LI trigger signal by passing the analog channel sums through a discriminator 
to the CTP. The coincidence trigger, ZDC_A_C, is a logical AND of the one-sided triggers, ZDC_A and 
ZDC_C. The requirement of a single neutron on both sides is effective in rejecting against photo- 
nuclear collisions, which are typically asymmetric and are a large background to the most peripheral 
collisions. 



3.7 Inner Detector 

The ATLAS inner detector occupies a cylindrical volume around the detector center spanning 
±3512 mm in the z direction and 1150 mm in radius. Space-point measurements are provided 
by three separate sub-detectors over the radial extent of the system, each divided into barrel and 
end-cap modules. A slice of the detector in the R — z plane showing the various components is 



shown in Fig. 3.16 Additional figures showing perspective views of the inner detector system in 



the barrel and end-cap region are shown in Fig. 3.17| 



The silicon pixel detector is composed of three cylindrical layers beginning at a radial distance of 
50.5 mm and three end-cap disks per side. It contains 1744 identical pixel sensors which are 250 yum 
thick oxygenated n-type wafers with readout pixels on the n + -implanted side of the detector. Each 
sensor contains 47232 pixels for a total of approximately 80 million readout channels (the total 
number of channels is reduced by the ganging of adjacent pixels on front-end chips). The pixels 
have a nominal size of 50 fim in R — (ft and 400 /im in z providing measurements with intrinsic 
inaccuracies of 10 /jm in R — (ft for all modules and 115 /im in z and R for the barrel and disk 
modules respectively. 

The silicon microstrip detector (SCT) consists of four cylindrical layers and nine disks per side 
each composed of small angle stereo strips positioned beyond the pixel detector. The barrel covers 
the radial region 299 < R < 514 mm, and \z\ < 805 mm. The end-cap modules are positioned 
251 < R < 610 mm and 810 < \z\ < 2797 mm. The strips are composed of two 6.4 cm daisy-chained 
sensors, 285 ± 15 /*m thick, with a strip pitch of 80 /im. In the barrel these strips are oriented 
parallel to the beam direction, measuring R — eft, while in the end-caps they are trapezoidal wedges 
and are oriented radially. 
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Figure 3.16: Schematic view of the inner detector showing the positions of the various modules in 
R and z. 

The Transition Radiation Tracker (TRT) is the outermost component of the tracking system. 
The barrel spans radial region 554 < R < 1066 mm and \z\ < 780 mm. The end-cap wheels span 
615 < R < 1106 mm, positioned between 827 < \z\ < 2744 mm. It is composed of polymide 
drift tubes, 4 mm in diameter and 71.2 cm in length. The tube wall is composed of two 35 //m 
multi-layer films. The cathode is a 0.2 /xm Al film on the inner surface. The anode wires are 31 p 
gold-plated tungsten, positioned at the nominal center of the tube. To ensure stable operation 
the wires are required to have an offset with respect to the tube center of less than 300 fim. The 
tubes are filled with a gas mixture of approximately 70% Xe, 27% CO2 and 3% O2, held slightly 
over-pressure. Photons from transition radiation of high energy electrons can be distinguished from 
normal ionization tracking signals on a per tube basis using separate high and low thresholds on 
the readout electronics. The wires are divided into three sections. They are read out at each 
end, with an inefficient center section of approximately 2 cm, which contains a mid-wire support 
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Figure 3.17: Perspective view of the inner detector showing the barrel (left) and end-cap (right) 
regions. 

structure. The TRT barrel is composed of 73 layers of tubes oriented along the beam direction, 
while the end-caps consist of 160 tube planes. Charged tracks will cross at least 36 tubes except in 
the barrel/end-cap transition region (0.8 < \rj\ < 1.0) where this minimum drops to 22. 
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Chapter 4 



Jet Reconstruction 



4.1 Subtraction 

For a calorimetric measurement, the input to the jet reconstruction procedure is the measured E^ 
distribution, which contains both the jet signal and contributions from the underlying event, 



The magnitude of the underlying event can vary over several orders of magnitude depending on 
the collision geometry, thus an accurate procedure requires that a background estimation be per- 
formed on an event-by-event basis. Correlations in the underlying event, specifically t>2, must be 
incorporated into the subtraction procedure as well. The background energy density is estimated 



2n in azimuth in strips of constant r/ with width 0.1. To account for variation of detector response 
in the different longitudinal sampling layers, the background determination and subtraction is 
performed per layer. 

The possibility exists for the jets to bias the determination of p. If a jet's energy is included in 
the background calculation, the background is overestimated, resulting in an over-subtraction being 
applied to the jet. This effect, called the self-energy bias, will bias the jet energy after subtraction 
by approximately 10% and will distort the V2 as well. To address this, the jet reconstruction used 



dE^ 1 _ dE^ E dEp 



(4.1) 



drjdcf) drjdcf) drjdcf) 




is the average -E/p density taken over the full 



(4.2) 
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in this analysis employed an iterative procedure where jets were reconstructed and identified as jet 
seeds. B was constructed with the seeds excluded from the determination of p. This reduced the 
bias in the background and the subtracted jets. These jets were then be used as seeds to iterate 
the procedure, by constructing new estimates of p and vi- 

4.1.1 Seeds 

In clustering algorithms, such as anti-A^, the entire r\ — (f> distribution is tessellated by reconstructed 
jets. While some of these jets likely correspond to a true jet signal, the clustering of many of the 
background jets is due to the underlying event. If a set of criteria can be chosen to partition the 
full set of jets output by the reconstruction into one that characterizes the background distribution 
and one that does not, the latter can be used as jet seeds. A simple -E/p threshold cannot be used 
to define the seeds at this stage since the subtraction has not yet been performed. For the first step 
in the iterative procedure, a discriminant was constructed by considering the distribution of tower 
.E/r's inside a jet and taking the ratio of the maximum and mean of this distribution: 

max(£t° wer ) , N 

p= ( 4-> • (4 - 3) 

A cut of D > 4 identifies jets with a dense core spanning a few calorimeter towers, which are 
typical of the true signal. This requirement fails when calorimeter noise allows the denominator 
of the discriminant to become arbitrarily small resulting in a large discriminant even when the 
distribution maximum is also small. To prevent this from occurring, a threshold on D was be 
applied in conjunction with a minimum requirement on the tower distribution maximum or mean. 
This cut was taken to be a minimum value of max(^ owcr ) > 3 GeV. 

In subsequent stages of the iteration procedure the jets that were the output of the previous 
iteration step were reused as seeds. Since these jets had already been subtracted, an E^ cut of 
25 GeV could be applied to define the seeds. In addition to these jets, track jets with > 10 GeV 
were also added to the list of seeds. 
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4.1.2 Background Determination 

The average density was computed by considering all calorimeter cells and calculating an average 
transverse energy density per cell. For a sampling layer i the average is defined as 



pi(v) 



En 



-y 

Ni y> Ar]jA(j)j 1 + 2v 2 i cos(2(^- - * 2 )) 



(4.4) 



The average was only taken over cells that are in the same pseudorapidity bin as rj and are not 
associated with a seed. In the discriminant-based approach, the seed-associated cells were defined 
as the constituents of the jets that pass the discriminant cut. In all other cases they were defined 
as being within AR < 0.4 of a seed position as shown in Fig. |4.1| This exclusion radius was used 
regardless of the R value used in the anti-/ct algorithm. The ^-dependent demodulation factors were 
required since the average was not performed over the full (f) range and thus elliptic flow effects do 
not cancel. 



V 



V 



Figure 4.1: An example of how seeds are excluded from the p (left) and t>2 (right) calculations. The 
positions of the seeds are indicated by white dots. The red regions are excluded from the averages 
and the blue regions are included. 
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Figure 4.2: ^2 computed from the second layer of the electromagnetic barrel as a function of the 
leading jet's angle with respect to the event plane ^2- Two centrality bins (0-20% and 40-60%) are 
shown before (solid circle/triangle) and after (hollow circle/triangle) the iteration step. Without 
removing the seeds from the calculation, the i>2 is enhanced in events where the leading jet is aligned 
with the event plane and reduced when the leading jet is out of plane. The correlation is removed 
by the iteration step. 

4.1.3 Elliptic Flow Correction 

The event plane angle was determined using the first layer of the FCal using the same procedure 
as in Refs. |211j . This angle was then used to determine the V2 per calorimeter sampling layer i, 



v 2i = (COS 2(0- tf 2 ))t 



(4.5) 



where the brackets denote an Ei weighted average over all cells in layer i. 

As was the case for p, it is possible for the jet signal to bias the V2- This effect is worsened 
by the fact dijet signals have the same ir symmetry as the flow modulation. Jet seeds can be used 
to exclude regions containing potential jets from the averaging and eliminating the bias. However, 



excluding small regions around a seed as described in Section 4.1.2 cannot be done. In a strip 
of constant rj, excluding some of the cells from the full azimuthal interval will cause an anti-bias; 
the V2 for events where jets are aligned with the event plane will be reduced. To prevent this 
overcompensation, the entire 77 interval was considered to be biased by the jet. Such rj intervals 
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were flagged using the jet seeds and excluded from the average in Eq. 4.5 This exclusion is shown 



schematically in Fig. 4.1 The effect of jets on V2 as measured in the second sampling layer of the 



barrel electromagnetic calorimeter in MC simulations is shown in Fig. 4.2 Without excluding jets 
from the 1)2 determination (solid markers), the calculated 1)2 is larger in events where the leading 
jet is aligned with the event plane. This bias is removed once the jet seeds are excluded from the 
calculation as indicated by the open markers in the figure. 

4.1.4 Subtraction 

Once the background and flow parameters were determined, a subtraction was applied at the cell 
level. For a cell at position rj, in layer i, the corrected E^ was determined from the TclW Et vicl 

E T = Ef w - p i (r ? )A 7 ?A0[l + 2v 2i cos(2(0 - * 2 ))] . (4.6) 

Each cell was considered as a massless four-vector using the E^ after subtraction. The jet's kine- 
matics were set as the four-vector sum of the constituent cells. 

This subtraction was applied to all cells except dead cells that receive no correction from 
neighboring cells or the associated trigger tower. The latter correction was only applied if the 
trigger tower was above a 2 GeV noise threshold. Thus for the purposes of subtraction, dead cells 
eligible for this correction were checked event-by-event. 

4.1.5 Calibration 

Jet collections were constructed using two different calibration schemes. The EM+JES scheme uses 
the cell energies, after subtraction, at the EM scale and applies a final, multiplicative calibration 
factor to the jet energy to calibrate to the full hadronic scale. In the GCW scheme, a calibration is 
done at the cell level, which designed to minimize variation in hadronic response to jets. However, 
a multiplicative calibration factor is still required to account for the mean hadronic response. In 
this scheme the p determination and subtraction use calibrated cell energies^ The GCW also 



1 The V2 was determined using cells at the electromagnetic scale. For jets using the GCW scheme, the E™ w in 



Eq. 4.6 was first demodulated by the flow weight 1 + 2i>2; cos(2(</> — ^2)), calibrated and then re-weighted by the 



same factor. 
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accounts for energy lost due to material in the cryostat by calculating a correction using the energy 
deposition in the third electromagnetic and first hadronic sampling layers after subtraction. 

After subtraction the multiplicative jet energy calibration factor was applied to correct for the 
overall jet energy. The factors depend on both the jet's energy and rj, and different constants 
are used for the GCW and EM+JES schemes. The numerical inversion constants were developed 



specifically for heavy ion jets and are discussed in Section 4.3.2 



4.1.6 Track Jets 

Jets reconstructed with charged particles provide an invaluable cross-check on the calorimetric jets. 
The anti-/ct algorithm with R = 0.4 was run using charged particles as input. Instead of performing 
a subtraction, particles were required to have a minimum of 4 GeV to be included in a track 
jet in addition to the following track selection requirements, where impact parameter variables are 
defined with respect to the primary reconstructed vertex: 

• do, zqsiiiO < 1.0 mm, 

do z sin 6 

• , < 3.0, 

&d Cosine 

• AWi > 2, 

• A^sct > 8. 

4.1.7 Jet Collections 

With the exception of the track jets, all the results presented here use jet collections obtained 
with the iterative method using a single iteration step applied to both the background and flow. 
Seeds were constructed by running anti-fct with R = 0.2. The discriminant procedure was used to 
determine the background, and the initial V2 values were used with this background to produce 
subtracted, flow corrected jets. The list of seeds for the iteration step was constructed from jets 
with E>? > 25 GeV and all track jets with E^ > 10 GeV. This list of seeds was used to define new 
background and V2 values. These were used to construct iterated jet collections with four different 
values of R: R = 0.2, 0.3, 0.4 and 0.5. For each of these values of R, two versions of the collections 
exist corresponding to the two calibration schemes: GCW and EM+JES. The 25 GeV cut defining 
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the seeds was chosen based on the observation that the fake rate for R = 0.2 jets above this pp is 
not too high to cause biases in the subtraction when a fake jet is used as a seed. 



4.2 Monte Carlo Sample 



The studies presented here use a Monte Carlo sample of minimum bias HIJING Pb+Pb events with 
embedded PYTHIA pp events at ^/snn = 2.76 TeV. Events from both generators were produced 
and the detector response to each event was independently simulated using a full GEANT [2121 1213] 
description of the ATLAS detector |214j . The signals were combined during the digitization stage 
and then reconstructed as a combined event. For each HIJING sample a set of PYTHIA samples 
(referred to as J samples) was produced, each with a fixed range set on the p™ in and p™ ax in 
the PYTHIA hard scattering. A single HIJING event was overlaid on an event from each of the 
different J samples, with no reuse of HIJING events within a J sample, i.e. A^ JING = iV e J vt . The 
different J samples were then combined using a cross-section weighting obtained from PYTHIA to 
build a combined sample with good counting statistics over a wide range of jet px- The definitions 



of these samples and the associated cross sections are shown in Table 4.1 



J 


pf n [GeV] 


p^ x [GeV] 


a [nb] 


1 


17 


35 


157000 


2 


35 


70 


7090 


3 


70 


140 


258 


4 


140 


280 


5.85 


5 


280 


560 


0.0612 



Table 4.1: Definitions of PYTHIA pp samples used in embedding. For each, J value samples were 
produced with the same number of events to ensure high statistical sampling for a jet pp out to 
500 GeV. 



As HIJING does not contain a mechanism for simulating elliptic flow, an afterburner was applied 
to modulate the azimuthal angles of the HIJING particles using apT. 1 and centrality-dependent 
parametrization of existing V2 data |215j . Samples used in this study use a parametrization derived 
from 2010 Pb+Pb data at the LHC. 
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The goal of the embedding procedure is to study the combined effect of the underlying event and 
detector response on the jet signal. Truth jets were constructed from the generator-level particle 
content of the PYTHIA event only. These jets were then compared to the output of offline jet 
reconstruction to extract metrics of jet performance and information for correcting measurements 
for these effects. 

The HIJING events themselves contain jet production at a rate consistent with binary scaling. 
A per-event jet yield for these jets is shown in Fig. |4.3| A subtlety arises in the embedding procedure 
when the particles from a hard scattering in the HIJING background overlap with those from the 
PYTHIA jets. The reconstructed jet signal will contain contributions from both PYTHIA and 
HIJING, however the matching truth information will only contain the jet energy from PYTHIA. 
While this problem affects truth jets of all energies, it is most easily seen in samples of Jl events 
where the PYTHIA jets have relatively modest energies. In a high-statistics sample, some of the 
HIJING events contain jets in excess of 100 GeV, resulting in many events where these jets overlap 
in r)-<f) with a PYTHIA jet of much lower pp of order 10 GeV. Given the large weight of the Jl 
sample, without correcting for the HIJING contamination the most probable way to get a 100 GeV 
jet in the full study is through this mechanism. The associations between 10 GeV truth jets and 
100 GeV reconstructed jets significantly distort studies of jet performance and the response matrix. 

Running jet reconstruction on the HIJING particles leads to many of the challenges of back- 
ground separation/subtraction that are faced in the full offline reconstruction. Uncertainties in 
identifying the truth jet energy would reduce the effectiveness of any performance studies. It is 
possible for multiple hard scatterings in the same heavy ion event to produce jets that overlap in 
the calorimeter. This is likely to be a small effect and would be best addressed by extending the 
embedding process to multiple overlapping PYTHIA jets, where the identification of truth energy 
is clear. The procedure used in this study to remove HIJING contamination is to run the anti-/c t 
algorithm with R = 0.4 on HIJING particles with pr > 4 GeV. The energies of these jets were 
not used in the performance study, but instead any PYTHIA jet found within AR < 0.8 of a 
HIJING jet with > 10 GeV is removed from the sample. The fraction of jets removed during 
this procedure in each and centrality bin is recorded and the final distributions are re-weighted 
by the inverse of this fraction. 

The overlaid events must have the same vertex position at the generator level. This constraint 
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Figure 4.3: A per event jet yield as a function of E^ for jets produced with HIJING. Jets were 
reconstructed from a 75,000 event minimum bias sample running the anti-fot algorithm on particles 
with pt > 4 GeV. 



requires that multiple samples be produced at different, fixed z vertex positions, rather than the 
usual case of a single production sampling a continuous distribution. The positions of the samples, 



and the relative number of events produced per position is shown in Fig. 4.4 The vertex x and y 
positions were fixed for all samples to match the mean beam-spot position from the 2010 Pb+Pb 
run. A total of one million HIJING events were generated along with one million PYTHIA events 
in each of the five different J samples. 
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Figure 4.4: Distribution of 2 vertex position in combined MC sample corresponding to one million 
events for each J sample. Transverse position is fixed at x = 0.1352 mm, y = 1.1621mm. 

4.3 Corrections 



4.3.1 Self and Mutual Energy Biases 



As mentioned in Section 4.1 , a self-energy bias can occur when there is an incomplete match between 
the jets used as seeds and the final jets after iteration. To address this possibility, a variable was 
computed during reconstruction assuming the jet failed to be excluded from the background and 
thus biased its own energy subtraction. At the analysis stage, any jet not associated with a seed 
has its energy corrected. The iteration step reduces this problem considerably as all high energy 
jets (£t > 25 GeV) are excluded from the background and are unaffected. However for jets around 
this threshold the self-energy bias correction must be included. 

The above procedure assumes the bias is binary: either a jet is biased or it is not. However, if 
the region around a seed does not contain the entire geometric area spanned by the jet, a differential 
bias can occur determined by the region of non-overlap. Additionally, jets that are not excluded 
from the background may cause over-subtraction in another jet in the same 77 interval (mutual 
energy bias). The underlying event can produce calorimetric signals that appear as jets, but are 
in fact uncorrelated with hard particle production in the event. The contribution and systematic 



removal of such fake jets is discussed extensively in Section 4.3.3 The set of jets capable of biasing 
the background (or being biased themselves) must be consistent with the definition of fake jets. 
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The energy bias correction was applied by considering the constituent towers of jets that pass 
fake rejection and calculating a bias for each strip in r] by summing the energies of those towers 
unassociated with a seed (AR > 0.4). Each jet was then corrected at the constituent level for 
biases introduced by constituents of other jets that were not excluded by the seeds. 

4.3.2 Numerical Inversion 

Numerical inversion is a method to derive the 77 and /^-dependent jet calibration. The method 
described below follows the procedure used in the jet energy scale determination used in pp jet 
measurements |216j . The main motivation to derive calibration constants for heavy ion jet recon- 
struction is that no such constants exist at this time for R = 0.2, 0.3, and 0.5 in pp. Additionally, 
some aspects of the reconstruction differ slightly between pp and Pb+Pb, most notably the issue of 
noise suppression. Thus the application of numerical inversion constants derived using heavy ion 
jets will be more appropriate than the existing R = 0.4 and R = 0.6 pp constants. 

The calibration constants are derived using 60-80% peripheral events from PYTHIA+HIJING 



Monte Carlo samples described in Section 4.2 The basic quantity used in the derivation of the 
calibration constants is the response, 

* = ■ ( 4 - 7 ) 

.C/rp 

where -^Tclio * s reconstructed Et at the EM scale and £^ mth is the E^ of the truth jet. The response 
was evaluated for all calorimeter jets which match a truth jet and track jet simultaneously. The 
matching condition required a separation of AR < 0.3 between the calorimeter jet axis and truth 
(track) jet axis. The calorimeter jet was required to be isolated, that is no other calorimeter jet 
within AR = 2.5R must be present (R is the parameter of the anti-Ait algorithm). 

The response was evaluated in 12 bins in pseudorapidity of the reconstructed calorimeter jet, 
r/det) 011 the interval |??det| < 2.8 and 23 bins in the £^ mth . For £^ mth in bin j, the quantities (TZ) j 
and (^Tcaio)j were determined using a Gaussian fit and a statistical mean respectively. The set of 
points ((-Sxcalo)i' (^)i) was n * with the function: 

N max 

W^J = X> (ln(^ alo )) 4 , (4.8) 
i=0 

where N max is between 1 and 4. The actual value of N max is optimized to provide the best x 2 /NDF. 
A separate set of points and fit parameters were determined for each T^et- The calibrated E^, 
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I Cenlrality 60-80%, -0.8 < eta < -0.3 | 



I Cenlralily 60-80%, -0.3 < ela < 6~] 




I Centrality 60-80%, < ela < 0.3 | 





Figure 4.5: The i^xcilo dependence of the response with the logarithmic fit. 



£: ™a!o JES > was then defined by 



E; 



EM+JES 
Tcalo 



F EM 
^Tcalo 

-t'calibl^TcaloJ 
EM 



(4.9) 



An example of the response as a function of ^xcaio * s shown in Fig. 4.5, where the response 



is well described by the functional form of Eq. 4.8 A similar level of agreement is found for all 



r/det bins. Figure |4.6| shows the rj dependence of the jet energy scale after applying the derived 
calibration constants. A detailed closure test is discussed in Section 14.4.11 



4.3.3 Fake Rejection 

Fluctuations of the underlying event that are unassociated with a hard scattering can create lo- 
calized high E'y regions in the calorimeter that may be reconstructed as jets. In the most central 
collisions these fake jets contribute significantly to the jet spectrum up to 80 GeV. As most of this 
signal is caused by accumulated energy from soft particles, requiring a physics signature consistent 
with hard particle production in the neighborhood of a jet can reject against this background. 
To minimize the contamination of the sample of reconstructed jets by fake jets different methods 
for fake jet identification and rejection were studied. The general requirement of a coincidence of 
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Figure 4.6: r/ dependence of the jet energy scale after applying the derived calibration constants. 



the reconstructed calorimeter jet with the track jet or an electromagnetic cluster was found to be 
efficient at removing fakes without a severe impact on the overall efficiency. For the duration of 
this thesis, the fake rejection, unless otherwise noted, is defined as the requirement that a jet have 
a track jet or EM cluster with px > 7 GeV within AR < 0.2 of the jet's axis. The effectiveness of 



this rejection procedure is discussed in Section 4.3.4 



Fig. 4.7 compares the jet reconstruction efficiency in central (0-10%) and peripheral (60-ou7oj 
collisions in two cases: before applying the fake rejection (solid symbols) and after applying the fake 
rejection (open symbols). The overall jet reconstruction efficiency is reduced by the fake rejection 
requirement but still remains near 50% for 30 GeV jets. 

Corrections derived from MC for efficiency loss due to fake rejection (moving from solid to open 



markers in Fig. 4.7) will be sensitive to how accurately the MC describes this effect. The accuracy 
of the MC result can be checked in the data by comparing the ratio of the spectrum after fake 
rejection to the total spectrum, which shows the survival fraction of jets at a given that pass 
the rejection criteria. This comparison must be made with caution and only in regions where a 
significant fake rate in the data is not expected because the embedded MC sample is not expected 
to accurately describe the fake rate. Since every event in the embedded MC sample contains a 
PYTHIA hard scattering event the relative rates of true jets (from the PYTHIA signal) to fake 
jets (from soft particles in the HIJING background) will not be correct. A HIJING-only sample is 



required for an evaluation of the absolute fake rate, a procedure which is discussed in Section 4.3.4 
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Figure 4.7: The jet reconstruction efficiency in central (black) and peripheral (red) collisions both 
before (solid) and after (open) fake rejection for R = 0.4 jets. 



A comparison of the survival fraction^] in data and MC is shown in Fig. 4.8 The fake rate is 
expected to be negligible in peripheral collisions and good agreement is seen between data and MC 
in the survival fraction in the 60-80% centrality bin (red markers). This indicates that efficiency 
loss of true jets due to the fake rejection requirement is well described by the MC, and that a 
MC-based correction for this effect is appropriate. For > 100 GeV, fakes are also expected to 
be a negligible contribution to the spectrum in the 0-10% bin. In this region, where the efficiency 
approaches its asymptotic value, the agreement between data and MC is better than 0.5%. 

High-px tracks and clusters are expected to be associated with particles within jets, however 
it is possible for reconstructed jets to pass fake jet rejection due to combinatoric overlap between 
the jet and the track jet/cluster. This combinatoric rate is the mechanism by which fake jets pass 
rejection, and a data to MC comparison of this rate can be used to construct a more accurate 
estimate of the fake rate in data from the residual fake rate obtained from the MC study. This 



4.7 



One 



2 It should be emphasized that this distribution is not trivially related to the efficiency shown in Fig. 
important difference is that the efficiency is shown as a function of pt of the truth jets. The survival fraction shown 
in Fig. |4.8| is evaluated as a function of reconstructed px, since the data is being used, which has the effects of upward 
feeding from finite resolution. 
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Figure 4.8: The effect of the fake rejection data and MC in the 0-10% and 60-80% centrality bins 
as defined by survival fraction, the ratio of the spectra after and before fake rejection. 

rate was evaluated by considering the AR distribution of all track jets, electrons and photons in an 
event with respect to a jet's axis. To simplify geometric effects from splitting of nearby jets, only 
isolated jets were considered. These jets were required to be the highest pt jet in a AR < 1 cone 
about their axes. In the MC sample jets were additionally forced to match a truth jet (AR < 0.2) 
and be isolated from any HIJING jets (AR > 0.8 for £^ IJING > 10 GeV). Comparisons between 
the different objects (track jets, electrons and photons) and between data and MC are shown 



for R = 0.4 jets in central and peripheral collisions in Fig. 4.9 The distributions show strong 
correlation for matches at small AR with a more diffuse peak at AR ~ ir associated with the dijet 
structure of the events. The remainder of the distribution is expected to be determined purely by 
geometry, thus the trivial geometric enhancement has been removed by dividing by R. With this 
effect removed, the matching rate is expected to be constant away from the jet and dijet peaks, 
which is estimated by the background extrapolations shown in the figures. These distributions have 
been normalized per jet, such that the background distribution integrated over the matching region 
can be interpreted as a per jet probability of combinatoric match, P, via 

P= F(— — — W<iAi^0.02x (J-J- ™L) , (4 . 10 ) 

where the quantity in parentheses is estimated from fitting the distribution in the region 1.5 < AR < 
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Figure 4.9: Data to MC comparison of AR distributions between track jets, electrons and pho- 
tons and R = 0.4 jets in central (top) and peripheral (bottom) collisions. Distributions have been 
normalized per jet and have radial Jacobian removed to facilitate extraction of combinatoric back- 
grounds (dashed lines). 



2 and extrapolating underneath the correlation peak at AR < 0. Estimates of these probabilities in 



central collisions are compared for data and MC in Fig. 4.10| The rates are nominally independent 
of R, which is expected as the jet definition does not explicitly enter into the matching. The MC 
are higher than the data by factor of 2-3, which is qualitatively consistent with observed single 
particle suppression factors, caused by a phenomenon (quenching) present in the data but not the 
MC. 
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Figure 4.10: Comparison between combinatoric matching probabilities between data and MC for 
different R values. Rates are generally independent of R and are a factor of 2-3 larger in MC, 
consistent with the absence of quenching. 



4.3.4 Residual Fake Rate 

The residual per-event rate of fake jets (underlying event fluctuations) after fake rejection was 
determined from an analysis of a HIJING-only MC sample. Each reconstructed jet passing the fake 
rejection criteria but not within AR < 0.2 of a HIJING jet with E^ uth > 10 GeV was considered a 
fake jet. The residual fake rate is defined as the ratio of these unmatched jets to the total number 
of reconstructed jets passing fake rejection and is a function of E^ co . 

Unlike the truth jet reconstruction in PYTHIA-embedded events, the HIJING truth jet recon- 
struction was run with a > 4 GeV cut on the input particles to the reconstruction algorithm. 
Because of this, a hard-scattered parton may fragment in a way that is not reconstructed by the 
HIJING truth jet reconstruction clS Si jet with E^ uth > 10 GeV. A jet reconstructed fr om these 
fragments will fail to match a truth jet causing the fake rate calculation to come out anomalously 
high. Therefore, it is necessary to calculate the truth jet reconstruction efficiency of a truth jet if 
the reconstruction had been performed with a > 4 GeV cut on the constituents. After correcting 
the fraction of jets that match to a HIJING truth jet by the truth jet reconstruction efficiency, the 
difference between the efficiency corrected rate and the full jet rate is the residual fake rate. 

Only reconstructed jets within < 2.0 are considered for the residual fake rate analysis. 
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Reconstructed jets that pass fake jet rejection are corrected for the self-energy bias as described in 



Section 4.3.1 The analysis is repeated for three different fake jet rejection schemes: 

• pt > 7 GeV track jet or E T > 7 GeV cluster within AR < 0.2, 

• pt > 7 GeV track jet or E T > 9 GeV cluster within AR < 0.2, 

• pt > 7 GeV track jet within AR < 0.2. 

To calculate the truth jet reconstruction efficiency resulting from the pi > 4 GeV constituent 
cut, a subset of the Jl and J2 PYTHIA-embedded Monte Carlo samples were used (400, 000 events 
z = 2.87 mm only). PYTHIA truth jets were matched to reconstructed jets with AR < 0.2. 
To determine if this truth jet would be reconstructed by the definition used in the HIJING only 
sample, the summed energy of PYTHIA truth particles with pi > 4 GeV within a AR < 0.4 
radius of the truth jet axis was required to be greater than 10 GeV. Because > 4 GeV truth jets in 
the underlying HIJING event can combinatorially overlap with a PYTHIA truth jet, the truth jet 
reconstruction efficiency is actually higher for a PYTHIA truth jet at the same £^ mth in central 
than in peripheral collisions. To include this effect, the E^ of HIJING jets within AR < 0.4 of the 
truth jet's associated reconstructed jet's axis was also added the PYTHIA truth particle sum. 

As previously noted, a significantly higher- E^ hard scattering in HIJING may combinatorially 
overlap with an embedded PYTHIA truth jet. To prevent such scenarios from contaminating the 
efficiency calculation, the same HIJING contamination removal condition is applied. Even with 
this constraint, there are cases when a HIJING jet overlaps a PYTHIA jet. However, this effect 
of such overlaps is not obvious due to truth jet reconstruction inefficiency. Therefore, jets with 
E v ^ co - E^ uth > 30 GeV (about ~ 2-2.5 times the RMS energy from underlying event fluctuations 
in central events for R = 0.4 jets) were also excluded. In central events and peripheral events, 
~ 75% and ~ 99.5% of jets, respectively, are properly isolated from contamination from HIJING 
hard scatterings. 

To determine the efficiency in each E^ co bin, the J samples were combined in a weighted average, 

^ej(I? T ™)^a J N$(E™°) 
e(E^°) = ^ t ^— , (4.11) 

j iv evt 
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Figure 4.11: Truth jet reconstruction efficiency as a function of E^ co in 0-10% events (left) and 
60-80% events (right). The efficiency using three different fake jet rejection schemes is shown. 

where N^ vt is the number of generator events in that centrality selection, a J is the generator cross- 
section for the J sample given in Table 4.1 and N^ et (E^ co ) is the number of jets at the given 



reconstructed E^ co . Figure 4.11 shows the truth jet reconstruction efficiency in the most central 
and most peripheral events for each of the three fake jet rejection schemes. 

The truth jet reconstruction efficiency was then used to determine the residual fake rate in 
central events from the HIJING-only data set. Reconstructed jets that pass fake rejection were 
tested to see if they match to a HIJING hard scattering. This is accomplished by considering the 
summed E^ of all HIJING truth jets with E^ > 4 GeV within AR < 0.4 of the reconstructed jet. 
Reconstructed jets where this sum exceeds 10 GeV were considered a match, with the unmatched 
yield defined as the remainder after the matched is subtracted from the total. The total per-event 
yield in the 0-10% centrality bin along with the matched and unmatched is shown on the left side 



of Fig. |4.12 The matched distribution was then corrected for the truth jet inefficiency and the 



difference between the total yield and this corrected, matched yield provides an estimate of the 



absolute fake rate is shown on the right side of Fig. 4.12| 

For reconstructed jet at a given p^, the probability that this jet is unmatched defines the 
residual fake rate. This quantity was calculated by normalizing corrected, unmatched yield by the 



total yield, and is shown in Fig. 4.13 for the three fake jet rejection schemes in 0-10% centrality 



bin. For E^ co > 40 GeV, this rate is 4-5% for the fake jet rejection scheme used in this analysis: 
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Figure 4.12: Per-event yields of reconstructed in 0-10% central events using the HIJING-only 
sample. The total yield (black) is the same in both figures. The unmatched contribution is what 
remains after subtracting the matched from the total. The left plot shows the matched (blue) 
and unmatched (red) as extracted from the MC. The right figure shows the matched (green) after 
correction for truth jet reconstruction inefficiency and the residual fake yield (purple). 

matching to track jet or cluster with > 7 GeV. This rate drops to 1-2% for the schemes with 
stricter requirements, although these were not used in the jet spectrum analysis due to their effect 
on the efficiency. If a fake jet passes the rejection it is due to a combinatoric overlap between 
an underlying event fluctuation reconstructed as a jet and a single track or cluster also from the 
underlying event. The rate for this behavior was estimated in the previous section, and the MC 
rate was found to exceed that of the data by a factor of 2-3. The residual fake rate extracted from 
the MC should be reduced by this factor when describing data. Thus the absolute fake rate is 
estimated to be 2-3% for E^ > 40 GeV in the most central 10% of events. 



4.4 Performance 

To evaluate jet performance using the Monte Carlo samples, "truth" jets were matched to re- 
constructed jets by requiring an angular separation between truth and reconstructed jet, AR = 
\J At] 2 + A(p 2 < 0.2. If multiple reconstructed jets were matched to a truth jet, the reconstructed 
jet with the smallest AR was selected. For the results presented in this section, performance plots 
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Figure 4.13: Residual fake rate in 0-10% events (as a fraction of the total reconstructed jet spectrum) 
shown for three different fake jet rejection schemes. 



are presented for all jets with < 2.8 using the EM+JES calibration scheme. 

4.4.1 Jet Energy Scale 

For each matched truth jet, the E^ difference, 



AE T = E T rS co 



iritruth 



(4.12) 



was calculated. The distribution of AE^ was then evaluated as functions of E^ uth , ^ tmth anc [ 
centrality. To evaluate the jet energy scale (JES) as a function of E^ mth the AE^ distributions 
were each fit by a Gaussian. At low E^ uth ( E^ uth < 50 GeV for R = 0.4 jets) the AE T distributions 
are affected by a truncation of reconstructed jet Et at 10 GeV. For finite bin size in E^ uth , this 
results in a distortion of the shape of the AE^ distribution at low E^ nth . Therefore the AEt 
distribution was fit in the region which is unaffected by this truncation. Above 50 GeV the effect 
of the truncation is minimal as indicated in the left plot of Fig. ??. 

Upper plots and 



The evaluation of the JES as a function of E^i uth is shown on Fig. 



4.14 



lower two plots of Fig. 4.14 compare the JES for different jet definitions in the 0-10% and 6O-01170 
centrality bins. Lower right plot then shows the JES in seven centrality bins for R = 0.4 jets. 
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At low £^ uth the accuracy of the evaluation of the JES is limited by precision of fitting of the 
truncated Gaussian distribution. In these cases, to the left of vertical lines are provided in the 
figures, an upper limit on the JES (as indicated by the arrows on the error bars) is estimated. The 
difference between central and peripheral collisions is at most 2% for R = 0.2 — 0.5 collections. The 
one percent non-closure in R = 0.2 jets is due to the fact that the derivation of numerical inversion 
constants requires isolated jets, but the JES evaluation includes non-isolated jets as well. 
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Figure 4.14: Jet energy scale vs E^ uth for different jet definitions compared between central and 
peripheral collisions. 



Comparison of JES as a function of rj tmth is shown in Fig. 



4.15 



The JES as a function of 



n 



truth 



was evaluated by calculating the mean of the AE't /Et distribution of all reconstructed jets 



corresponding to truth jets with E^ > 90 GeV. 

The behavior of the JES as a function of E^ uth can be summarized by fitting the distribution 
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Figure 4.15: Jet energy scale vs ^ truth for different jet definitions compared between central and 
peripheral collisions. 
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Figure 4.16: Non-closure in JES calibration obtained from fits. 



with a linear function: 

AE T /E T (E^ uth ) = aE^ uth + b . (4.13) 
The a term is typically small, and the constant b term represents the non-closure in the JES 



calibration. The values of this constant term for all centralities and R values are shown in Fig. 4.16 
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4.4.2 Jet Energy Resolution 



The jet energy resolution (JER) can be denned using the width of the Gaussian fit discussed in 
the previous section. Figure |4.17| shows the evaluation of JER in different bins of centrality and 
for different jet definitions. The JER increases both with increasing centrality and increasing jet 
radius. Validation studies of the JER using data are presented in Sec. |5.3.2 
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Figure 4.17: Jet energy resolution vs E^i uth for different jet definitions compared between central 
and peripheral collisions. 

An alternative definition of the JER uses the standard deviation of the AEt/E^ distribution 
in bins of E^ uth ,i.e. 



a(AE T /E T ) = 




AE T 

iritruth 



(4.14) 



Figure 4.18 shows the evaluation of the difference between the two JER definitions for R = 0.4 jets 
and 0-10% central collisions. The Gaussian-fit JER is slightly smaller for E^ uth < 50 GeV, which 
can be attributed to the asymmetric AE^ distribution resulting from the minimum E^ threshold 
on reconstructed jets. 
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Figure 4.18: Comparison of JER evaluated using the Gaussian fit and as a standard deviation of 
the AE^/En? distribution. 

4.4.3 Jet Reconstruction Efficiency 

The jet reconstruction efficiency is defined as: 



e(E% uth ) 



/\ ^ymatch 



(4.15) 



' AATtruth ' 

where AiV truth and AiV match are the total number of truth jets and number of truth jets with 
a matching reconstructed jet in a given E^ uth bin. The reconstruction efficiencies are shown in 



Fig 4.19 The R = 0.4 and R = 0.5 jet reconstruction efficiency reaches 95% by 60 GeV in the 
most central (0-10%) Pb+Pb collisions and 30 GeV in the most peripheral (60-80%) collisions. 
The last 5% increase in the efficiency takes place near 70 GeV for all centrality bins. For the 
R = 0.2 and R = 0.3 jets, the efficiency reaches 95% below 40 GeV for all collision centralities and 
is approximately 1 above 50 GeV. 

Figure 4.20 shows the jet reconstruction efficiency a function Q f ^ tmth for jets with Ey > 40 GeV. 



The smaller radii have saturated the efficiency turn on at this energy and show no r/ variation. The 
R = 0.4 and R = 0.5 jets are not full efficient at this energy in the most central collisions and show 
a slight reduction in the efficiency between the barrel (\r]\ < 1.5) and end-cap (1.5 < \r]\ < 3.2) 



regions. 



CHAPTER 4. JET RECONSTRUCTION 



105 



Efficiency |f/| < 2.1 

• R=0.2 

• R=0.3 











1 E T [GeV; 






10 2 r 


[GeV] 










E T [GeV 




1 E T [GeV 


Elficienc 






















-•- 






































Pff ir*iorw»v Inl ^91 


































■ R=0.4 










0-10% 








10-20% 












20 - 30 % 


































• R=0.5 


0.4 
0.3 

1 1 
| 0.9 
S 0.8 






















































































































0.7 
0.6 








































; 30 - 40 % 








; ; 40 - 50 % 












50 - 60 % 




60 -80 % '"""] 


0.5 
0.4 


































































0.3 



































































1 E T [GeV] 



1 E T [GeV] 



Figure 4.19: Jet reconstruction efficiency vs E^ uth in the EM+JES scheme for jets with |7y truth | < 2.1 
for all R values and centralities. 
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Figure 4.20: Jet reconstruction efficiency vs r/* in the EM+JES scheme for jets with E^ > 
40 GeV for all R values and centralities. 
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Chapter 5 

Data Analysis 

5.1 Data Set 

The data presented in the following section were recorded by ATLAS during the 2010 Pb+Pb 
run. During this run a total integrated luminosity of approximately C = 7/ib -1 was recorded and 
determined to be of sufficient quality for physics analysis. An additional sample of C ~ l^ib^ 1 was 
recorded with the solenoidal magnetic field turned off; this data set will not be considered here. 

The events analyzed were required to satisfy a set of minimum bias event selection criteria. 
These were chosen to require a physics signature consistent with an inelastic Pb+Pb collision. 
Additional conditions were imposed to reject photo-nuclear and non-collision background. The 
requirements were as follows: 

1. The event was recorded because it satisfied a minimum bias trigger. Events were required to 
have fired one of the following triggers after prescale and veto (events dropped due to dead 
time): 

• MBTS coincidence: L1_MBTS_N_N, where N = 1, 2, 3 or 4. 

• ZDC coincidence: L1_ZDC_A_C or L1_ZDC_AND. 

2. Good MBTS timing. The A and C side times were required to satisfy |AtMBTs| < 3 ns. 

3. ZDC coincidence. Regardless of how the trigger selection was made the event must have 
one of the ZDC coincidence triggers before prescale, a logical OR between L1_ZDC_A_C and 
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L1_ZCD_AND. This was found to be essential for rejecting against photo-nuclear events. 
4. A good reconstructed vertex. 

Applying these event selection criteria, combined with selecting luminosity blocks consistent with 
stable running and detector operation resulted in a sample of 50 million events. These event 
selection criteria were used in all heavy ion analyses within ATLAS using the 2010 data. 



5.2 Centrality Determination 



The centrality determination used the total transverse energy as measured by the ATLAS FCal 
system, EE™ . The pseudorapidity coverage of this system, 3.2 < \rj\ < 4.9, is well separated from 
the central region of the detector, meaning that specific, centrality-dependent physics processes will 
not bias the centrality determination. The Si?^ Cal was found to be strongly correlated with the 
total energy deposited in the rest of the calorimeter, making it an excellent indicator of global event 
activity. The correlation between total energy in the electromagnetic barrel, \rj\ < 2.8, and TiE^P 8,1 



is shown in Fig. 5.1 
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Figure 5.1: Minimum bias EM barrel vs TiEJtP 8,1 correlation. 



The minimum bias S£'S Cal distribution has a distribution typical of other centrality variables 
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such as the total number of charged particles, as well as Glauber quantities N co \\ and iVp ar t. This 



distribution is shown in Fig. 5.2 A two-component model as given by Eq. 2.67 was used to fit this 



ATLAS 

2.76 TeV_ 



L int = 200 mb" 1 




FCal EE T [TeV] 

Figure 5.2: Minimum bias E£^ Cal distribution, with the labeled centrality bins corresponding to 
10% of the total integral of the distribution. 



distribution. An additional check was performed by taking an n-fold convolution of the pp TiE^P 3,1 
distribution based on a Glauber MC, which agreed well. A complete analysis of the 2010 centrality 
determination is given in Ref. |217| . The minimum bias event selection cuts were determined to 
introduce a 2% inefficiency, almost entirely from the MBTS timing cut. The distribution of A£mbts 
clS cl function of TjE^P^ is shown in Fig. 5.3l 

Under the assumption of 98% efficiency, the E£J^ Cal distribution was divided into 10% intervals, 
except for the range 60-100%, which was divided into two intervals of 20% each, in the fashion 



discussed in Section 2.5.1 In each bin the corresponding events in the Glauber MC were averaged 



to obtain iV co u and Apart- These values are shown in Table 5.1 This table also shows the ratios 
of Aeon factors, R^i which are used directly in the Rqp determination. The uncertainties were 
evaluated by considering the effects of varying ingredients into the Glauber calculation such as the 
nucleon-nucleon cross section, the two Woods-Saxon parameters as well as a 2% variation of the 
inefficiency. The errors on A co n in different centrality bins are correlated so the uncertainties on 
i^ccJ? were constructed appropriately. 
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Figure 5.3: A£mbts distributions in bins of T,E^ is shown on the left. The structure of the 
satellite peaks is clearly visible, as is the substantial broadening of the timing width at low £i?^ Cal . 
The Gaussian sigma of these distributions as a function of Si?^ Cal is shown on the right. 
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Table 5.1: Centrality bins, iV co u, -R^Ji vames an d their fractional error evaluated directly from the 
Glauber Monte Carlo [217] . 

5.3 Validation 
5.3.1 Jet Energy Scale 

A number of studies were performed using data to check the accuracy of the jet energy scale 
calibration procedure which was derived from MC. 
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5.3.1.1 EM+JES vs GCW Comparison 

The jet reconstruction was performed with both the EM+JES and GCW calibration schemes. A 



comparison of spectra using each of these schemes is shown for R = 0.4 jets in Fig. 5.4 In general 
good agreement is found between both schemes. 
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Figure 5.4: spectra for EM+JES (black) and GCW (red) calibration schemes for R = 0.4 jets and 
the ratio in central (left) and peripheral (right) collisions. Error bars in the ratio are constructed 
by considering which spectrum has the largest relative error and assuming the ratio has the same 
relative error. 



5.3.1.2 Track Jet Matching 

One of the methods to validate the performance of the jet finding, in particular the jet energy 
scale, is to match the calorimeter jets to track jets and compare their momenta. Such a comparison 
can be done both in the MC and in the data. Calorimeter jets were matched to track jets with 
Pt > 7 GeV, with a matching requirement of AR < 0.2 between the calorimeter and track jet axes. 
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Fig. 5.5 shows the calorimeter jet pp as a function of a corresponding track jet pp both for the 
data and MC for the 0-10% and 60-80% centralities. For calorimeter jets with pp > 50 GeV, the 
mean calorimeter jet pp, (pp al °), was computed as a function of the track jet ppand is overlaid on 
top of the two-dimensional distributions. 
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Figure 5.5: Two-dimensional distribution of p<r of reconstructed calorimeter jets versus pt of corre- 
sponding track jets for 0-10% central (left) and 60-80% peripheral (right) events for anti-&t R = 0.4 
jets. Upper plots show MC, lower plots show data. Markers indicate the average reconstructed 
calorimeter jet pr as a function of the track jet for calorimeter jets with pp > 50 GeV. 



The centrality dependence of this energy scale was assessed directly by constructing ratios of 



(p^ l °) distributions in central and peripheral collisions. Figure 5.6 compares this ratio as computed 
in data and MC. The ratios are roughly constant for track jets with p^ > 50 GeV, away from the 
minimum calorimeter jet pp threshold, and each distribution is fit to a constant in this region. The 
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values and errors of these fit parameters are shown in Fig. 5.7 for different centrality bins. The 
differences in the fit constants are evaluated directly in Fig. |5.8| The comparison indicates that the 
centrality dependence of the calorimetric response differs by no more than 2% between data and MC. 
The jet energy scale calibration is sensitive to the particle composition and fragmentation of the 
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.6: The ratio {p™ lo ) CC nt/ '(pf l °) mso as a function of p^ ack for both data (black) and MC 
R = 0.4 jets. Only jets with p^ al ° > 50 GeV are included in the average. A fit to a constant 
for p!p ck > 50 GeV. 



jet. In particular it was found in pp events that the response was lower for broader jets |216j . One 
concern is that quenching effects could result in broader jets introducing a centrality dependence to 
the JES: a lower response in central collisions relative to peripheral. However, such an effect is not 
consistent with the results of the in-situ study as there is no evidence for the (Px alo ) C ent/(PT al °)60-80 
to be systematically below one. 



CHAPTER 5. DATA ANALYSIS 



113 



1.06 
1.04 
1.02 
1 

0.98 
0.96 
0.94 
0.92 



R = 0.2 

• Data 



♦ MC 



■10% 10-20% 20-30% 30 - 40 % 40-50% 50-60% 



1.1r — 
1.08 E- 
1.06 z- 
1.04 
1.02 E- 
1 E- 
0.98 E- 
0.96 E- 
0.94 E- 
0.92 E- 

0.9E 



B = 0.3 
• Data 



♦ MC 



•10% 10-20% 20-30% 



- 40 % 40 ■ 50 % 





















1.08 














o 


1.06 














LL 


1.04 
















1.02 
1 




}. }• j. 


!• 










0.98 








f 








0.96 




R = 0.4 












0.94 
0.92 
0.9 




• Data ♦ MC 















■10% 10-20% 20-30% 


30 ■ 40 % 


40 ■ 50 % 


50 ■ 60 % 



1.08 E- 
1.06 E- 
1.04 E- 
1.02 E- 

1 E- 

0.98 E- 
0.96 E- 
0.94 E- 
0.92 E- 
0-9 o 



R = 0.5 
• Data 



♦ MC 



Figure 5.7: Fit constants for (p^ lo } C cnt/ (p'r lo }m~80 f° r data (black) and MC (red) in different 
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Figure 5.9: Comparison of the per-event standard deviation of summed -E/p for 7x7 groups of 
towers between data and the HIJING+GEANT MC simulated events as a function of FCal 



5.3.2 Jet Energy Resolution 
5.3.2.1 Fluctuations 

To evaluate the accuracy of HIJING in describing the underlying event fluctuations seen in the 
data, the per-event standard deviation of the summed E^ for all N x M groups of towers has 
been evaluated. The numbers M and N were chosen so that the area of the group of towers had 
approximately the same nominal area as jets reconstructed with a given distance parameter R, and 



values are given in Table 5.2 A complete systematic study is presented Ref. |218j : examples of 
that study are shown here. Figure |5.9| shows the event-averaged standard deviation as a function 
of the EEf Cal in the event for 7x7 groups of towers, which have approximately the same area as 
R = 0.4 jets. The same quantity is shown in both data and MC, with the MC EE™ scaled up 
by 12.6% to match the energy scale of the measured distribution. This procedure is equivalent to 
matching data to MC using fractions of the minimum-bias cross-section. The re-scaled MC results 



in Fig. 5.9 agree reasonably well with the results from the data over the entireX££ Cal range. Based 



on this data-driven test of the MC, it can be concluded that the HIJING+GEANT simulations of 
underlying event fluctuations can differ from the data by at most 10%. 
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9x9 



Table 5.2: Correspondence between the size of a jet and the size of groups of towers. 



5.3.2.2 Jet Energy Resolution Validation 

The knowledge of the size of the underlying event fluctuations can be used to validate the internal 
consistency of the jet energy resolution ( JER) estimates. The jet energy resolution can be expressed 
in the following form 



a(E T )/E T = a/y/E T © b / E T © c, 



(5.1) 



referred to as the stochastic, noise and constant terms respectively. The stochastic term represents 
fluctuations in the calorimeter shower sampling while the constant term relates to effects with 
AEt ~ Et such as dead regions. For the intrinsic energy resolution a(E^), the noise does not 
depend on the energy of incident particles and scales as l/E^ for the relative energy resolution 
ct(£'t)/-E'T- The constant b that quantifies the size of this term is the standard deviation of the 



noise energy and was evaluated in the fluctuation analysis described in Sec. 5.3.2.1. The JER was 



fit with the functional form of Eq. 5.1 with b taken from the fluctuation analysis and a and c 



determined by the fit, and a sample of these fits is shown in Fig. 5.11 The b values were adjusted 
slightly to account for the small differences in area between the rectangular tower groups and the 
nominally circular jet size, with area nR 2 . The fit parameters a and c are found to be independent 
of centrality, confirming that the centrality dependence of the JER is related to the underlying 
event fluctuations in a quantifiable way. The noise term extracted from the fluctuation analysis is 



shown in the left panel of Fig. 5.10 as a function of the jet area for different centralities, along with 
a fit =p APi. 



1 Although not directly relevant to jet performance, for all the centralities a fit value of p\ fa 0.59 was obtained. 
This hints at some structure in the fluctuations that should be investigated 
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jet area Centrality [%] 



Figure 5.10: The noise term obtained from the fluctuations analysis as a function of the jet area 
together with a power-law fit, b(A) = PqA Pi is shown on the left. The b values for fixed R value as 
a function of centrality are shown on the right. 
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Figure 5.11: The jet energy resolution in central and peripheral collisions for different jet sizes. Jet 



energy resolution is fit with the fixed noise term obtained from the fluctuation analysis. 
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5.3.3 Data Overlay 

To supplement the data-driven checks discussed in the preceding sections, samples using MC jets 
from PYTHIA overlaid onto Pb+Pb data were also studied. These samples were generated by 
combining minimum bias Pb+Pb events with PYTHIA dijet samples using the pileup overlay 
framework |219j and performing reconstruction on the combined signal. The PYTHIA events are 



produced in different J samples as discussed in Section 4.2 For each event the PYTHIA generation 
and subsequent GEANT simulation is run with conditions matching data, including vertex position. 
The data events are taken from a dedicated MinBiasOverlay stream in the 2011 running which 
uses a LI ZDC coincidence trigger and records the data without zero suppression. 

The validation of these samples is not yet complete, and high statistics samples have not yet 
been produced, however aspects of the ongoing validation studies provide additional support of 
the data-driven checks presented in the preceding sections. The overlay samples used embedded 
J3 and J4 PYTHIA jets. These results are only meaningful over the jet E^ range for which 
the J3 and J4 samples are the dominant contribution in the total cross section weighted sample 
(80 < Et < 300 GeV). The JES as a function of £t for R = 0.4 jets is compared between 0- 



10% and 60-80% centrality bins in Fig. 5.12 This comparison indicates no dependence of the 
JES non-closure on centrality to the level of statistical precision allowed by the sample (~ 1%). 
It is consistent with the estimates of the relative centrality-dependent JES non-closure derived 



in Section 4.4.1, and summarized in Fig. 4.16 Comparisons were made of different performance 
variables between these overlay samples and the J3 and J4 embedded PYTHIA+HIJING samples 



in described in Section 4.2 A comparison of the JER between these two samples in central and 
peripheral collisions is shown in Fig. 5,13| Excellent agreement is found between the overlay and 
HIJING samples, consistent with the data/MC agreement found in the fluctuations analysis |218j. 



A similar comparison of the jet reconstruction efficiency is shown in Fig. 5.14 Above 50 GeV the 
peripheral efficiencies show good agreement, and there is a slightly greater efficiency (~ 5%) in the 
HIJING sample in central collisions than in the overlay sample. As the Jl and J2 samples dominate 
at lower p^, the larger disagreement there is indicative of significant differences between the data 
and MC. 
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Figure 5.12: JES in the overlay sample for R = 0.4 jets in the 0-10% (black) and 60-80% (red) 
centrality bins. No systematic deviation is observed in the non-closure between the two centralities. 

5.3.4 Jet Kinematics 

As an evaluation of the jet in terms of pt and y is sometimes preferred to Ei and r]. For high-E^r 
jets, it is expected that px and Ei will be very similar as will r\ and y. This was verified directly 
in data. The two-dimensional correlations in E^ — Pt and tj — y in different centrality bins were 
computed in fine bins and are highly diagonal. The means of these distributions, (pp) vs and (y) 
vs rj, are compared in central and peripheral collisions for R = 0.4 jets in Fig. |5.15| No significant 
deviation from the behavior p^ ~ E^ and y ~ rj is observed. 

5.4 Unfolding 
5.4.1 Observables 

The goal of this measurement is to study jet quenching by examining the centrality dependence 
of the single inclusive jet spectrum. Collisions with different centralities have different degrees of 
geometric overlap between the colliding nuclei. This geometric effect translates into different per 
collision effective nucleon luminosities. Thus comparisons between different centrality bins should 
be scaled by iV co ii, the effective number of colliding nucleon pairs in a given centrality bin, to remove 
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Figure 5.13: Comparison of the JER in HIJING (solid) and the overlay sample (hollow) for R = 0.4 
jets in two centrality bins, 0-10% (black) and 60-80% (red). 
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Figure 5.14: Reconstruction efficiency in the overlay sample for R = 0.4 jets in the 0-10% (black) 
and 60-80% (red) centrality bins. No systematic deviation is observed in the non-closure between 
the two centralities. 



CHAPTER 5. DATA ANALYSIS 



120 




50 100 150 200 250 300 

E T [GeV] 



Figure 5.15: Plots of (pr) vs Et (left) and (y) 
(red). The lines pt = E^ and y = r\ are shown i 
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the geometric enhancement. The "central to peripheral" ratio, Rcp, is defined as 

-Cj 3-ji 
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Arccnt at ccnt ' 
7V coll iV evt 



N con 



pcriph a t pcriph 
iVovt, 



dp 3 



(5.2) 



The common phase space factors cancel, so the relevant quantity to extract from the data is 
the per event yield in pt bin i, 

Yi = ~ r^ 1 , (5.3) 



v evt dp l rj, 

in each centrality bin. The 60-80% centrality bin is used as the peripheral bin in all -Rcp an d the 
notation R^fi = / N®®^ 80 is used to describe overall geometric enhancement as determined 
by the Glauber Model. 

Jets from the 2010 Pb+Pb data set were reconstructed using the anti-fct algorithm with R = 



0.2 — 0.5 using the procedure described in Section 4.1 over the rapidity interval \y\ < 2.1. All 



jets failing the fake rejection criteria described in Sec. 4.3.3 are excluded from the analysis. These 



spectra are shown in Figs. 5.16 and 5.17 



This measurement is a sampling of the per event jet multiplicity distribution in a given p^ bin, 
which is given by the random variable and is sampled iV ev t times. The measured per event yield 
is the mean of this variable, 

1 1 



N, 



evt 



Yi 



1^ 71 j 



(5.4) 
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Figure 5.16: Per event yields before correction for R = 0.2 (left) and R = 0.3 (right) jets. 



A single event can produce multiple jets; this introduces a statistical correlation between bins 
and a full covariance matrix is required to describe the statistical uncertainty of the measurement 
described by 



where the quantity 



V;, 



Coy (Y^Yj) 



1 dN^ 



N, 



cvt 



1 



1 N cvt 



(5.5) 



(5.6) 



-^evt dp l T dp 3 T dp^dpj* ^cvt k= 
must also be recorded on a per event basis in addition to Yi. The measured covariances are shown 



(5.7) 



in Figs. 5.18 and 5.19, as the variance- normalized correlation p 

= CovjY^Yj) 
Pij ^V^^VaT^ ' 

The matrix is strongly diagonal, indicating a weak correlation between different bins. However, 
the unfolding procedure will correct for bin migration and will introduce additional correlations 
between bins, resulting in a non-trivial covariance matrix. 
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Figure 5.17: Per event yields before correction for R = 0.4 (left) and R = 0.5 (right) jets. 



5.4.2 Method 



The Monte Carlo is used to correct the raw spectrum for the finite jet energy resolution introduced 
by underlying event fluctuations and detector effects. For each truth jet, the nearest matching 
reconstructed jet with AR < 0.2 is found and the pair (p^ uth ,pJj? co ) is recorded to build the response 
matrix A. Then the input truth, Tj, and output reconstructed, Ri, distributions are related via a 
linear transformation 

i 

where Tj is the full truth pp spectrum, including truth jets that had no matching reconstructed 
jet. If the Monte Carlo accurately describes the underlying event and detector response, a similar 
bin migration is assumed to apply to the measured spectrum, hi as well: 



(5.9) 
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Figure 5.18: Spectrum correlations for EM+JES R = 0.4 jets. 



Xi is the desired spectrum corrected for bin migration. The goal of an unfolding procedure is to 
invert this equation in a controlled fashion and obtain the true spectrum. The full inverse of the 
response matrix may possess properties introducing undesired features on the data, in particular 
the inversion may amplify statistical fluctuations. Additionally, the system of equations is over- 
determined and modifications to the linear system, motivated by the specific application, must 
be applied to have a well-determined solution. This study addresses these problems by using a 
singular value decomposition (SVD) based approach to regularize the unfolding as presented in 
Ref. |22U| . The implementation of the SVD method in the RooUnfold software package [22 1 J was 
used to perform the unfolding [^J 

The first modification to Eq. |5.9| is to recast the system in terms of a linear least squares problem, 
where each equation in the system is inversely weighted by the uncertainty in the measured spectrum 



2 RooUnf old-1 .1.1 was used with minor additions to allow for checks with continuous values of the regularization 
parameter. 
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Figure 5.19: Spectrum correlations for EM+JES R = 0.2, 0.3, 0.4 and 0.5 jets for the 0-10% 
centrality bin. 



Abi 



. , . min . (5.10) 
4 ' Ah } 

* \ 3 

Thus equations with larger statistical uncertainty have less significance in the solution. In this 
framework the response matrix, A'^ = Aij/Abi, is expressed in terms of its SVD 



A' = USV T . 



(5.11) 



where X! is a diagonal matrix containing the singular values, Si, and U and V are matrices com- 
posed of left and right singular vectors which each form orthonormal bases. For small singular 
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values, 1/ Si can become large and statistical fluctuations in the data become magnified resulting in 
distortions in the unfolded distribution. To address this effect, the singular values are regularized 



by imposing smoothness requirements on the unfolded spectrum by expanding Eq. 5.10 to include 
the simultaneous minimization of the curvature of the spectrum, 

(Ax - b) T (A'x - b) + r(Cx) T (Cx) = min , (5.12) 

where bi = bi/ Abi, Cx is the second derivative of x expressed in terms of finite differences and r is 
a regularization parameter. The behavior for arbitrary r is obtained by rewriting Eq. |5.9| as 

A'C- 1 Cx = b, (5.13) 



and decomposing A'C as in Eq. 5.11 The form of C used in [220J and |221j assumes x is filled 
using bins of fixed width. In this analysis, variable bins are used so that lnp l T — Inp 1 ^ 1 is constant, 
and thus Cx is a finite difference approximation of d ^p T y ■ This measure of curvature is more 
appropriate given the steeply falling nature of the jet spectrum. 

The behavior of the system under statistical fluctuations can be understood by analyzing the 
behavior of the data expressed in the basis supplied by the SVD 

d = U T b. (5.14) 

The d{ are normalized by the error and statistical fluctuations are of order unity. The values of i 
for which \di\ 3> 1 are the statistically significant equations in the linear system, and the number of 
such equations is the rank of the system k. The remaining singular values correspond to statistically 
irrelevant equations and are regularized via 

— A -J*— , r = sl. (5.15) 

The final modification is to express the solution as a set of weights relative to another distribu- 
tion x ini , 

Xi = Wi/Xi ni i , (5.16) 
absorb the x; n i into the response, Aij = A^/ximj, solve for the weights 

(Aw - b) T (Aw - b) + r(Cw) T (Cw) = min , (5.17) 
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and multiply back by Xj n i to obtain x. 

For well-determined systems, these operations have no effect on the solution. However when 
the system is overdetermined this is no longer true, and X} n ; must be chosen carefully. In Ref. [220 
it is argued that the truth spectrum, T, used to generate A is the proper choice of Xi n i. 

In this study it was found that the turn on of the efficiency, which results from a combination of 
reconstruction efficiency and fake rejection, introduced additional curvature and thus larger values 
of Wi at low pt- This has the potential to propagate instabilities in the unfolding at low pr where 
the corrections are the least well known, throughout the entire spectrum, including regions where 
no efficiency corrections are required. Therefore the efficiency correction was removed from the 
unfolding problem by constructing A to not include inefficiency. Then the unfolded spectrum x is 
divided by e to correct for efficiency in the final step after the unfolding is complete. To construct 
3,'i n i the MC truth spectrum was re-weighted by a smooth power-law to remove any statistical 
fluctuations in the MC then multiplied by the efficiency, 

Xim(pr) = xoz(pt)Pt 5 ■ (5A8) 
5.4.3 Response Determination 

The SVD method relies on the curvature constraint to minimize the effects of statistical fluctuations 
on unfolding. This means that any non-smoothness in the response matrix or x; n i due to statistical 
fluctuations in the Monte Carlo sampling will distort the unfolding. To address this issue, the 
response matrices obtained from the Monte Carlo study were subjected to a series of smoothing 
procedures before begin used in the unfolding. 

For each J sample, a two-dimensional histogram, Mij, was populated from the (p^ nth ,p^ co ) 
of each truth-reconstructed matched jet pair. The p^ uth distribution, without requiring a recon- 
structed jet match, is was also recorded in a histogram Tj. The My and Tj from the different J 



samples are then combined using the cross section weighting discussed in Section 4.2. The response 
matrix was then constructed from these two distributions 

A l3 = , (5.19) 



which satisfies Eq. 5.8 if there is no inefficiency One-dimensional distributions, obtained by project- 



ing Mij in fixed p^ nth slices, are typically Gaussian-like with integrals equal to the total efficiency 
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-Nbins 


a 


Bin low 


Bin high 


pl£ w [GeV] 


p T igh [GeV] 


Truth 


35 


1.15 


5 


27 


12 


320 


Reco 


60 


1.09 


20 


44 


34 


270 



Table 5.3: Binning for truth and reconstructed spectra and response matrices showing the minimum 
and maximum pi[ ranges and the corresponding histogram bins. 



in that bin. The tails of these distributions are often produced by cases where a single truth jet 
is reconstructed as a pair of split jets and vice versa and can display behavior that is qualitatively 
different than the rest of the distribution. While this represents a small effect on the bin migra- 
tion, the changes in shape may cause problems with the additional smoothness constraints in the 
SVD method. Thus contributions to each slice less than 0.1% were removed and the rest of the 
distribution rescaled to keep the integral the same. The one dimensional slice of this distribution 



after smoothing for R = 0.4 jets in the 0-10% centrality bin are shown in Fig. 5.20 Checks were 
performed to test the sensitivity to this cutoff that showed no significant change in the final results 
with this cut removed. 

The binning used for the spectra and response matrices cover the range 7 < < 1000 GeV. 
The bin boundaries for the i th bin, p^. and p^-, were chosen such that P% i jp^ i = a is independent 
of i and fixed by the number of bins. The unfolding was performed on an interval, [Px w ,Px gh ]> °^ 
this range chosen based on requirements of the unfolding procedure. To allow proper bin migration 
in both directions the range of p<r values allowed for the unfolded spectrum must exceed the range 
of the reconstructed values on both sides (plj? w1:ruth < p!^ wreco and ^?!^ lstl truth > p T lghreco ). The 
unfolded result should then only be used over the reconstructed pt range, possibly excluding the 
most extremal bins on either edge of this range due to bias. The linear system is under-determined 
unless there are at least as many reconstructed bins as truth (unfolded) bins on their respective 



ranges. A bin selection satisfying these requirements is summarized in Table 5.3 



5.4.4 Statistical Uncertainties 

The measured covariance is used to generate statistical uncertainties in the unfolded spectrum 
by using a pseudo-experiment technique. A pseudo-experiment consists constructing a new input 
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Figure 5.20: ID projections at fixed p^ uth of the response matrix after smoothing for EM+JES 
R = 0.4 jets in the 0-10% centrality bin. A vertical line at 34 GeV is shown to indicate the minimum 
value of reconstructed used in the unfolding. 



spectrum, which is then unfolded. The procedure uses the measured covariance of the data to 
generate a new spectrum consistent with correlated Poisson fluctuations of the measured spectrum. 
1000 pseduo-experiments are performed, and the statistical covariance of the distribution over all 
such pseudo-experiments is a measure of the effect of the statistical uncertainty of the measured 
spectrum on the unfolded result. 

The effect of statistical fluctuations in the MC used to obtain the response matrix is computed 
in a similar fashion. In this case each pseudo-experiment uses a different response matrix, which 
is a possible statistical fluctuation of the input response matrix. These fluctuations are generated 
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by using the error on the matrix element to generate a new value of the matrix element consistent 
with a Poisson fluctuation of the original. The total statistical covariance is taken to be the sum of 
the covariance matrices obtained from these two procedures, which is equivalent to summing the 
uncertainties in quadrature. 

The unfolding method systematically associates features of the unfolded spectrum with statis- 
tical fluctuations and damps their contributions to the result through regularization. It is expected 
that increases in the regularization parameter may artificially reduce the sensitivity of the unfolding 



procedure to statistical fluctuations in the input spectrum. This effect is shown in Fig. 5.21, where 



the statistical errors on the Rcp have been evaluated using the pseudo-experiment technique for 
values of k from 4 to 10, with the statistical errors on the measured spectrum before unfolding 
shown for reference. The maximal value of k was chosen to be 10, which is approximately half the 
number of bins in the reconstructed sample with finite counts. This corresponds to the Nyquist 
frequency of the system: the highest k value that can be included without introducing aliasing 
effects in the signal decomposition and reconstruction. 
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Figure 5.21: Statistical errors on the unfolded Rqp evaluated using the pseudo-experiment method 
for 4 < k < 10 in various centrality bins for R = 0.4 jets. Statistical errors on the raw spectrum 
before unfolding are shown as a red curve. 
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5.4.5 d Distributions and Refolding 

The choice of the regularization parameter is driven by the distribution of d values as defined in 



Eq. 5.14 The coefficient d{ is the projection of the data (error normalized) onto the direction 



associated with the ith largest singular value. For values of \d{\ ^> 1, the corresponding equation 
in the now diagonal linear system is statistically well determined. Values of \di\ of order unity 
are consistent with statistical fluctuations in the data that select out a particular singular value. 
The contributions from these equations, which become large when the matrix is inverted, must be 
regulated. From these considerations, the value of k is chosen where the d distribution transitions 
from statistically relevant equations to equations consistent with statistical fluctuations. The d 



distributions for the R = 0.4 jets are shown in Fig. 5.22 k = 5 was chosen for all centrality bins 
and R values. 




2 4 6 8 10 12 14 16 18 2 4 6 8 10 12 14 16 18 2 4 6 8 10 12 14 16 1t 



Figure 5.22: d distributions for EM+JES R = 0.4 jets in various centrality bins, with regularized 
distributions overlaid with different values of k. 



Choosing k to be too small (or equivalently choosing r too large) will over-regularize the result 
by removing relevant information from the response matrix. Thus the choice of k must be made in 
conjunction with a refolding test, where the unfolded spectrum is multiplied by the full, unregulated 
response matrix and compared to the input spectrum. To make this comparison more rigorous both 
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in justifying the choice of k and providing a mechanism for analyzing the systematics, Eq. |5.15 



was 



extended to consider continuous r values. The refolded spectrum is compared to the original 
measured spectrum, and the degree of closure of the operations is quantified by a x 2 , 



,2 



X 



(Ax-b) T B _1 (Ax-b), (5.20) 



where B is the measured covariance of the data. Since x is the solution to the unfolding problem 
constructed from b, the fluctuations on b and the refolded spectrum Ax are not statistically 
independent. This means that this x 2 1S n °t a sum °f squares of independent variables with unit 
normal distribution, and does not follow a x 2 probability distribution. Despite this limitation, this 
variable, especially its behavior as a function of r, is still useful in describing the refolding closure 
so long as it is not misinterpreted. 



The distribution of the refolded x 2 as a function of r is shown in Fig. 5.23 The growth 
in x 2 a t large r in the neighborhood of the smallest k values shows that such regularizations 
remove enough information to prevent closure in the refolding. Small values of r maintain a 
faithful refolding procedure, with low x 2 > at the expense of amplifying statistical fluctuations in 
the unfolded spectrum. 

Under-regularization will allow statistical fluctuations to become amplified in the unfolding. 
Once an appropriate choice of Xj n i is made, the curvature of the weights distribution can be used 
to quantify the extent to which this effect occurs. There are two sources of curvature in the weight 
distribution: statistical fluctuations and differences between the unfolded spectrum and Xi n j. The 
optimal regularization point occurs when the statistical fluctuations have been suppressed, but 
differences between and the unfolded spectrum are not smoothed away. As the sources of 
curvature are changing near this value of r it is expected that the curvature distribution will change 
shape in this region. The curvature as a function of r is shown overlaid with the x 2 distributions 



in Fig. 5.23 Regions of flatness in the x 2 overlap in r with the qualitative changes in the curvature 
distribution. 

The procedure for extracting statistical uncertainties on the unfolded spectrum is discussed in 
Section |5.4.4| This addresses the question of how sensitive the unfolded spectrum is to statistical 
fluctuations on the measured spectrum for a given regularization procedure (i.e. choice of fcorr). It 
is possible that were a different statistical fluctuation of the data observed, a different regularization 
procedure could have been chosen, which is more sensitive to statistical fluctuations. This would 
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occur if different fluctuations had different \ 2 an d curvature distributions as functions of r. To 
check the stability of these distributions, the pseudo-experiment method was again used where 
these two distributions were generated for many possible fluctuations of the data. The stability of 
the two distributions under statistical fluctuations indicates that choice of regularization would not 



change under statistical fluctuations, and further justifies the procedure described in Section |5.4.4 
as a full measure of the statistical uncertainty. 
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Figure 5.23: x 2 an d curvature distributions as a function of the regularization parameter r for 
EM+JES R = 0.4 jets in various centrality bins. In the discrete case, values of r are fixed to be the 
singular values, = s\. The alternating colors indicate a range of r values between those allowed 
in the discrete case, with the r corresponding to a fixed k value occurring at the boundaries of the 
colored regions. 



5.5 Systematic Uncertainties 

Systematic uncertainty on the final, unfolded result is estimated both from ingredients from the 
MC into the unfolding and the unfolding procedure itself. For some sources, the uncertainties in 
Rqp are point-wise correlated, e.g. in i?cp vs - Pt changes in the JER cause all points in the curve 
to move up/down together. Such systematics are indicated by continuous lines drawn around the 



CHAPTER 5. DATA ANALYSIS 



133 



Systematic 


PT 


Centrality 


R 


Regularization 
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Efficiency 
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X 




JES 
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X 


X 


JER 


X 




X 


iVcoll 


X 


X 


X 



Table 5.4: Systematics on i?cp with point-wise correlations when plotted as a function of different 
variables are denoted by an x . 



values of the Rqp- Which systematics show point- wise correlations depend which variable the Rqp 



is being plotted against and are summarized in Table 5.4 As the Rcp is a ratio between two 
different centrality bins, errors that are correlated in centrality are taken from the variation of the 
Rcp due to that systematic. Errors on the -Rcp due to sources that are uncorrelated in centrality 
are taken as independent errors on the central and peripheral spectra and combined in quadrature. 
Errors that show no point-wise correlation are shown by shaded boxes and indicate the maximal 
variation of the points independent of one another due to such systematics. The relative error of 



each individual contribution to the systematic error on -Rcp vs Pt is shown in Figs. 5.24 and 5.25 



5.5.1 Regularization 

The systematic uncertainties associated with the unfolding procedure were estimated by considering 
the sensitivity of the procedure to the regularization parameter. The unfolding was performed with 
a different value of k and the variation in the unfolded spectra were taken as systematic errors. 
The errors on the central and peripheral spectra were then combined to determine an error on 
the -Rcp- Independent systematics were generated for both increased and decreased regularization 
by unfolding with k = 4 and k = 6. 

An alternative estimate of the systematic uncertainty associated with under-regularization was 
also tested and gave similar results as regularizing with the k = 4. This evaluation was motivated 
by linear least squares problems where the systematics are commonly estimated by considering 
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Figure 5.24: Contributions of various systematic uncertainties on Rqp as a function of jet for 
all centrality bins expressed as a relative error for the R = 0.2 (top) and R = 0.3 (bottom) jets. 
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Figure 5.25: Contributions of various systematic uncertainties on Rqp as a function of jet 
all centrality bins expressed as a relative error for the R = OA (top) and R = 0.5 (bottom). 
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variation of the x 2 by 1 about the minimum value 

A X 2 = xl-X 2 = X 2 /NDF. (5.21) 

The unfolding problem is not a strict minimization of the \ 2 as X 2 monotonically increases with 
r. Thus for the purposes of estimating the effect of increasing r, the system is effectively at a 
minimum and the above estimate was used to determine a r value. The spectra were unfolded with 
this t and the variation in the spectra were used as above to generate a systematic on the Rcp- 

5.5.2 x ini and Efficiency 

One advantage of the SVD procedure is that the unfolded result is less sensitive to the choice of 
a^ini than the truth spectrum used in a bin-by-bin unfolding procedure. To test this sensitivity the 



power in Eq. 5.18 was varied ±0.5. 

In peripheral collisions the fake jet rate is expected to be vanishingly small above 40 GeV. The 
fraction of jets remaining after fake rejection as a function of in peripheral data is consistent 
with the same quantity in all centralities in the MC as shown in Fig. |4.8| There is a small deviation 
at the highest pt where the MC does not fully saturate. However, this disagreement is present in 
all centrality bins and divides out in the Rqp- No systematic error was taken from this aspect of 
the efficiency. 

The uncertainty on the jet reconstruction efficiency prior to any fake rejection was estimated by 
comparing efficiencies in the HIJING sample and the overlay. The uncertainty in the ratio of central 
to peripheral efficiencies was taken conservatively to be 5% for < 100 GeV for all centralities 
and R values. 

5.5.3 Energy Scale 

Systematic uncertainties in the jet energy scale can arise both from disagreement between the MC, 
where the numerical inversion constants are derived, and data, as well as systematic shifts in the 
energy scale between different centralities associated with the underlying event subtraction. The 
former source only affects the absolute jet energy scale uncertainty which divides out in the Rqp- 
This is true to the extent to which the changes in absolute scale do not distort the unfolding 
differently in different centralities. 
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R 


- 10 % 


10 - 20 % 


20 - 30 % 


30 - 40 % 


40 - 50 % 


50 - 60 % 


0.2 


0.5 % 


0.5 % 


0.5 % 


0.5 % 


0.5 % 


0.5 % 


0.3 


1.0 % 


0.5 % 


0.5 % 


0.5 % 


0.5 % 


0.5 % 


0.4 


1.5 % 


1.0 % 


0.5 % 


0.5 % 


0.5 % 


0.5 % 


0.5 


2.5 % 


1.5 % 


1.0 % 


0.5 % 


0.5 % 


0.5 % 



Table 5.5: Energy scale shifts relative to peripheral used to generate systematic uncertainties. 



Relative changes in the energy scale between central and peripheral collisions are estimated by 



combining two studies. The difference in closure in the MC JES analysis, as presented in Sec. 4.4 



provides sensitivity to the effect of the background subtraction on the JES. Additionally, the in-situ 



JES validation comparing the track-calorimeter jet energy scales in Sec. 5.3.1.2 provides further 
constraints as it is sensitive to differences in jet fragmentation between the two centralities. 

Using these studies as input, systematic uncertainties were estimated on the energy scale for each 



centrality relative to the peripheral and are shown in Table 5.5 The MC samples were processed 
adding an additional shift, pt - > Pt + Apx, to the pt of each reconstructed jet when filling out the 
response matrices. A constant shift was applied above 70 GeV. Below the shift value was a linear 
function of pt, increasing to twice the high-px value at 40 GeV, 

Apr = [ / (l + Z2 3cF) Pt < 70 GeV ^ 
PT 1 / p T > 70 GeV 



where / is the value given in Table 5.5. The data was unfolded using the new response matrices 
resulting in a modified spectrum used to compute a new Rqp- Variations from this i?cp from the 
original in each px bin were taken as symmetric errors due to this systematic. 

5.5.4 Energy Resolution 

To account for systematic uncertainties coming from disagreement between the jet energy resolution 
in data and MC, the unfolding procedure was repeated with a modified response matrix. This 
response matrix was generated by repeating the MC study but with modifications to the Apx for 
each matched truth-reconstructed jet pair. A "detector" systematic was constructed to account 
for uncertainty between the data and MC without considering the effects of the underlying event. 
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Figure 5.26: N X M tower fluctuation distributions and fits. Values of iV and M associated with 
jet R values are indicated in the legend. 



This procedure follows the recommendation for JER uncertainty in 2010 jet analyses in pp. The 
JetEnergyResolutionProvider tool [222j was used to retrieve the fractional resolution, Cj^^ as 
a function of jet pt and rj |223j . The jet pr was then smeared by 

Pt - >■ Pt x 0jKa ) , (5.23) 

where M(0, 1) is the unit normal distribution. 

A separate systematic uncertainty considering the effect of disagreement of underlying event 



fluctuations in the data and MC was also determined. As discussed in Section 5.3.2.2, MC is 



internally consistent between the N x N tower RMS E^ distributions and the JER as a function 



of Et as expressed in Eq. 5.1 The b term is sensitive to fluctuations in the underlying event, while 



the a and c terms are independent of centrality. This is shown in Fig. 5.11, where fits of the JER 
with this form only differ significantly in the value of the b parameter. Additionally, the underlying 
event fluctuations, as measured by the RMS of the £/p distribution of N x N tower groups, is 
consistent with the b value obtained from the fit of the JER for jets with the same area. 

The MC fluctuations, which give b as a function of FCal E^ are fit with the functional form 



fe(^ Cal ) 



(5.24) 
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R 


N x M 


/[%] 


£§, [ GeV] 


a 


P 


7 [ GeV] 


0.2 


3x4 


2.5 


94.4 


0.622 


0.0083 


181.0 


0.3 


5x5 


2.5 


81.9 


0.641 


0.0133 


89.3 


0.4 


7x7 


5.0 


75.1 


0.680 


0.0168 


50.5 


0.5 


9x9 


7.5 


75.4 


0.712 


0.0176 


38.8 



Table 5.6: Fit parameters for background fluctuations in N x M towers as a function of ££^ Cal . 

and the values of the fit parameters E^, a, j3 and 7 for the different R 



which is shown in Fig. 



5.26 



values are presented in Table 5.6 along with maximal estimates of the MC/data deviation of these 
distributions given by the quantity /. 

As in the detector systematic, a new response matrix was constructed. The Apr was rescaled 
by an amount 



Apt — > Apt 



^JER(ft(l+fl)) 

ojerO) 



(5.25) 



b is determined from the £i?£ Cal and g is chosen to account for differences in the fluctuations 



between data and MC based on the values of / shown Table 5.6 In most central 10% of collisions 
the tower Ex distributions in the MC are broader than the data, so the effect of the fluctuations 
was symmetrically reduced by using g = —f. This systematic was applied as a one-sided error 
to the most central bin. Any other differences in the RMS results from the data having larger 
downward fluctuations, but smaller upward fluctuations. To account for this effect the g = 2/ was 
chosen for Apx < and g = —f for Apx > 0. This was applied as an asymmetric error to all but 
the most central bin. In no region of centrality does the data have larger upward fluctuations than 
the MC. 



5.5.5 N C( ,u Uncertainties 



The uncertainties in the iV co u values obtained from the centrality analysis are shown in Table 5.1 
The errors reported for all centrality bins are partially correlated because individual variations in 
the sensitive parameters in the Glauber calculation cause all N co \\ values to increase or decrease 
together. To evaluate the systematic uncertainties in the ratios of N co \\ between two centrality 



bins, iCn 



-^coffV-^coii 8 °' t na * appear in the expression for i?cp we have directly evaluated how 
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this ratio varies according to variations of the sensitive parameters in the Glauber Monte Carlo. 
We obtain from this study separate estimates for the uncertainty in the N co n ratios; the results are 



provided in Table 5.1 



5.5.6 Consistency Checks 

In addition to the systematic variations discussed previously, cross checks were performed by varying 
aspects of the analysis. These checks serve to test self-consistency of the procedure, not to quantify 
sensitivities to quantities expected to be uncertain. One such example is the fake rejection scheme, 
where it is natural to consider whether the default scheme imposes a fragmentation bias on the jets 
included in the spectrum and the Rqp- An alternative fake rejection scheme was assessed where 
the track jet matching requirement was replaced by a single track match within AR < 0.2 of the 
jet axis, where the track is subject to the following selection criteria: 

• pt > 4 GeV 

• do, zosin6 < 1.5 mm 

l _do_ go sin 9 < 5 q 

a d ' °"z sin B 

• Np ixe i > 2 

• N S ct > 7 

These requirements are significantly looser than those required on the tracks input into the track 
jet reconstruction. Although these looser matching criteria may lead to increased sensitivity to 
fake jets at low pp> as long as the fake rate is not dramatically enhanced, the lower p^ threshold 
(4 GeV vs 7 GeV for track jets), should address any concerns regarding a fragmentation bias. An 
identical analysis was performed using measured spectra and response matrices generated from 
this alternate scheme. The ratio of the Rcp values to those from the default scheme is shown in 



Fig. 5.27 with the statistical errors taken from errors generated on the nominal -Rcp through the 
pseudo-experiment method. The Rqp values show good agreement to the level of a few percent. 
The residual difference, an enhancement at low p^ that decreases with centrality, is consistent with 
the expected slight increase in fakes. 
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Figure 5.27: Ratio of Rcp values in the alternative fake rejection scheme (single track and cluster) 
to the nominal scheme (track jet and cluster) for R = 0.4 jets in various centrality bins. Statistical 
errors are taken from the denominator only. 



The GCW calibration scheme has a reduced sensitivity to jet width and flavor content relative 
to the EM+JES [216J. A comparison between these two calibration schemes was performed and 
the ratio of the Rqp obtained from the two schemes prior to unfolding is shown for R = 0.4 jets in 



Fig. 5.28 The distribution shows good agreement, indicating that the calorimeter response is not 
significantly altered for quenched jets. 

Consistency checks were performed to investigate any biases introduced by the unfolding pro- 
cedure or the choice of binning and pt range. The reconstructed spectrum from the MC sample 
was treated like the data and unfolded. The unfolded spectrum was then compared directly to the 
MC truth, and the ratio is taken to quantify the closure of the unfolding procedure. Since the final 
physics result is the Rqp, the bias on this quantity was also investigated. For each centrality bin 
the unfolded spectrum was divided by the unfolded spectrum in the 60-80% centrality bin to form 
an Rqp (with no N co \\ factor). The analogous quantity was constructed from the MC truth and 
the ratio of the i?cp' s was taken. Finally, a refolding test was preformed. When the unfolding was 
performed on the data, only this last check was possible. Any relationship between the unfolding 
non-closure in MC truth and MC refolding provides a context for interpreting the refolding in data 
as a potential bias on the unfolded data result. 
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Figure 5.28: The ratio of R C p for EM+ JES and GCW R = 0.4 jets before unfolding for data (black) 
and MC (red). Neither distribution shows a systematic difference between the two calibration 
schemes. 



The results of this test are shown in Fig. 5.29, regularized with k = 5, for the R = 0.4 jets. The 



MC truth non-closure shows some structure at the 5% level which is common to all centrality bins 
and R values (not shown). This structure is associated with the finite regularization parameter, and 
the effects can be removed by increasing k. Since the MC sample has much higher statistics than 
the data, the optimal value of k is likely to be much higher in the MC than in the data, and k can 
be increased without introducing strong effects from statistical fluctuations in the reconstructed 
spectrum. To facilitate the comparison with data, the value of k was chosen to be near the 
values appropriate for the data, and thus the MC unfolded spectrum is slightly over-regularized. 
Nevertheless the comparison shows good closure and insignificant systematic bias in the procedure. 
The i?cp shows even less bias as some of the structure is common to all centrality bins. The 
refolded reconstructed spectrum shows similar structure, suggesting that this check will be sensitive 
to residual bias in the data. 
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Figure 5.29: Results of the unfolding self-consistency check performed on the MC for R = 0.4 jets 
in different centrality bins. The ratio of unfolded to MC truth for both the spectrum and i?cp is 
shown in red and black respectively. The refolding on the reconstructed spectrum is shown in blue 
to facilitate a comparison with the analogous quantity in the data. 
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Chapter 6 

Results 



6.1 Jet Spectra and Rqp 



6.1.1 Unfolded Spectra 



The unfolded spectra are shown in Figs. 6.1 and 6.2 including the effect of the efficiency correction. 
The ratios of these spectra to the raw distribution as well as the refolding closure test are shown 



in Figs. 6.3 and 6.4 The unfolding procedure alters the measured statistical correlations between 



bins resulting in large off-diagonal covariances. These covariances were obtained using the method 



described in Section 5.4.4 and are shown in Figs. 6.5 and 6.6 
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Figure 6.1: Per event yields after unfolding and efficiency correction for R = 0.2 (left) and R 
(right) jets. 
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Figure 6.2: Per event yields after unfolding and efficiency correction for R = 0.4 (left) and R 
(right) jets. 
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Figure 6.3: Results of unfolding in centrality bins presented as various ratios to the raw spectrum for 
the R = 0.2 (top) and R = 0.3 (bottom) jets. The ratio of the unfolded spectrum both before and 
after efficiency correction to the input spectrum are shown in red and black respectively. The ratio 
of the refolded spectrum to the input distribution is shown in blue with the error bars indicating 
the relative error on the data. For all centrality bins this ratio is near unity, indicating good closure 
in the unfolding procedure. 
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Figure 6.4: Results of unfolding in centrality bins presented as various ratios to the raw spectrum for 
the R = 0.4 (top) and R = 0.5 (bottom) jets. The ratio of the unfolded spectrum both before and 
after efficiency correction to the input spectrum are shown in red and black respectively. The ratio 
of the refolded spectrum to the input distribution is shown in blue with the error bars indicating 
the relative error on the data. For all centrality bins this ratio is near unity, indicating good closure 
in the unfolding procedure. 
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Figure 6.6: Spectrum correlations after unfolding for R = 0.2, 0.3, 0.4 and 0.5 jets for the 0-10% 
centrality bin. 
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6.1.2 R C P 

The jet Rqp as a function of px are shown in Figs. 6.7 - 16. 10 with the black bands and shaded boxes 



representing the correlated and uncorrelated errors composed of the sources indicated in Table 5.4 



A subset of the results for the R = 0.2 and R = 0.4 jets are presented in a slightly different fashion 



in Figs. 6.11 and 6.12 In these figures the px-dependences of selected centrality bins is compared 
in different vertical panels. These results indicate that jets are suppressed in the most central 
collisions by approximately a factor of two relative to peripheral collisions. For all jet radii, the 
suppression shows at most a weak dependence on jet px- The correlated systematic error bands 
indicate that there may be a slight increase in Rqp in the range 40 < px < 100 GeV, although 
this is not present in all centralities and radii. For some bins, the .Rcp shows a mild decrease with 
increasing px at the highest px values, this effect is also marginal within the systematic errors, and 
will be discussed further in Chapter [7j 

The Rqp as a function of centrality, as expressed by the average number of participants in each 
centrality bin, is shown for different bins in jet px in Figs. |6.13 - 6,16| Note that the composition 
to the two types of systematic errors changes when plotting Rqp vs centrality as opposed to Rqp 
vs px- In general the iVp art dependence shows a weak dependence on jet px- The N paTt dependence 
for the lowest px bins, shows a linear drop in Rqp in the most peripheral collisions before flattening 
out. This dependence is slightly different than in the highest pt bins, where the dependence is 
more linear throughout the full N par t range. 

Finally the R dependence of the Rcp is evaluated directly. This comparison is shown for the 



0-10% centrality bin for different px bins in Fig. 6.17 A comparison of this dependence for different 



centrality bins but at fixed px, 89 < px < 103 GeV, is shown in Fig. 6.18 These results indicate 
an at most weak reduction in the suppression with increasing jet radius. There is essentially no 
reduction in suppression in the 50-60% bin, where there is a slight increase in the -Rcp with R in 
the most central bin. 

An alternative evaluation of the dependence of the Rcp on jet radius is provided in Fig. |6. 19 



which shows the ratio of Rqp values between R = 0.3, 0.4, and 0.5 jets and R = 0.2 jets, Rqp/Rqp, 
as a function of px for the 0-10% centrality bin. When evaluating the ratio, there is significant 
cancellation between the assumed fully correlated systematic errors due to i^jf , JES, JER, and 
efficiency. The results in the figure indicate a significant dependence of the Rqp on jet radius for 



CHAPTER 6. RESULTS 



151 



Pt iS 100 GeV in the most central collisions. For lower values, Rq P /Rqp for both R = 0.4 
and R = 0.5 differ from one beyond the statistical and systematic uncertainties. The maximal 
difference between the R = 0.2 and R = 0.5 jets occurs px ~ 60 GeV, where Rgp/R$$ exceeds 1.3. 
The central values for R = 0.3 jets also differ from one, but that difference is not significant when 
the systematic errors are accounted for. 
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Figure 6.7: Jet i?cp vs pt for R = 0.2 jets in different centrality bins. 
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Figure 6.8: Jet Rqp vs pt for R = 0.3 jets in different centrality bins. 
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Figure 6.9: Jet i?cp vs Pt f° r R = 0.4 jets in different centrality bins. 
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Figure 6.10: Jet i?cp vs pt for R = 0.5 jets in different centrality bins. 
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Figure 6.11: Jet i?cp vs |?t for R = 0.2 jets in different centrality bins. 
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Figure 6.12: Jet i?cp vs |?t for R = 0.4 jets in different centrality bins. 
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Figure 6.13: Jet Rqp vs -/V part for R = 0.2 jets in different bins. 
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Figure 6.14: Jet Rep vs -/V part for R = 0.3 jets in different px bins. 
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Figure 6.15: Jet Rqp vs -/V part for R = 0.4 jets in different bins. 




Figure 6.16: Jet Rqp vs -/V part for R = 0.5 jets in different px bins. 
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Figure 6.17: Jet i?cp vs -R f° r different bins in the 0-10% centrality bin. 



CHAPTER 6. RESULTS 



159 



Q_ 

O 



1.5 
1 

0.5 



Ok 



;89< p T < 103 GeV 


Pb+Pb V S A,A,= 
| L dt = 7 |ib" 1 


2.76 TeV 




: 50 - 60% 









1 

0.5 


1 

0.5 


1 

0.5 












♦ i : 


: 30 -40% 














— • — 




— ♦ — 
: 10-20% 


• 







: 0-10% 

















ft = 0.2 



R = 0.3 



ft = 0.4 



ft = 0.5 



Figure 6.18: Jet -Rcp vs R f° r different centrality bins for jets with 89 < < 103 GeV. 
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Figure 6.19: Ratios of i?cp values between R = 0.3,0.4, and 0.5 jets and R = 0.2 jets as a 
function of p^ in the 0-10% centrality bin. The error bars show statistical uncertainties obtained by 
propagating the statistical uncertainties on the individual Rcp values. The shaded boxes indicate 
partially correlated systematic errors. The solid lines indicate systematic errors that are fully 
correlated between different pt bins. 
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6.1.3 Energy Loss Estimate 



If a simple model of energy loss is assumed, the pp-dependence of the Rcp can be isolated. Following 
the discussion in Ref. |224j . an Rcp that is independent of pp can be the consequence of energy 
loss that is proportional to the jet's original px- The number of jets measured at a final pp after 
energy loss, p T , will be populated by jets that were produced with an original momentum of p T , 



P T = (1 - Sloss)p T 



(6.1) 



If the unmodified (i.e. initial) jet spectrum is described by a power law which is independent of 
centrality, 

dN A 

(6.2) 



dN A 
d-PT pp™ 

The final measured spectrum is modified by the transformation from p^ to p^: 



dN 
dp { T 



cent 



dN dpi? 
dp 1 ^ dp { T 



(6.3) 



cent 



If the energy loss is defined relative to peripheral collisions then the initial and measured spectra 
are equal, 

dN dN A 

(6.4) 



dpt. 



_ dN _ A 

periph <^Pp p^ 



In central collisions the measured spectrum is related to the initial by, 



dN 
dpL 



A 



cent 



(1 " SoaO-^W 



Rcp — 
dp { T 



(6.5) 



periph 



Therefore the Rcp and lost energy fraction, «Si oss , are related by 



Rcp = (1 - Si oss ) n 1 



'loss 



I -Rcp 



l/(n-l) 



(6.6) 



The centrality and R dependence of S\ oss have been extracted from the measurements made here 
of i?cp f° r Pt > 70 GeV, where the assumptions of power law spectrum and pp-independent Rqp 
are appropriate. The power, n, was obtained from fitting the pp spectra in the 60-80% centrality 



bin, and the extracted fit constants and errors are shown in Table 6.1 



In each centrality bin the Rcp was fit to a constant, Rqp, by minimizing the x function, 



X 



E E^Pi - Rcp^HRcpj - Rcp) 



(6.7) 
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R 


n ± 5 n 


k±5 k 


S±6 S 


2 


5.49 ± 0.07 


0.90 ± 0.10 


8.21e-04 ± 5.87e-05 


3 


5.64 ± 0.07 


0.88 ± 0.14 


9.13e-04 ± 4.92e-04 


4 


5.72 ± 0.06 


0.88 ± 0.20 


7.46e-04 ± 4.38e-04 


5 


5.78 ± 0.06 


0.74 ± 0.12 


1.25e-03 ± 7.81e-04 



Table 6.1: Values of fit parameters and associated errors for different jet radii. The n values were 
obtained from a power law fit of the jet spectrum in the 60-80% centrality bin. k and S were 
obtained from fitting the iV part dependence of S\ oss . 



Here Cij is the full pr covariance matrix of Rqp , which contains contributions from statistical corre- 
lations, uncorrelated systematic errors and correlated systematic errors. The statistical covariance 
was obtained in Chapter [6] For uncorrelated systematic errors, the covariance matrix was taken as 
diagonal Cfj = afa^5ij, and for fully correlated errors the covariance matrix is = crfa^. The 



solution to the x minimization, i?cp is expressed analytically in terms of the full covariance as 



RcP = ~y r j • ( 6 - 8 ) 



An uncertainty on this quantity was determined by varying -Rcp such that the \ 2 P er degree of 
freedom differed from its minimum value by one. This analysis was performed in each centrality bin 



for all R values with the resulting evaluation of S\ oss shown in Fig. 6.20 The grey boxes indicate 



the combined error due to uncertainties in both i?cp an d n with the horizontal size determined 



by the uncertainty on iVp ar t given in Table 5.1 the points have been artificially offset to facilitate 
comparison. The smaller radii show a nearly linear growth in S\ oss with iVp art . However the rate of 
growth diminishes with increasing R value indicating that some of the lost energy may be recovered 
with the larger jet definition. To provide a quantitative estimate of the iVp art dependence, the Si oss 
for each R value was fit with the functional form, 

Sioss = SN part k , (6.9) 

with the fit values of S and k and their associated fit errors shown in Table IBTTI and the fit functions 



shown with the points on Fig. 6.20 
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Figure 6.20: S\ function of -/Vp ar t for different jet radii, with fits using the form, S\ oss — 

SN patlt k • The grey boxes indicate the combined error due to uncertainties in both R C p and n with 
the horizontal size determined by the uncertainty on N paTt ; the points have been artificially offset 
to facilitate comparison. 
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6.2 Dijet Asymmetry 

The observation of a centrality-dependent dijet asymmetry was the first indication of jet quenching 
at the LHC and was reported by ATLAS in Ref. pp. That publication was a snapshot of an ongoing, 
more extensive asymmetry analysis which is one of the main results of this thesis. The results pre- 
sented here improve on the first published result by using the increased statistics of the full collision 
data set and benefiting from improvements in the jet reconstruction performance. These results 
were presented at Quark Matter 2011 |225l 1226] and are presented more thoroughly in Ref. [227J. 
It should be noted that the jet reconstruction procedure has evolved between the completion of the 
asymmetry analysis and the i?cp analysis. In fact, these improvements were a significant techni- 
cal undertaking designed to improve the performance and enhance the precision and pt reach of 
the Rcp measurement. The reconstruction procedure and performance as documented in previous 
chapters coincides with the Rqp analysis. While main features of the background subtraction (seed 
finding, flow correction and iteration step) were also present when the asymmetry analysis was 
performed, the details differ slightly (flow iteration step, E^ cut on seeds). Those details and the 
performance of that version of the reconstruction are documented Ref. [227] , and generally will be 
omitted here. One important distinction is the lack of jet energy scale calibration constants for 
R 7^ 0.4 jets at the time of the asymmetry analysis. An analysis of the asymmetry for R = 0.2 
jets was performed using the R = 0.4 calibration constants. The result is a ~ 5% non-closure in 



the JES, which drops to ~ 10% at the lowest E^s. This effect, shown in Fig. 6.21 for the 0-10% 
and 60-80% centrality bins, is independent of centrality and largely cancels in the asymmetry. The 
R = 0.4 jets do not suffer from this problem. 

Following the procedures used in Ref. [1] the per-event dijet asymmetries are calculated accord- 
ing to 

-C/T 1 + 2 

where En is the transverse energy of the leading (highest Et) jet in the event and £t2 is the trans- 
verse energy of the highest £/p jet in the opposite hemisphere (A<p > ir/2 with Acf) the azimuthal 
angle difference between the jets). Both jets are required to satisfy the \rj\ < 2.8 requirement. 
Events for which the second highest E^ jet in the event fails the r] or Acf) requirements are not 
excluded, but rather the highest jet that does satisfy these criteria is used. All presented results 
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Figure 6.21: The JES performance for R = 0.2 jets in the 0-10% (left) and 60-80% (right) centrality 
bins with the older version of jet reconstruction that was used in asymmetry analysis. As no R = 0.2 
calibration constants were available at the time, the R = 0.4 constants were used, resulting in a 
5 — 10% non-closure in the energy scale. 



have the requirement E^2 > 25 GeV; events for which the second jet (passing the above selections) 
fails this requirement are excluded from the analysis. 



Figure 6.22 shows the asymmetry distribution for R = 0.4 jets with Eri > 100 GeV for six bins 
of collision centrality. Also shown for each centrality bin are results obtained from pp measurements 
at 7 TeV [228J and the asymmetry distributions obtained from the HIJING+PYTHIA MC samples 
for the corresponding centrality bin using the above described analysis procedures. The features 



observed in the original dijet asymmetry publication are seen in Fig. 6.22, namely good agreement 
between Pb+Pb data, pp data, and MC results in the more peripheral (40-60% and 60-80%) 
centrality bins and an increasingly strong modification of the asymmetry distributions for more 
central collisions. The MC Aj distributions show modest broadening from peripheral to central 
collisions due to the increased underlying event fluctuations. However, the modifications seen in 
the data are much stronger than those seen in the MC for which the underlying event fluctuations 



are shown to be consistent with those in Pb+Pb data in Fig. 5.9 



CHAPTER 6. RESULTS 



166 



% , Centrality 0-10% y_ Centrality 10-20% jl Centrality 20-30% j 



4 

"D 3.5 



4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 



ATLAS Preliminary 
Pb+Pb \[s^ = 2.76 Te\T" 

L int = 7^b-1 




Centrality 30-40% 



E T1 > 100 GeV 
E T2 >25 GeV 
R = 0.4 



Centrality 40-60% 



Centrality 60-80% 




~i 1 i 1 i 1 r 



□ hijing+pythia. 

• Pb+Pb Data 
o p+p Data 




0.2 0.4 0.6 0.1 



0.2 0.4 0.6 0.8 



0.2 0.4 0.6 0.8 



Figure 6.22: Dijet asymmetry for R = 0.4 jets in six centrality bins in events with a leading jet 
with Ei > 100 GeV. A comparison to HIJING with embedded PYTHIA dijet events (yellow) and 
ATLAS yfs = 7 TeV pp data (open circles) is shown. 



Figure 6.23 shows the measured asymmetry distributions for R = 0.2 jets with E^i > 100 GeV 
compared to results from MC simulations. An analysis of pp data with R = 0.2 is not yet available. 
The R = 0.2 asymmetry distributions show the same general features as seen in the R = 0.4 results 



shown in Fig. 6.22 However, the asymmetry distributions generally extend to slightly larger Aj 
values for R = 0.2 jets than for R = 0.4 jets. Another difference between the R = 0.4 and R = 0.2 
results can be seen in the 40-60% centrality bin where a modification of the asymmetry distribution 
is seen for R = 0.2 jets and not for R = 0.4 jets. 

One of the most important results presented in the original dijet asymmetry letter was the 
observation that the dijet Acj) distribution remains mostly unchanged in all centrality bins while 
the asymmetry distribution is strongly modified in central collisions. Figures [6~24] and |6~25] show 
updated results for the R = 0.4 Acp distributions for dijets with En > 100 GeV using the same 
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Figure 6.23: Dijet asymmetry for R = 0.2 jets in six centrality bins in events with a leading jet 
with Ei > 100 GeV. A comparison to HIJING with embedded PYTHIA dijet events (yellow) is 
shown. 



dijet selection procedure described above. Both linear and logarithmic vertical scales are shown. 
These figures confirm the original observation that the dijet A(f> distributions are unmodified with 
the caveat that a small combinatoric contribution is observed in the 0-10% and 10-20% centrality 
bins. Such a combinatoric contribution was also seen in the original measurement. The dijet Acfi 



distributions for R = 0.2 jets with i?T2 > 100 GeV are shown in Figs. 6.26 and 6.27 with linear and 
logarithmic vertical scales, respectively. The R = 0.2 A(f> distributions show no modification and a 
much smaller combinatoric contribution. 

The evolution of the asymmetry with the -E/ri threshold is shown in Figures |6.28| and 6.29 



which contain distributions of dijet Aj for three different centrality bins, 0-10% (top), 30-40% 
(middle), 60-80% (bottom), for three different ranges of energies (see figures) for the leading jet. 
The peaking at larger values of Aj becomes more pronounced when E^i is restricted to a lower 
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Figure 6.24: Dijet Acj) distributions for R = 0.4 jets in ^/snn = 2.76 TeV Pb+Pb collisions having 
a leading jet with > 100 GeV in six centrality bins. A comparison to HIJING with embedded 
PYTHIA dijet events (yellow) and ATLAS yfs = 7 TeV pp data (open circles) is shown. 



range and that peaking becomes less pronounced as E^i increases. The peaking is particularly 
prominent for R = 0.2 jets with 75 < E^i < 100 GeV where the peak at Aj ~ 0.55 corresponds 
to -&T2 ~ 25 GeV - the minimum value for £'t2 allowed in the analysis. With increasing Eti 
the peak occurs at approximately the same E^2 value, but larger Aj. The Aj range accessible 
in the measurement is limited by the presence of low E^ real or false soft jets in the event that 
are selected as the second jet and by inefficiencies in the jet reconstruction at low Er- Because of 
the smaller jet size, the R = 0.2 jets have better efficiency at low E^ and fewer fake jets. These 
arguments would explain why the asymmetry distribution extends to larger values of Aj in the 



100 < Eti < 125 GeV bin. To complete the survey of dijet asymmetry measurements, Fig. 6.30 
shows the full centrality dependence of the Aj distributions for R = 0.2 jets with E^ > 75 GeV. 
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Figure 6.25: Dijet distributions plotted with a logarithmic vertical scale for R = 0.4 jets in 
y / SNN = 2.76 TeV Pb+Pb collisions having a leading jet with > 100 GeV in six centrality bins. 
A comparison to HIJING with embedded PYTHIA dijet events (yellow) and ATLAS yfs = 7 TeV pp 
data (open circles) is shown. 
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Figure 6.26: Dijet Acj) distributions for R = 0.2 jets in ^/snn = 2.76 TeV Pb+Pb collisions having 
a leading jet with > 100 GeV in six centrality bins. A comparison to HIJING with embedded 
PYTHIA dijet events (yellow) is shown. 
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Figure 6.27: Dijet distributions plotted with a logarithmic vertical scale for R = 0.2 jets in 
- v /snn = 2.76 TeV Pb+Pb collisions having a leading jet with > 100 GeV in six centrality bins. 
A comparison to HIJING with embedded PYTHIA dijet events (yellow) is shown. 
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Figure 6.28: Dijet asymmetry for R = 0.4 jets in three centrality bins: 0-10% (top), 30-40% 
(middle), 60-80% (bottom) for -^/snn = 2.76 TeV Pb+Pb collisions with a leading jet with 
100 < Eti < 125 GeV (left), 125 < E Tl < 150 GeV (middle), 150 < E T1 < 200 GeV (right). 
For all of the plots, comparison to HIJING events with embedded PYTHIA dijets for the same 
conditions on reconstructed jets as the data is shown (yellow). 
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Figure 6.29: Dijet asymmetry for R = 0.2 jets in three centrality bins: 0-10% (top), 30-40% 
(middle), 60-80% (bottom) for ^/snn = 2.76 TeV Pb+Pb collisions with a leading jet with 
75 < E T1 < 100 GeV (left), 100 < E T1 < 125 GeV (middle), 125 < E T1 < 150 GeV (right). 
For all of the plots, comparison to HIJING events with embedded PYTHIA dijets for the same 
conditions on reconstructed jets as the data is shown (yellow). 
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Figure 6.30: Dijet asymmetry for R = 0.2 jets in six centrality bins in events with a leading jet with 
£"t > 75 GeV. A comparison to HIJING with embedded PYTHIA dijet events (yellow) is shown. 
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Chapter 7 



Conclusion 



This thesis has presented results for the jet i?cp an d dijet asymmetry. The asymmetry distributions 
indicate that there are a class of events where there is a significant momentum imbalance between 
the two leading jets in the event. The extent of this imbalance gradually increases with centrality 
implying that the momentum balance between dijets is modified in the presence of the medium. This 
is strongly suggestive of the interpretation that the asymmetry is directly probing the differential 
energy loss of two partons traversing different in-medium path lengths. Furthermore, these dijets 
show an angular correlation that is consistent with vacuum dijet production. When this result 
was first released, the lack of modification of this angular correlation was considered surprising, 
especially given the copious energy loss implied by the asymmetry. 

The -Rcp results presented in this thesis show that jets are suppressed in central collisions 
by approximately a factor of two relative to peripheral collisions. This suppression is almost 
completely independent of ■ The approximate value of Rqp ~ \ and the slight increase with 
are qualitatively consistent with predictions using the PYQUEN MC event generator [229]. As the 
angular size of the jets is increased, by increasing the R parameter, the Rcp also increases. This 
i?-dependence of Rqp is weak and is a marginal effect above the measurement uncertainties. For 
lower energy jets, 50 < pr i$ 80 GeV, an increase in Rqp with R that is greater than that for higher 



energy jets was predicted in Ref. |183j . shown in Fig. 2.29 While the measurement uncertainties 
limit a strong conclusion, this trend is present in the data as well. 

The px-dependence can be compared with that of the single particle .Raa measured at the 



LHC shown in the summary plot in Fig. 7.1. The most recent RHIC results showed a constant 
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Figure 7.1: Single hadron Raa from various experiments combined with multiple theoretical calcu- 
lations. The LHC results extending to high p^ show a clear reduction in the suppression that was 
discovered at RHIC. Figure adapted from Ref. [230J. 

suppression Raa ~ 0.1 at the highest available pt- The LHC results have extended to higher px 
where the -Raa is found to rise up to ~ 0.5, which is in the same regime as the jet Rcp shown here. 
While it has been argued in this thesis that there is no requirement that the single particle and 
jet suppression factors be related, the plateau in the px-dependence exhibited by both observables 
suggests some common feature of the high energy loss regime. 

Since the anti-fct jets have highly regular geometry, the R parameter can be safely interpreted 
as a cone radius. Changing the jet radius has the effect that medium induced radiation emitted at 
angles between the old and new radii will be recovered by the jet and the energy loss will be reduced. 
An additional effect is the change in the spectral shape due to the change in definition, which is 
present even in the case of no quenching. Therefore, evaluations of the energy loss as a function of jet 
radius must be convolved with this second effect to be comparable to the experimentally observed 
spectra or i?cp- This is normally the case in theoretical evaluations of the Rqp or -Raa where the 
per parton energy loss must be convolved with the production spectrum. Equivalently, without 
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the analogous deconvolution, a comparison between the Rcp at two radii (i.e the experimental 
observation in this thesis) is indirectly sensitive to the energy recovery. 

A quantitative extraction of the energy loss was attempted in Ref. [224] with single particle -Raa 
measurements from PHENIX. The procedure used the unmodified inclusive single particle spectrum 
to remove the dependence on the spectral shape under the ansatz of an energy loss scenario where 
Aj*t oc pt- In the previous chapter, this procedure was applied to fully reconstructed jets with 
Pt > 70 GeV, with the extracted S\ oss values providing a basis of comparison between the different 



R values. The results of this exercise, summarized in Fig. 6.20 indicate an R dependence of the 
energy loss that is significant beyond the estimated uncertainty. At fixed A r par t, the differences in 
<Sioss between the two R values can be interpreted as the fractional recovered energy in the annular 
region between the two radii. The R = 0.2 and R = 0.3 values remain relatively close, however 
increasing the radius to R = 0.4 shows a reduction in the energy loss fraction. The marginal 
decrease in S\ oss is even greater when moving from R = 0.4 to R = 0.5. This trend indicates that 
the energy recovery grows faster than linearly with R. The relative magnitudes of S\ oss evolve with 
iVpart, indicating that the energy recovered by expanding the jet radius increases with S\ oss . 

To establish this behavior quantitatively, the iV par t dependence of S\ oss was fit with the form, 



•Sioss = SN pSLj: t k , with the fit results given in Table 6.1 This can be compared with the single hadron 
analysis performed by PHENIX. This analysis extracted larger values for S\ oss than those obtained 
here, which is consistent with the general expectation that some of the energy radiated by the 
parton may be recovered when using jets. Furthermore, the S\ oss was found to scale approximately 
as A'part 2 / 3 , which is consistent with the 1? energy loss models discussed in Section 



2.5.4 



The k 



values calculated here for full jets are much larger, with the R =0.2, 0.3 and 0.4 jets having k ~ 0.9. 
The R = 0.5 shows a slower iV par t 3 / 4 growth. 

The i?cp is only sensitive to nuclear modification effects relative to a peripheral reference 
spectrum in this case the 60-80% centrality bin. To quantify the full extent of these effects a pp 
reference spectrum should be used. The LHC recorded pp data at 2.76 TeV during early 2011, 
and a full Raa measurement using this data is the obvious next step in a systematic study of 
jet quenching at the LHC. A jet triggered sample of approximately 150 nb -1 was recorded by 
ATLAS, which provides ample statistics for the jet spectrum measurement out to pt — 250 GeV. 
As this reference spectrum will contain better statistics than the 60-80% Pb+Pb spectrum, the Raa 
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measurement can be performed in bins of rapidity as well. Since the Raa is sensitive to the full 
extent of the quenching, this measurement will likely focus on the maximal suppression observed 
in the 0-1% or 0-5% centrality bins. 

To more firmly associate the jet suppression observation with the quenching interpretation, 
suppression from the NPDF must be constrained through future measurements. The jet i?cp' s tend 
to exhibit a decrease with increasing p^, particularly in the most central bins. While this effect 
is marginally significant beyond systematic errors, a possible explanation could be the sampling of 
the NPDFs in the EMC region. This scenario must be investigated; any detailed measurements of 
the quenching require the decomposition of the total observed suppression into contributions from 
cold nuclear matter effects and a quenching component. This can be checked with a measurement 
of the direct photon Raa, which is sensitive to hard scattering production rates determined by the 
NPDFs, but not to quenching. Additionally the rates for Z production can be used to constrain cold 
nuclear matter effects, as the Z — > £1 decay modes should also be immune to the quenching. Both 
ATLAS [231J and CMS [232J have published Z i?cp values, however these results are statistics 
limited and cannot provide rigorous constraints on the NPDFs. This could be improved with 
results from the higher statistics 2011 run. The upcoming p+Pb run at the LHC should provide an 
opportunity to perform full jet reconstruction in these collisions for the first time. ATLAS should 
not only be able to provide precision measurements for determination of the NPFs, but also measure 



the impact parameter-dependence of the nuclear modification effects discussed in Section 2.5.2 

In the case of photon identification, a significant background comes from 7r°'s in jets. One of 
the primary techniques at rejecting this background in analyzing pp collisions is to require that the 
photon be isolated; this is typically accomplished by placing a maximum threshold on the energy 
deposited in a cone of fixed radius about the reconstructed photon axis, excluding the energy of 
the photon. The subtraction procedures developed for jet reconstruction in heavy ion collisions 
discussed in this thesis have been adapted by ATLAS to be applied to the isolation cones in the 
photon reconstruction. 

As discussed previously, dijet observables provide complementary information to the single 
inclusive measurements as they are sensitive to the relative energy loss of the two jets. Furthermore, 
both gamma-jet and Z-jet correlations can be used in analogy with the asymmetry to determine 
single jet energy loss, and are considered topics of future interest. The asymmetry results presented 
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here show that for systems with a 100 GeV leading jet, a significant fraction of events were found 
to have the sub-leading jet correlated in Acf) ~ tt and have -E/p2 ~ 25 GeV. This means that the 
quenching is strong enough to regularly cause energy loss of 75% of the jet's energy. Furthermore, 
the relatively strong angular correlation indicates that this energy loss must be incurred without 
significantly altering the direction of the jet. Whether this feature is captured by various energy 
loss models is not clear, as most models used the soft eikonal limit, which does not fully describe 
the transverse deflection of the leading parton. 

The measured angular correlation can be used to test predictions of modified acoplanarity 
distributions through elastic scattering [147- 149J . By rescaling the i^-dependence by the leading 



jet pt, this quantity, defined in Eq. |2.77[ can be related to the normalized dijet angular correlation, 

1 dN 1 da _ dP dP 
N~d0 K (7o(pr)# ~ PT d^ ~ d(K v / PT ) ' 
This quantity shows almost no modification in heavy ion collisions, and significantly less modifica- 
tion than the predicted acoplanarity distributions, which generally show a reduction of 25-50% in 
the maximum of the correlation relative to pp. 

The measured angular correlation can be used to constrain the magnitude of transverse mo- 
mentum kicks, (fc t 2 ), imparted to the jet through random elastic scattering in the medium. For a 
dijet system with back-to-back jets with momentum px before energy loss, each jet can undergo 
random small angle multiple scattering resulting in an RMS angular deflection of 



(h 2 



a jet ^<A0/ ct )^^— . ( 7 - 2 ) 
The angular deflection of each jet is independent of the other so the effect on the dijet distribution 
is an additional kick of Odijet = v2<7jet- An upper limit on <7j e t can be placed by sampling the 
dijet Acf) distribution obtained from the MC including additional Gaussian smearing with different 



values of Odijet and testing to see which values are consistent with the data. Figure 7.2 shows the 
data along with the nominal MC distribution and several smeared variants of the MC distribution 
for R = 0.4 jets with E^i > 100 GeV in the 0-10% centrality bin. The distributions are scaled to 
have unit value in the highest A</>dij e t bin to facilitate comparisons of the shapes. The Oj e t = 0.015 
distribution shown in green is clearly broader than the data, as is the <7j e t = 0.01 distribution shown 
in blue except in one bin. If the former is taken as an estimate of the upper bound on cTjet and 
100 GeV is used as an estimate of the jet £t, this corresponds to a limit of y (k t 2 ) < 1.5 GeV. 
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Figure 7.2: The A^ijet distribution for data (black), MC (filled yellow) and MC with additional 
Gaussian smearing using three different a values. The distributions are scaled to have unit value 
in the highest A^^ijet Dm - The errors represent the statistical errors in the data or MC sampling. 



The fact that the jet suffers essentially no angular deflection means that there is a single 
well-defined jet axis. For events with incomplete nuclear overlap the collision zone is oblong and 
there are a variety of path lengths the jet could take through the same medium. Furthermore, if 
the overlap is assumed to have elliptical geometry, the orientation and eccentricity of the overlap 
region are accessible on a per event basis through the flow observables ^ and V2 respectively. If 
coordinates are chosen such that the impact parameter vector is aligned with the y axis, the jet's 
linear trajectory is given by y = cot Aipx + y' — cotAipx', where Aip is the angle between the 
jet axis and the event plane, with the primed coordinates indicating the jet production location. 



The geometry of this scenario is shown in Fig. 7.3 The production point of the jet is unknown 



on a per event basis, but the average of such positions could be taken over many collisions with 
the same flow characteristics. In this geometrical picture, two-jet differential and single inclusive 
jet observables are sensitive to different moment-like quantities in the energy loss. Consider the 
variable £, which is some quantity sensitive to the energy loss. In a given event this should be a 
function of the production position, initial parton kinematics and collision geometry. The average 
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Figure 7.3: Elliptical geometry of collision zone, with the impact parameter direction chosen to 
align with the y-axis. The jet production point is located by the vector x'. Jets travel path lengths 
L\ and L<i through the medium which are denoted by blue and green lines respectively. The jet 
axis is oriented with angle Atp with respect to \&2- 

of this variable over events of similar eccentricity and centrality should be related to the differential 
energy loss, integrated over all different jet production points inside an ellipse by the eccentricity, 

<0 e (AVO = j ciVp(x')£(x', s, Aip) , (7.3) 

where p is the transverse density of possible jet production points. The dependence of £ on x' comes 
from the path length travelled in the medium, i.e. £(x', e, A^) = £(L(x', e, Aip)). The function L 
can be found by determining the coordinates of the jet-ellipse intersection and computing the path 
length. The coordinates of where the jet leaves the medium are defined by the solutions to the 
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simultaneous set of equations for the jet axis and the ellipse, 

x 2 + T ^ = R 2 A , (7.4) 

where e is the eccentricity and Ra is the nuclear radius which forms the semi-major axis of the 
ellipse. The eccentricity is assumed to be determined by centrality, but it may be possible to use 
per event v% values as an additional handle to constrain the range of eccentricity values in a class 
of events. Using the jet reconstruction techniques described in Chapter [4] the Er-mtegrated flow 
harmonics and total can be measured by excluding energy deposition associated with jets on a 
per event basis. This should allow better understanding of non-flow contributions to the harmonics 
as well as the fraction of total E^ and particle production that comes from hard processes. 

When considering the single inclusive jet energy loss, £ is sensitive to L directly. For example, 
in the BDPMS-Z scheme the energy loss for a single jet is approximately AE^ oc L 2 , using this 
variable for £, 

(A£ T } £ (AV0 oc J d 2 x'p(x')L 2 (x, e, A^) = 1?{A^, e) . (7.5) 

This is contrasted with the case where the dijet differential energy is the energy loss observable, 
A£/xi — A&T2 oc AL 2 . This quantity has a different dependence on the collision geometry since 
the path lengths of the two jets are correlated on a per event basis. Thus when averaging over 
similar events the quantity, 

(AE T1 - A£ T 2) e (AV>) oc J d 2 x'p(x')AL 2 (x,e,A^) = AZ^A^e), (7.6) 

provides an independent handle on the geometrical dependence of the energy loss. If put on a more 
rigorous theoretical footing, calculations of these two variables in an energy loss framework would 
provide a testable prediction of the path length dependence. The three-jet rate is high enough at 
the LHC that a similar approach could be applied to events with three jets. These have additional 
per event geometric correlations and would provide independent moments to check the energy loss. 

Relating a measured asymmetry to the partonic level energy loss can be challenging due to 
the intrinsic width of the vacuum asymmetry distribution. Furthermore, the distinction between 
multi-jet events and splitting due to hard radiation is not always clear. Some of these issues may 
be resolved by adding an additional step to the experimental analysis where the final jets are used 
as input to another round of clustering with a jet finding algorithm. For example, R = 0.2 jets 
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are robust against underlying event fluctuations, but are prone to more splitting than larger jet 
definitions. A set of R = 0.2 jets could be used as input to a less geometrically restrictive finding 
such as the k t or Cambridge- Aachen algorithm with larger R values, to build super-jet objects to 
be used in asymmetry-like analyses. 

Jet structure and properties should also provide details on the quenching mechanism. Measure- 
ments of jet shapes as well as the transverse and longitudinal distributions of particles within jets 
are sensitive to the angular pattern of the medium-induced radiation. A preliminary analysis of the 
jet fragmentation has been presented by ATLAS [227J, and a more complete study using the 2011 
data set is presently underway. This analysis should provide constraints on quenching mechanism 



as it relates to the medium modified fragmentation functions discussed in Section 2.5.4 Although 
experimentally challenging, the analysis of jet fragmentation can be applied to reconstructed pho- 
tons. The fragmentation photons provide an even cleaner probe of the angular radiation pattern; 
unlike the hadronic fragments the photons do not continue to interact with the medium after being 
emitted. 

It is hoped that the measurements presented here will inspire theoretical calculations in the 



context of the various paradigms discussed in Section |2.5.4| to consider the angular dependence 
of the energy loss and its effect on Rqp- Comparisons can be made both in terms of radiation 
recovery between the different radii as well as the total amount of large angle radiation (R > 5) to 
see whether this can give the expected total suppression of Rqp ~ 0.5. Any plausible quenching 
scenario must simultaneously satisfy both of these observations. 
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